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Abstract 

We analyze the one-dimensional periodic Kardar-Parisi-Zhang equation in the 
language of paracontrolled distributions, giving an alternative viewpoint on the sem¬ 
inal results of Hairer. 

Apart from a basic existence and uniqueness result for paracontrolled solutions 
to the KPZ equation we perform a thorough study of some related problems. We 
rigorously prove the links between KPZ equation, stochastic Burgers equation, and 
(linear) stochastic heat equation and also the existence of solutions starting from 
quite irregular initial conditions. We also build a partial link between energy solu¬ 
tions as introduced by Gongalves and Jara and paracontrolled solutions. 

Interpreting the KPZ equation as the value function of an optimal control prob¬ 
lem, we give a pathwise proof for the global existence of solutions and thus for the 
strict positivity of solutions to the stochastic heat equation. 

Moreover, we study Sasamoto-Spohn type discretizations of the stochastic Burg¬ 
ers equation and show that their limit solves the continuous Burgers equation pos¬ 
sibly with an additional linear transport term. As an application, we give a proof 
of the invariance of the white noise for the stochastic Burgers equation which does 
not rely on the Cole-Hopf transform. 
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1 Introduction 

The Kardar-Parisi-Zhang (KPZ) equation is the stochastic partial differential equation 
(SPDE) 

2z?/r(t, x) = x)) 2 + £(t, x), x € R, t ^ 0, (1) 
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where h : M + x R —> R is a real valued function on the real line, 2z? = dt — A denotes the 
heat operator, D = d/dx and dt = d/d t are the spatial respectively temporal derivatives 
and £ is a space-time white noise: the centered Gaussian space-time random distribution 
with covariance 

HS,(t,x)£(s,y)] = 5(t, - s)5(x - y), t,s^ 0, x,yeR. 

The KPZ equation was introduced by Kardar, Parisi, and Zhang [KPZ86] as an SPDE 
model describing the large scale fluctuations of a growing interface represented by a 
height field, h. Based on non-rigorous perturbative renormalization group arguments 
they predicted that upon a particular “1-2-3” rescaling and centering, the height field 
h (or at least its finite dimensional distributions) must converge to a scale invariant 
random field i?k pz (the KPZ fixed point ) obtained as 

-f^kpz (t,x) = lim A/r(A 3 i, X 2 x) — c(X)t. (2) 

A—>-oo 

According to the general renormalization group (RG) understanding of dynamic critical 
phenomena a large class of similar interface growth mechanisms must feature the same 
large scale statistical behavior, namely their height fields h(t, x) converge, upon rescaling 
(and maybe centering and a suitable Galilean transformation to get rid of a uniform 
spatial drift) to the same KPZ fixed point Pfkpz- 

Proving any kind of rigorous results about the convergence stated in (2) (called 
sometimes the strong KPZ universality conjecture ) is a wide open problem for any con¬ 
tinuous SPDE modeling interface growth and in particular for the KPZ equation; see 
however [SS10, ACQ11] for recent breakthroughs. On the other side there has been 
a tremendous amount of progress in understanding the large scale behavior of certain 
discrete probabilistic models of growth phenomena belonging to the same universality 
class, mainly thanks to a special feature of these models called stochastic integrability. 
For further details see the excellent recent surveys [Corl2, Qual4, QS15]. 

A weaker form of universality comes from looking at the KPZ equation as a meso¬ 
scopic model of a special class of growth phenomena. Indeed by the same theoretical 
physics RG picture it is expected that the KPZ equation is itself a universal description 
of the fluctuations of weakly asymmetric growth models, a prediction that is commonly 
referred to as the weak KPZ universality conjecture. In this case the microscopic model 
possesses a parameter which drives the system out of equilibrium and controls the speed 
of growth. It is then possible to keep the nonlinear contributions of the same size as the 
diffusion and noise effects by tuning the parameter and at the same time rescaling the 
height field diffusively. It is then expected that the random field so obtained satisfies 
the KPZ equation (1). 

For a long time, the main difficulty in addressing these beautiful problems and ob¬ 
taining any kind of invariance principle for space-time random fields describing random 
growth has been the elusiveness of the KPZ equation itself. The main difficulty in making 
sense of the equation stems from the fact that for any fixed time tfi 0 the dependence of 
the solution h(t,x ) on the space variable x € T cannot be better than Brownian. That 
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is, if we measure spatial regularity in the scale of Holder spaces < ^ 7 (M), we expect that 
h(t, •) € < ^ 7 (M) for any 7 < 1/2 but not better, in particular certainly the quadratic 
term ( Dh(t,x )) 2 is not well defined. 

The first rigorous results about KPZ is due to Bertini and Giacomin [BG97], who 
show that the height function of the weakly asymmetric simple exclusion process (WASEP) 
converges under appropriate rescaling to a continuous random field h which they charac¬ 
terize using the Cole-Hopf transform. Namely they show that the random field w = e h 
is the solution of a particular Ito SPDE : the stochastic heat equation (SHE) 

d w(t,x) = Aw(t, x)dt + w(t, x)dB t (x), (3) 

where Bt(x) = £(s, x)ds is a cylindrical Brownian motion in L 2 (R). The regularity 

of the solution of this equation does not allow to determine the intrinsic dynamics of 
h, but the convergence result shows that any candidate notion of solution to the KPZ 
equation should have the property that e h solves the SHE (3). 

The key tool in Bertini and Giacomin’s proof is Gartner’s microscopic Cole-Hopf 
transform [Gar 88 ] which allows them to prove the convergence of a transformed version 
of the WASEP to the stochastic heat equation. But the microscopic Cole-Hopf transform 
only works for this specific model (or at least only for a few special models, see the recent 
works [DT13,CT15] for extensions of the Bertini-Giacomin result) and cannot be used 
to prove universality. 

Another notion of solution for the KPZ equation has been introduced by Gongalves 
and Jara [GJ14] who proved tightness for a quite general class of rescaled weakly asym¬ 
metric exclusion processes and provided a description of the dynamics of the limiting 
random fields by showing that they must be energy solutions of the KPZ equation. In 
particular, they showed that (after a subtraction) the quadratic term in ( 1 ) would make 
sense only as a space-time distribution and not better. This notion of energy solution has 
been subsequently interpreted by Gubinelli and Jara in the language of controlled paths 
in [GJ13]. The key observation is that the fast time decorrelation of the noise provides 
a mechanism to make sense of the quadratic term. Unfortunately, energy solutions are 
very weak and allow only for few a priori estimates and currently there is no hint that 
they are unique. In this respect Gongalves and Jara’s result falls short of proving weak 
universality for systems without microscopic Cole-Hopf transformation. 

A related problem is studied in the the recent work of Alberts, Khanin and Quas- 
tel [AKQ14] where a universality result is rigorously shown for discrete random directed 
polymers with weak noise, which converge under rescaling to the continuum random 
directed polymer whose partition function is given by the stochastic heat equation. 

This was the state of the art until 2011, when Hairer [Hail3] established a well- 
posedness theory for the periodic KPZ equation using three key ingredients: a partial 
series expansion of the solution, an explicit control of various stochastic terms, and a 
fixed point argument involving rough path theory. As a final outcome he was able to show 
very explicitly that solutions h to equation ( 1 ) are limits for e —>• 0 of the approximate 
solutions h e of the equation 

d t h £ (t,x) = Ah s (t,x) + (Bh £ (t,x)) 2 - C £ + £ e (t,x), .x€l, t^z 0 , (4) 


4 



where T = R/(27rZ) and A is the Laplacian with periodic boundary conditions on T, £ e 
is a Gaussian process obtained by regularizing £ by convoluting it with a smooth kernel 
p, that is £ e = p £ * £ with p £ (x) = e _1 p(x/e), and where the constant C e has to be 
chosen in such a way that it diverges to infinity as e —>• 0. While it had been known 
since Bertini and Giacomin’s work [BG97] that the solutions to (4) converge to a unique 
limit, the key point of Hairer’s construction is that using rough path theory he is able 
to give an intrinsic meaning to the product (D h(t, x)) 2 and to obtain good bounds for 
that operation on suitable function spaces, which ultimately allowed him to solve the 
equation. 

The solution of the KPZ equation was one of the stepping stones in the development 
of Hairer’s regularity structures [Hail4], and now it is of course also possible to solve the 
equation using regularity structures rather than rough paths; see for example [FH14]. 

Using regularity structures, Hairer and Quastel [HQ15] recently proved the univer¬ 
sality of the KPZ equation in the following sense: they consider a regularized version £ e 
of the white noise and an arbitrary even polynomial P and show that under the right 
rescaling the solution h e to 


dthe = A h e + y/eP(Dhe) + ie (5) 

converges to the KPZ equation with (possibly) a non-universal constant in front of the 
non-linearity. Let us also mention the recent work [HS15], where a similar problem is 
studied except that P(x) = x 2 but £ e is not necessarily Gaussian. 

In the same period Hairer was developing his theory, together with P. Imkeller we 
proposed the notion of paracontrolled distributions [GIP15, GIP14, GP15] as a tool to 
give a meaning to and solve a large class of singular SPDEs. The class of models which 
can be analyzed using this method includes the generalized Parabolic Anderson model 
in 2 dimensions and the stochastic quantization equation in 3 dimensions and other 
interesting singular SPDEs [CC13, CF14, ZZ14b, ZZ14a, BB15,CC15,BBF15,PT15]. 

After Hairer’s breakthrough, another set of tools which should allow to solve the KPZ 
equation was developed by Kupiainen [Kupl4], who builds on renormalization group 
techniques to make sense of the three dimensional stochastic quantization equation </>|. 

The main aim of the present paper is to describe the paracontrolled approach to 
the KPZ equation. While the analysis of the KPZ equation itself is quite simple and 
shorter versions of the present paper circulated informally among people interested in 
paracontrolled equations, we kept delaying a public version in order to use the full 
versatility of paracontrolled analysis to explore various questions around KPZ and related 
equations and more importantly to find simple and direct arguments for most of the 
proofs. 

Indeed, despite paracontrolled calculus being currently less powerful than the full 
fledged regularity structure theory, we believe that it is lightweight enough to be effective 
in exploring various questions related to the qualitative behaviour of solutions and in 
particular in bridging the gap between stochastic and deterministic PDE theories. 

We start in Section 2 by recalling some basic results from paracontrolled distributions, 
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which we then apply in Section 3 to solve the conservative Burgers equation 

Jgf u = D u 2 + D£ 

driven by an irregular signal £ in a pathwise manner. In Section 4 we indicate how 
to adapt the arguments to also solve the KPZ equation and the rough heat equation 
pathwise, and we show that the formal links between these equations (u = D h, w = e h ) 
can indeed be justified rigorously, also for a driving signal that only has the regularity 
of the white noise and not just for the equations driven by mollified signals. 

We then show in Section 5 that by working in a space of paracontrolled distribu¬ 
tions one can find a natural interpretation for the nonlinearity in the stochastic Burgers 
equation, an observation which was not yet made in the works based on rough paths 
or regularity structures, where the focus is more on the continuous dependence of the 
solution map on the data rather than on the declaration of the nonlinear operations 
appearing in the equation. Our result allows us to build a partial bridge between energy 
solutions and paracontrolled solutions, although the full picture remains unclear. 

Section 6 extends our previous results to start the equations in more irregular initial 
conditions. In the case of the KPZ equation we can start in any function for /3 > 0, 
but the linear heat equation can be started in wo € Bp : 2o whenever 7 < 1/2 and p € 
[1, 00 ]; in particular we can take wq to be a fractional derivative of order 7 of the Dirac 
delta. 

In Section 7 we develop yet another approach to the KPZ equation. We show that 
its solution is given by the value function of the stochastic control problem 

h{t,x) = supE z HiS) + J (T(s> 7 s)-°o-^H 2 )ds , ( 6 ) 

where £ (s) = £(t — s ) and under E,„ we have 

7 s = x + f v r dr + V2W S 
Jo 

with a standard Brownian motion W. Such a representation has already proved very 
powerful in the case of a spatially smooth noise £, see for example [EKMS00, BCK14], 
but of course it is not obvious how to even make sense of it when £ is a space-time white 
noise. Based on paracontrolled distributions and the techniques of [DD15,CC15] we can 
give a rigorous interpretation for the optimization problem and show that the identity ( 6 ) 
is indeed true. Immediate consequences are a comparison principle for the KPZ equation 
and a pathwise global existence result. Recall that in [Hail3] global existence could only 
be shown by identifying the rough path solution with the Cole-Hopf solution, which 
means that the null set where global existence fails may depend on the initial condition. 
Here we show that this is not the case, and that for any uj for which £(<u) can be enhanced 
to a KPZ-enhancement (see Definition 4.1) and for any initial condition ho £ with 
/3 > 0 there exists a unique global in time paracontrolled solution h(ui) to the equation 

h(t , x, id) = (D h(t, x, d)) 2 — 00 + £(f, x, id), h(0, u) = ho. 


6 



and that the L°° norm of h is controlled only in terms of the KPZ-enhancement of the 
noise and the L°° norm of the initial condition. 

A surprising byproduct of these estimates is a positivity result for the solution of 
the SHE which is independent of the precise details of the noise. This is at odds with 
currently available proofs of positivity [Mue91, MF14, CK14], starting from the original 
proof of Mueller [Mue91], which all use heavily the Gaussian nature of the noise and 
proceed via rather indirect arguments. Ours is a direct PDE argument showing that 
the ultimate reason for positivity does not lie in the law of the noise but somehow in its 
space-time regularity (more precisely in the regularity of its enhancement). 

Section 8 is devoted to the study of Sasamoto-Spohn [SS09] type discretizations of the 
conservative stochastic Burgers equation. We consider a lattice model un- [0, oo) x Ttv —>• 
R (where Tjv = (27rZ/lV)/(27rZ)), defined by an SDE 

d UN(t,x) = x)dt + (DwBjv(uiv(t),ujv(t)))(x)dt + d(Djv£^ 1 // 2 Wjv(t, x)) 

Ujv(0,x) = Uq (x), 

where A jy, are approximations of Laplacian and spatial derivative, respectively, Bn 
is a bilinear form approximating the pointwise product, (W)v(i, ®))teK+,xeTjv is an N- 
dimensional standard Brownian motion, and Uq is independent of Wj\f. We show that if 
(uq ) converges in distribution in c to~ l3 for some /3 < 1 to uq, then (un) converges weakly 
to the paracontrolled solution u of 

Jz fu = D u 2 + cDu + D£, u( 0) = Uq, (7) 

where c € R is a constant depending on the specific choice of A tv, Dtv, and B n- If A tv 
and Dtv are the usual discrete Laplacian and discrete gradient and Bn is the pointwise 
product, then c = 1/2. However, if we replace Bn by the bilinear form introduced 
in [KS91], then c = 0. It was observed before in the works of Hairer, Maas and We¬ 
ber [HM12,HMW14] that such a “spatial Stratonovich corrector” can arise in the limit 
of discretizations of Hairer’s [Haiti] generalized Burgers equation A fu = g(u)d x u + £, 
and indeed one of the motivations for the work [HMW14] was that it would provide a 
first step towards the understanding of discretizations of the KPZ equation. To carry 
out this program we simplify many of the arguments in [HMW14], replacing rough paths 
by paracontrolled distributions. A key difference to [HMW14] is that here the “para¬ 
controlled derivative” is not constant, which introduces tremendous technical difficulties 
that were absent in all discretizations of singular PDEs which were studied so far, for 
example [ZZ14a,MW14]. We overcome these difficulties by introducing a certain random 
operator which we bound using stochastic computations. 

As an application of our convergence result, we show that the distribution of m + 
2 " 1 / 2 7 /, where r] is a space white noise and m € R, is invariant under the evolution of 
the conservative stochastic Burgers equation. While this is well known (see [BG97] and 
also the recent work [FQ14] for an elegant proof in the more complicated setting of the 
non-periodic KPZ equation), ours seems to be the first proof which does not rely on the 
Cole-Hopf transform. 
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In Section 9 we construct the enhanced white noise which we needed in our pathwise 
analysis. We try to build a link with the Feynman diagrams from quantum field theory. 
The required bounds are shown by reducing everything to a few basic integrals that 
can be controlled by a simple recursive algorithm. In Section 10 we indicate how to 
adapt these calculations to also obtain the convergence result for the enhanced data 
in the lattice models. We also calculate the explicit form of the correction constant c 
appearing in (7). 

Acknowledgments. We are very grateful to Khalil Chouk, who pointed out that by 
working in more general Besov spaces than Soo,oo one can start the rough heat equation 
in the Dirac delta. We would also like to thank Rongchan and Xiangchan Zhu, who 
found a mistake in a previous version of the paper. 

The main part of the work was carried out while N.P. was employed by Univer- 
site Paris Dauphine and supported by the Fondation Sciences Mathematiques de Paris 
(FSMP) and by a public grant overseen by the French National Research Agency (ANR) 
as part of the “Investissements d’Avenir” program (reference: ANR.-10-LABX-0098). 

2 Paracontrolled calculus 

Paracontrolled calculus and the relevant estimates which will be needed in our study of 
the KPZ (and related) equations have been introduced in [GIP15]. In this section we 
will recall the notations and the basic results of paracontrolled calculus without proofs. 
For more details on Besov spaces, Littlewood-Paley theory, and Bony’s paraproduct the 
reader can refer to the nice recent monograph [BCD11], 

2.1 Notation and conventions 

Throughout the paper, we use the notation a < b if there exists a constant c > 0, 
independent of the variables under consideration, such that a ^ c-b, and we write a ~ b 
if a < 6 and b < a. If we want to emphasize the dependence of c on the variable x, then 
we write a{x) < x b(x). For index variables i and j of Littlewood-Paley decompositions 
(see below) we write i < j if there exists N e N, independent of j, such that i ^ j + N 
(or in other words if 2 * < 2 J ), and we write i ~ j if i < j and j < i. 

An annulus is a set of the form = {x G M. d : a ^ |x|^ 6 } for some 0 < a < b. A 
ball is a set of the form SB = {x € M. d : |x|^ b}. T = M/ (27 rZ) denotes the torus. 

If / is a map from Aclto the linear space Y, then we write f s j = f(t) — f(s). For 
/ G LP{ T) we write ||/(®)||^ g(T) := f T \f(x)\ p dx. 

Given a Banach space X with norm ||-||x an d T > 0, we write CtX = C([0,T],X) 
for the space of continuous maps from [0, T] to X , equipped with the supremum norm 
||-||c T x, an d we set CX = C(M+, A), equipped with the topology of uniform convergence 
on compacts. Similarly C*(M, X) will always be equipped with the locally uniform topol¬ 
ogy. For a € (0,1) we also define C^X as the space of a-Holder continuous functions 



from [0, T] to X , endowed with the seminorm ||/||c“X= su Po<s<tcrll/s,tll^/l^ — s|“, and 
we write CXX for the space of locally a-Holder continuous functions from to X. 

The space of distributions on the torus is denoted by @'(Y) or &. The Fourier 
transform is defined with the normalization 

&u[k) = u(k) = f e~ l< ' k ’ x ' > u(x)dx, k G Z, 

J T 

so that the inverse Fourier transform is given by ^~ l v(x) = 'Yf,.e il ' k,x ' ) v(k). We 

denote Fourier multipliers by <p(D)u = ^ r ~ 1 (p^ r u) whenever the right hand side is well 
defined. 

Throughout the paper, (y, p) will denote a dyadic partition of unity such that 
supp(/9(2 - *-)) FI supp(p(2“J'-)) = 0 for |i — j\> 1. The family of operators 
will denote the Lit.tlewood-Paley projections associated to this partition of unity, that 
is A_iu = ^" -1 (x^u) and A j = (p(2~^-)^u) for j ^ 0. We also use the notation 

Sj = J2i<j A i- We write Pj = P(^~ j ') f° r 3 ^ 0 and p -1 = y, and Xj = x(2 _J ’-) for 
j ^ 0 and Xj = 0 for j < 0. We define if^{k,i) = JJj Xj-i(k)pj(£) and ip 0 {k,£) = 
Pi(k)Pj(£)■ The Holder-Besov space -B£j )00 (T, M) for a € R will be denoted by 
c £ cx and equipped with the norm 

ll/lla = ll/ll B “ j00 = SU P (2* a '||A J ;/|| L ^ (T) ). 
i^—1 

If / is in ( ^ >a ~ £ for all e > 0, then we write / £ r io a ~ . For a € (0,2), we also define the 
space = C^ 2 L°° fl Ct ^ a , equipped with the norm 

ll/IU“= max { \\f\\ c a/2 Loo , ||/||c T y«} • 

The notation is chosen to be reminiscent of Af = dt — A, by which we will always 
denote the heat operator with periodic boundary conditions on T. We also write 2 zf“ = 

C \Jc L °° n C ^°- 

2.2 Bony—Meyer paraproducts 

Paraproducts are bilinear operations introduced by J. M. Bony [Bon81,Mey81] in order 
to linearize a class of nonlinear PDE problems. Paraproducts allow to describe functions 
which “look like” some given reference functions and to perform detailed computations 
on their singular behavior. They also appear naturally in the analysis of the product 
of two Besov distributions. In terms of Lit.tlewood -Paley blocks, the product fg of two 
distributions / and g can be (at least formally) decomposed as 

fg= Y Y A if A J9 = f ^9 + fyg + fog. 

j^-li^-1 

Here f <g is the part of the double sum with i < j — 1, and f >~ g is the part with 
i > j + 1, and fog is the “diagonal” part, where |i — 1. More precisely, we define 

j -2 _ 

f<g=9^f=Y Y A if A j9 and fog= ^ AifAjg. 
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This decomposition behaves nicely with respect to Littlewood-Paley theory. Of course, 
the decomposition depends on the dyadic partition of unity used to define the blocks 
and also on the particular choice of the set of pairs (i. j) which contributes 
to the diagonal part. Our choice of taking all (i,j) with \i — j\^ 1 into the diagonal part 
corresponds to a dyadic partition of unity which satisfies supp(/ 9 ( 2 _ *-))nsupp(/ 9 ( 2 _ - ? -)) = 
0 for \i—j\> 1. This implies that every term Sj-ifAjg in the series / -< g = Yhj Sj-ifAjg 
has a Fourier transform which is supported in an annulus and of course the same 
holds true for / g. On the other side, the terms in the diagonal part fog have Fourier 
transforms which are supported in balls. We call f ~<g and / >- g paraproducts , and fog 
the resonant term. 

Bony’s crucial observation is that f <g (and thus f F g) is always a well-defined 
distribution. In particular, if a > 0 and ft £ M, then (/, g) / -< g is a bounded bilinear 
operator from if" x if ^ 3 to Heuristically, f~<g behaves at large frequencies like g 
(and thus retains the same regularity), and / provides only a modulation of g at larger 
scales. The only difficulty in defining fg for arbitrary distributions lies in handling the 
resonant term fog. The basic result about these bilinear operations is given by the 
following estimates. 

Lemma 2.1 (Paraproduct estimates). For any ft £ M we have 

\\f<9b<p \\f\M\9b, ( 8 ) 

and for a < 0 furthermore 

\\f <g\U+p< a ,p WfWabb- ( 9 ) 

For a + ft > 0 we have 

\\f ° g\\a+/3<a,P WfWahb- ( 10 ) 

A natural corollary is that the product fg of two elements / £ if" and g € is well 
defined as soon as a + ft > 0 , and that it belongs to ^ 7 , where 7 = min{a, ft, a + ft}. 

2.3 Product of paracontrolled distributions and other useful results 

The key result of paracontrolled calculus is the fact that we are able to multiply para¬ 
controlled distributions, extending Bony’s results beyond the case a + ft > 0. We present 
here a simplified version which is adapted to our needs. Let us start with the following 
meta-definition: 

Definition 2.2. Let ft > 0 and a £ R. A distribution / £ c £ a is called paracontrolled 
by u £ if" if there exists f £ if ' 5 such that /** = / — f < u £ if" +/3 . 

Of course in general the derivative f is not uniquely determined by / and u, so more 
correctly we should say that (/, /') is paracontrolled by u. 
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Theorem 2.3. Let a,/? € (1/3,1/2). Let u £ c £ a , v £ and let (/,/') 5e para- 

controlled by u and ( g,g') be paracontrolled by v. Assume that uov £ c € 2a ~ x is given as 
limit of {u n ov n ) in where (u n ) and (y n ) are sequences of smooth functions that 

converge to u in c lo a and to v in c £ a ~ x respectively. Then fg is well defined and satisfies 

Wf9-f<9ha-l< (||/1l/3||w|U + ||/ t ||a+/3)(lb / ||/3||^IU-l + lb ti |U+/3-l) + ll/Vll/3||« 0 ^l|2a-l. 

Furthermore, the product is locally Lipschitz continuous: Let u £ c £ a , v £ c £ a ~ 1 with 
uov £ c £ 2a ~ x and let (/,/') be paracontrolled by u and ( g,g') be paracontrolled by 
v. Assume that M > 0 is an upper bound for the norms of all distributions under 
consideration. Then 

\\{fg-f<g)-{fg-f<g)h a -i< ( 1 +M 3 )[\\f' - f'\\p+\\g' - g'\\g+\\u-u\\ Q +\\v-v\\ a -i 

+ 11 / # - P\\a+P + \\g^ - gHa+P-l + \\u°V - UOv\\ 2a -l\. 

// /' = /' = 0 or g' = g 1 = 0, then M 3 can be replaced by M 2 . 

Remark 2.4. The proof is based on the simple result of [GIP15] that the commutator 

C(f , g, h) := (/ + g) o h - f(g o h) 

is a bounded trilinear operator from x < ^' a x ^" _1 to ^ 2a+/3_1 . 

We will write / • g instead of fg whenever we want to stress the fact that we are 
considering the product of paracontrolled distributions. For some computations we will 
also need a paralinearization result which allows us to control nonlinear functions of the 
unknown in terms of a paraproduct. 

Lemma 2.5 (Bony—Meyer paralinearization theorem). Let a € (0,1) and let F £ C 2 . 
There exists a locally bounded map Rp : c & a —> c € 2a such that 

F{f) = F\f)<f + R F {f) (11) 

for all f £ < tf a . If F £ C 3 , then R F is locally Lipschitz continuous. 

We will also need the following lemma which allows to replace the paraproduct by a 
pointwise product in certain situations: 

Lemma 2.6 (A further commutator estimate). Let a > 0, £ R, and let f,g £ < & a , 

and h € Then 


II f<(g<h) - (fg) -< h\\ a+ p< H/IUII0IUHI/9- 

When dealing with paraproducts in the context of parabolic equations it would be 
natural to introduce parabolic Besov spaces and related paraproducts. But to keep a 
simpler setting, we choose to work with space-time distributions belonging to the scale 
of spaces (C , ^ >a ) QgR . To do so efficiently, we will use a modified paraproduct which 
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introduces some smoothing in the time variable that is tuned to the parabolic scaling. 
Let therefore <p € C°°(M, R+) be nonnegative with compact support contained in M + 
and with total mass 1, and define for alii ^ — 1 the operator 

Qi : CV? -> CtfP, Qif(t ) = [ 2~ 2i ip(2 2i (t - s))f(s V 0)ds. 

J R 

We will often apply Qi and other operators on CtfP to functions / € Ct which we 
then simply extend from [0, T] to M + by considering /(• A T ). With the help of Qi, we 
define a modified paraproduct 


f^9 = ^2(QiSi-if)Ai9 

l 

for f,g € C (R + , It is easy to see that for f we have essentially the same 

estimates as for the pointwise paraproduct f ~<g, only that we have to bound / uniformly 
in time. More precisely: 

Lemma 2.7. For any f3 € M we have 

\\f^9(t)b<\\f\\c t L^\my, (12) 

for all t > 0, and for a < 0 furthermore 

\\f ^<g(t)\\a+i3< Wf\\ctV a \\9(t)b- ( 13 ) 

We will also need the following two commutation results. 

Lemma 2.8. Let a € (0, 2) and /3 € M. Then 

\\(f^<9-f-<g)(t)\\a+/3< ||/||if“lb(i)||/3 
for all t ^ 0. If a £ (0,1), then we also have 

II (■Z’if^g) - /-« (Sfg)) (0IU+/3-2 ll/ll^ t “lb(i)ll/3- 

For a proof see [GIP15]. As a consequence, every distribution on M + x T which 
is paracontrolled in terms of the modified paraproduct is also paracontrolled using the 
original paraproduct, at least if the derivative is in Jz?" and not just in 

Moreover, we introduce the linear operator / : C (M+, ^(T)) —>• C (R+, 3>'{T)) given 

by 

If(t) = [ P t - S f(s)ds, 

Jo 

for which we have standard estimates that are summarized in the following Lemma. 


12 



Lemma 2.9 (Schauder estimates). For a G (0,2) we have 


l|//IU~< WfWcrVc -2 (14) 

for all T > 0, as well as 

\\s ^ P s u 0 \\^< \\u 0 \\ a . (15) 

Combining the commutator estimate above with the Schauder estimates, we are able 
to control the modified paraproduct in 2*?^ spaces rather than just Ct c ^ ol spaces: 

Lemma 2.10. Let a G (0,2), 5 > 0 and let f G g G Ct and 2z ?g G Ct^ 01-2 ■ 
Then 

11 / \\f\\^f{\\9\\c T ^ a + W-^QWct^- 2 ) ■ 

Proof. We apply the Schauder estimates for ££ to obtain 


\\f^g\W« < ||/^<£f(0)|k-+||^(/^<£f)ll C T ^“- 2 

< \mf^g)-f^J?g\\ Ct^+S-2 + 11/^ {seg )II Ct^ q - 2 

< ll/lljsrf(Nlc T «’“+ W-^sWct^- 2 )- 

□ 

It is possible to replace ||/||_^i on the right hand side by ||/||c t l°°, see Lemma 2.16 
of [Fur 14]. But we will only apply Lemma 2.10 in the form stated above, so that we do 
not need this improvement. 

At the end of this section, let us observe that in (jSf^) spaces we can gain a small 
scaling factor by passing to a larger space. 

Lemma 2.11. Let a G (0, 2), T > 0, and let f G S£fjf - Then for all 5 G (0, a] we have 

11/lljf‘S ll/(o)IU+r < “- i > /2 ||/||j~. 

Proof. It suffices to observe that by interpolation \\f\\ c (a-s)/2^ s < ||/||jsf“- Cl 


3 Stochastic Burgers equation 

3.1 The strategy 

Instead of dealing directly with the KPZ equation we find it convenient to consider its 
derivative u = D h which solves the Stochastic Burgers equation (SBE) 

Ffu = Du 2 + D£. (16) 

In a later section we will show that the two equations are equivalent. The strategy we 
adopt is to consider a regularized version for this equation, where the noise £ has been 
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replaced by a noise that is smooth in space. Our aim is then to show that the solution 
u £ of the equation 

33u e = Du^ + D£ e 

converges in the space of distributions to a limit u as e —>• 0. Moreover, we want to 
describe a space of distributions where the nonlinear term Drt 2 is well defined, and we 
want to show that u is the unique element in this space which satisfies 33 u = Du 2 + D£ 
and has the right initial condition. 

That it is non-trivial to study the limiting behavior of the sequence (u e ) can be 
understood from the following considerations. First, the limiting solution u cannot be 
expected to behave better (in terms of spatial regularity) than the solution X of the 
linear equation 

= D£ 

(for example with zero initial condition at time 0). It is well known that, almost surely, 
X € CV' 1 / 2- but not better, at least on the scale of (C ( 3°) a spaces. In that case, the 
term u 2 featured on the right hand side of (16) is not well defined since the product 
(/, 9 ) fg is a continuous map from CV a x to CV“ A/3 only if a + (5 > 0, a 
condition which is violated here. Thus, we cannot hope to directly control the limit of 
(u 2 ) e as s —>• 0. 

What raises some hope to have a well-defined limit is the observation that if X e is 
the solution of the linear equation with regularized noise £ e , then {DX 2 ) £ converges to 
a well defined space-time distribution DJ 2 which belongs to JzfC^ 0- (the space of all 
distributions of the form S3v for some v € C"^ 0- ). This hints to the fact that there 
are stochastic cancellations going into DX 2 due to the correlation structure of the white 
noise, cancellations which result in a well behaved limit. Taking these cancellations 
properly into account in the full solution u is the main aim of the strategy we will 
implement below. 

In order to prove the convergence and related results it will be convenient to take a 
more general point of view on the problem and to study the solution map ^> r b e : (9, v) i->- u 
mapping well-behaved (for example smooth, but here it is enough to consider elements 
of C'(M+,C <1 ) and C 1 respectively) functions 6 and initial conditions v to the classical 
solution of the evolution problem 

ifu = Du 2 + D0, u( 0)=v, (17) 

which we will call the Rough Burgers Equation (RBE). 

To recover the original problem our aim is to understand the properties of <h when 6 
belongs to spaces of distributions of the same regularity as the space-time white noise, 
that is for 9 € C9o~ l / 2 ~. In order to do so we will introduce an n-tuple of distributions 
X(0) constructed recursively from the data 9 of the problem and living in a complete 
metric space X r i )e and show that for every bounded ball S3 around 0 in df r b e and every 
bounded ball 33 in there exist T > 0 and a continuous map T r b e : 33 x 33 —>• 

satisfying <h r b e (#,u) = T r b e (X(0),u)|[ 0j ;r] f° r aH (X(#),u) € 33 x 33 . This will 
allow us to obtain the local in time convergence of u £ = 3> r be(£e) <^(eD)uo) by showing 
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that the sequence (X(£ e )) e converges as £ ^ 0 to a well-defined element of X T ^ e : here ip 
is a compactly supported, infinitely smooth function with <p(0) = 1. That last step will 
be accomplished via direct probabilistic estimates with respect to the law of the white 
noise £. The described approach thus allows us to decouple the initial problem into two 
parts: 

1 . a functional analytic part related to the study of particular spaces of distributions 
that are suitable to analyze the structure of solutions to equation (17); 

2. a probabilistic part related to the construction of the RBE-enhancement X(£) € X T ^, e 
associated to the white noise £. 

3.2 Structure of the solution 

In this discussion we consider the case of zero initial condition. Let us first define a 
linear map J(f) for each smooth / on M + x T as the solution to 2z?J(/) = D/ with 
initial condition J(/)(0) = 0. Explicitly, we have 

= /(D/)(t) = f P t . s B(f(s))ds, 

Jo 

where (Ppt^o is the semi-group generated by the periodic Laplacian on T. Using the 
estimates on I given in Lemma 2.9, we see that J is a bounded linear operator from 
C T tf a to C T tf a+1 for all a G M and T > 0. 

Let us now expand the solution to the rough burgers equation (17) around the solu¬ 
tion X = J (9) to the linear equation PPX = D 9. Setting u = X + u ^ 1 , we have 

PPvP 1 = D(u 2 ) = D(A 2 ) + 2D (I^ 1 ) + D ((i^ 1 ) 2 ). 

We can proceed by performing a further change of variables in order to remove the term 
D(A 2 ) from the equation by setting 

u = X + J(X 2 ) + u? 2 . (18) 

Now uP 2 satisfies 

^uP 2 = 2D(XJ(X 2 )) + D{J(X 2 )J(X 2 )) 

+2B(XiP 2 ) + 2B(J(X 2 )vP 2 ) + B{(vP 2 ) 2 ). 1 

We can imagine to make a similar change of variables to get rid of the term 2D (XJ(X 2 )). 
As we proceed in this inductive expansion, we generate a certain number of explicit terms, 
obtained via various combinations of the function X and of the bilinear map 

B(f,g) = J{fg) = [ P- s B(f{s)g(s))ds. 

Jo 

Since we will have to deal explicitly with at least some of these terms, it is convenient 
to represent them with a compact notation involving binary trees. A binary tree r € T 
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is either the root • or the combination of two smaller binary trees r = (71T2), where the 
two edges of the root of r are attached to n and T2 respectively. The natural grading 
d:T —>• N is given by d(») = 0 and d{{r\T 2 )) = 1 + d[r 1) + d(j 2 )- However, for our 
purposes it is more convenient to work with the degree d(r), which we define as the 
number of leaves of r. By induction we easily see that d = d + 1. Then we define 
recursively a map A : T —>• C (M+, &) by 

X* = X, X( T1T2 ) = B{X T1 ,X T2 ), 

giving 

A V = H(A,A), X* = £(A,A V ), X^ = B(A,A*), A V = B(X V , A v ) 
and so on, where 

(••) = V, (V.)=V (.ty=\ (VV) = XX, 

In this notation the expansion (18)-(19) reads 

u = X + X v + v? 2 , (20) 

u^ 2 = 2A X ' + A V + 2 B(X, u? 2 ) + 2B{X V , vfi 2 ) + B(u^ 2 ,u^ 2 ). (21) 

Remark 3.1. We observe that formally the solution u of equation (17) can be expanded 
as an infinite sum of terms labelled by binary trees: 

u = ^c(r)A r , 
reT 

where c(r) is a combinatorial factor counting the number of planar trees which are iso¬ 
morphic (as graphs) to t. For example c(«) = 1, c(V) = 1, c(X r ) = 2, c{%) = 4, 
c(V) = 1 and in general c(r) = ^ to£ j- l( rir2 ) =r c(Ti)c(T2). Alternatively, truncating 
the summation at trees of degree at most n and setting 

u = J2 c{t)X t + u> n , 

rST ,d(r)<n 

we obtain a remainder u^ n that satisfies the equation 

u^ n = c(ri)c(T 2 )X (Tir2 )+ ^ c(T)B(X T ,u^ n ) + B(u >n ,u^ n ). 

ri,T2 : d(n) < n,d(T2) < n T:d(r)<n 

d([riT 2 )) ^ n 


Our aim is to control the truncated expansion under the natural regularity assump¬ 
tions for the white noise case. These regularities turn out to be 

lecr 1/2 ', A v €C^°", A*,A^ € C^ 1/2 ~, X v eCV 1 -. (22) 
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Moreover, for any / € C r <o a with a > 1/2 we have that B(X, f) G CV 1 / 2 : indeed, the 
bilinear map B satisfies 


\\B(f,g)\\c T vs< II/IIct^IIs'IIct^/ 3 ( 23 ) 

for any T > 0, where 6 = min(a,/3,a + f3) + 1, but only if a + f3 > 0. Equation (21) 
implies that u^ 2 has at least regularity (1/2)—. In this case B{v^ 2 ,u^ 2 ) is well defined 
and in C"^ 3 / 2- , but B(X,u^ 2 ) is not well defined because the sum of the regularities 
of X and of u^ 2 just fails to be positive. To make sense of B(X,u^ 2 ), we continue the 
expansion. Setting u^ 2 = 2X ^ + u^ 3 , we obtain from (21) that u^ 3 solves 

v ? 3 = 4X^ + X V + 2 B{X, u? 3 ) + 2H(X V , v? 2 ) + B{u> 2 , v? 2 ). 

At this stage, if we assume that u^ 3 ,u^ 2 € C'if 1 / 2 ', all the terms but 4X^ + 2 B(X,u^ 3 ) 
are of regularity indeed 

B(X v ,v? 2 ) € CV 1 ", B{v? 2 ,u> 2 ) g CV 3/2 “, 

and we already assumed that X ^ . This observation implies that u^ 3 has the 

form 

i? 3 = 2 B{X,v? 2 ) + 4X^ + Ctf 1 -, 

by which we mean u^ 3 — 2 B(X,u^ 2 ) — 4X^ € C^ 1- . Of course, the problem still 
lies in obtaining an expression for B(X,u^ 2 ) which is not well defined since our a pri¬ 
ori estimates are insufficient to handle it. This problem leads us to finding a suitable 
paracontrolled structure to overcome the lack of regularity of u^ 2 . 

3.3 Paracontrolled ansatz 

Here we discuss formally our approach to solving the rough Burgers equation (IT). The 
rigorous calculations will be performed in the following section. 

Inspired by the partial tree series expansion of u, we set up a paracontrolled ansatz 
of the form 

u = X+ X V + 2X*+ u Q , u Q = v! + (24) 

where the functions u',Q and vfi are for the moment arbitrary and later will be chosen 
so that u solves the RBE (17). The idea is that the term v! -4< Q takes care of the less 
regular contributions to u of all the terms of the form B(X, /). In this sense, we assume 
that Q G CV 1//2 ~ and that u' € CV 1 / 2- nC' 1 / 4 ~L°°. The remaining term iJ will contain 
the more regular part of the solution, more precisely we assume that u$ € C"^ 1- . Any 
such u is called paracontrolled. 

For paracontrolled distributions u, the nonlinear term takes the form 
Du 2 = D(XX+2XX v +X v X v +4XA^)-|-2D(Xu <2 )-|-2D(X v (u < ^-|-2X^))+D((u < ^-|-2X^) 2 ). 
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Of all the term involved, the only one which is problematic is D (Xifi), since by our 
assumptions 2D (X v (u® + X^)) + D((u^ + A^) 2 ) are well defined and 

D(XX + 2XX V + X V X V + AXX*) = Jf(X v + 2X* + X V + 4x\ 

But using the paracontrolled structure of u, the term T)(Xifl) is defined as D(A • u9) 
provided that QoX € CV 0_ is given. Under this assumption we easily deduce from 
Theorem 2.3 that 

D(Xtt Q ) -u^DIe C^~ x ~. 

In other words, for paracontrolled distributions the nonlinear operation on the right 
hand side of (17) is well defined as a bounded operator. 

Next, we should specify how to choose Q and which form v! will take for the solution 
u of (17). To do so, we derive an equation for vr: 

=^(u-X-X v - 2X X ' -u^Q) = 4D( AA^) + 2u Q -X DA -u! -«££Q + CV" 1 ', 

where we used the commutator estimate Lemma 2.8 to replace J£{u' -« Q) by v! -« ££Q. 
If we now assume that X^ o X e C*^ 0- , then we obtain 

2 = (2 u Q + AX^) -« BX - v! + Ctf- 1 -, 

which shows that we should choose Q so that JfQ = TAX and v! = 2ifi + 4A^. While 
this equation looks weird, it can be easily solved by a fixed point argument. However, in 
the next section we take a different approach and start with the solution to the equation 
driven by regular data (of which we know that it exists), and use it to derive a priori 
estimates which allow us to pass to the limit. 

3.4 A priori estimates on small time intervals 

We are now ready to turn these formal discussions into rigorous mathematics and to 
solve Burgers equation on small time intervals. For simplicity, we first treat the case 
of initial conditions that are smooth perturbations of the white noise. In Section 6 we 
indicate how to adapt the arguments to treat more general initial conditions. 

Throughout this section we fix a € (1/3,1/2). First, let us specify which distributions 
will play the role of enhanced data for us. We will usually consider regularized versions 
of the white noise which leave the temporal variable untouched and convolute the spatial 
variable against a smooth kernel. So in the following definition we are careful to only 
require a smooth X, but not a smooth 9. 

Definition 3.2. (RBE-enhancement) Let 

A rbe C CV"" 1 x CV 2 "" 1 x x J^ 2a x j^ 2a x C'T 2 "" 1 
be the closure of the image of the map 0 r b e : Jzf'C/M, C'°°(T)) —>• T rbe given by 

0 rbe (0) = (x(0),x v (e),x\d),x\e),xV(e),(Qox)(0)), ( 25 ) 
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where 


SfX(0) 

= D0, 

^X v (0) 

= D (X{9f), 

SfX*(0) 

= D (X(0)X^(0)), 

&x\9) 

= D(A*(0)o X(0)), 

A£X^{9) 

= D (X V (9)X V (9)), 

if Q{9) 

= D X{0), 


( 26 ) 


all with zero initial conditions except 


X(9) for which we choose the “stationary” initial 


condition 


X(0)(O) = [ P- S m{s)ds. 
J — oo 


We call © r be($) the RBE-enhancement of the driving distribution 9. For T > 0 we define 
Af r b e (T) = ^rbe|[o,T] and we write ||X||^ rbe ( T ) for the norm of X in the Banach space 
C T tf a - 1 X Cr^ 2 "" 1 x ifp x if|“ x A£ 2a x C T ^ 2a ~ l . Moreover, we set d* rbe(T) (X,X) = 


For every RBE-enhancement we have an associated space of paracontrolled distribu¬ 
tions. 


Definition 3.3. Let X € T r b e and 5 > 1 — 2a. We define the space £^ be = °f 

distributions paracontrolled by X as the set of all (u, u', v$) € CV " _1 x i ? 5 x i ?“ +<5 such 


that 

u = X + X v + 2X X ' W-KQ + ii 11 . 


For 6 = a we will usually write SW = ^b e - For ever y F > 0 we set ^rbe( T ) = ^rbel[0,T], 
and we define 


ll W H^rbe( T ) _ ll u, Hif| + ll wfl llif“ +4 - 

We will often use the notation u® = u! ~¥.Q + uK Since u is a function of u' and u we 
will also write (v!,v$) € fi? r be> and if there is no ambiguity then we also abuse notation 
by writing u € £^ r be- We call u' the derivative and vfi the remainder. 


Strictly speaking it is not necessary to assume ifi € i? a+<5 , the weaker condition 
vfi € C’t£ > ° 1+s would suffice to define Du 2 . But this stronger assumption will simplify 
the notation and the uniqueness argument below, and it will always be satisfied by the 
paracontrolled solution to Burgers equation. 

Our a priori estimate below will bound ||w||^ rbe (T) in terms of the weaker norm 
IMI^h (T) f° r suitable 5 < a. This will allow us to obtain a scaling factor T £ on small 
time intervals [0, T], from where we can then derive the locally Lipschitz continuous 
dependence of the solution u on the data (X, uq). Throughout this section we will work 
under the following assumption: 


Assumption (T,M). Assume that 9,9 G J§?C(R, C* 00 ^)) and uo,uq € to 2a , and that u 
is the unique global in time solution to the Burgers equation 


££u = Drt 2 + D£, u(0) = X(0) +uq. 


(27) 
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We define X = 0 r b e ($), v9 = u — A — A v — 2A ^, and v! = 2vP + AX**, and we set 
u t = U Q _ u > q Similarly we define X, u, zfi, u', uK Finally we assume that T, M > 0 
are such that 

max {||«o ha, ll^olha, ll X IU rbe (T) 5 IWU rbe (T)} ^ M - 
Lemma 3.4. Make Assumption (T,M) and let 5 > 1 — 2a. Then 

\\J?J - ifA V - A^x\ CT ^ a - 2 +\\u'\\^< (1 + M 2 )(l + \\u\\% s (T) ). (28) 

If further also ||«||^ rbejX (T), IHI^ rbe ,x( T )^ M , then 

\\(^J - 2^A V - 4ifA^) - (AFu^ - 2^A V - A^x\\ CT ^ a -2+\\u' - v! \\<g« (29) 

< M 2 (^ rbe(T) (X,X) + (\\u' -u'Wjrt + h* -uHjq?))- 

Proof. Let us decompose ifiJ into a sum of bounded multilinear operators: 

Jgfu # = Du 2 - DA 2 - 2DA V A - JS?(u' -« Q) 

= if A V + 2D[u Q A -u Q PX} + 2[D (u q <X)-u q P DA] 

+ 4if A^ + ADiX* y A] + 4[D(A Xr -< A) — A* X DA] + [u y DA - u' -« DA] 

+ D[2A V (2A X ' + u Q ) + (fix'** + u Q ) 2 ] - [if (u' -« Q) - u ' (ifQ)] , 

where we used that u' = 2rt^ + 4A^. From here it is straightforward to show that 

\\J?J - ifA V - 4if x\ Cr ^ a -2 < ||X|U rbe(T) (||^||c T ^ + ||«'|k T “) + \\u Q \\ 2 CtL - 

(30) 

+ dl x llA: rbe (r)+ll x ll|’ rbe (r))( 1 + \\ u ^\\c T v a + s +\\ u '\\cT , tf 5 )- 

Now \\u q \\ Ct ^< \\u'\\^s\\X\\ XrheiT) +\\J\\cT^ by Lemma 2.10 (observe that ||Q||^< 

11X11 x T i, e (T) by the Schauder estimates, Lemma 2.9), and therefore the estimate (28) 
follows as soon as we can bound (l+||X|| i y rbe (' r ))||tt'|| by the right hand side of (28). For 

this purpose it suffices to use the explicit form v! = 2v9 +4A^ and then u® = u! -4< Q+u ** 
in combination with Lemma 2.10. 

The same arguments combined with the multilinearity of our operators also yield (29). 

□ 


The main result of this section now follows: 

Theorem 3.5. For every (X, uo) € A r b e x ( tf 2a there exists T* € (0,oo] such that for all 
T <T* there is a unique solution (u,u',u^) € ^ r b e (T) to the rough Burgers equation (27) 
on the interval [0, T]. Moreover, we can choose 

T* = infjt ^ 0 : IMb rbe (t)= oo}. 
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Proof. Let M = 2 max{||«o||2a, ||^||;t rbe (i)}- Consider a sequence (O n ,v,Q) C _§fC(R, C°°)x 
c £ 2a such that X n = 0 r be(#n) approximates X|[ 0 ,t] i n ^rbe(X') for all T > 0 and (uq) ap¬ 
proximates uo in ^ 2 “. Without loss of generality assume that 11^0 ||2aV||X n || i% - rbe ( 1 )^ M 
for all n. Denote by u n the solution to Burgers equation driven by X n and started in Uq 
and define Un, u' n , uh as in Assumption ( T,M ). In particular, for all n we are in the 
setting of Assumption (1 , M). 

For hxed n, the remainder has the initial condition uh{ 0) = Uq G c € 2a , so that 
Lemma 3.4 and the Schauder estimates (Lemma 2.9) yield 

ll' u n||® rbe (r)^ (1 + M )(1 + ||'lt n ||^ / s) 

for all r G (0,1]. Therefore, we can apply Lemma 2.11 to gain a small scaling factor: 

ll“«ll^ rbe (r)~ (1 + M 2 )(l + T a S \Wn\\% rbe ( T ))- (31) 

Now ||M n ||^ rbe ( T )^ ||'Un||^ rbe (i) for ah r G (0,1], and therefore there exists a universal 
constant C > 0 such that ||^ra||^ rbe (r)^: C( 1 + M 2 ) for all sufficiently small r ^ r*. 
However, a priori t* may depend on stronger norms of 0 n and not just on ||X n ||^ rbe ( 1 ). 
To see that we may choose t* only depending on M, it suffices to note that r i->- 
||n n ||^ rbe ( T ) is a continuous function (which follows from the fact that uh and u' n are 
actually more regular than «Sf 2 “ and 2z?" respectively), to select the smallest t* with 
IKHsw(t*)= C(l + M 2 ), and to plug ||« n |k rbe (r*) into (31)- 

In other words ||'U n ||^ rbe ( T )^ C( 1 + M 2 ) for all n, so that the second statement of 
Lemma 3.4 in conjunction with the Schauder estimates for the Laplacian shows that 

Hi — u h\\^ a + \\ u 'm ~ I 1 + M2 )(Hbe{T)(^rn^n) + 11 ~ Uq 112o) 

for all m,n (possibly after further decreasing r = r(M)). 

So if 

® = {(©rbe((?),«o) : 0 G SfC(R,C°°(J)),v G ||0 rb eWIU rbe (i), M} , 

and if : SB —> 2z?5 M j x ^‘Hi) I s the ma P which assigns to ( 6,v ) the local solution 
(u r , uf ) of (27), then T has a unique continuous extension from to its closure in 
A r be(l) x ^ 2 ", which is just the ball be(l) x c & 2a ) of size M in this space. By 

construction, u = (u' ,w) = 4'(X, uq) solves Burgers equation driven by X and with 
initial condition uq if we interpret the product Du 2 using the paracontrolled structure 
of u. To see that u is the unique element in ^ r b e,x{r(M)) which solves the equation, it 
suffices to repeat the proof of estimate (29). 

It remains to iterate this argument in order to obtain solutions until the blow up time 
T* . Assume that we constructed a paracontrolled solution u on [0, r] for some r > 0 
and consider the time interval [r, r + 1], Denote the solution on this interval by ft, and 
write X(t) = X(t + r). The initial condition for u is ft(0) = u(r) G so that we are 
not in the same setting as before. But in fact we only needed a smooth initial condition 
to obtain an initial condition of regularity 2a for the remainder, and now we have 

u # (0) = u(0) - X(0) - X V (0) - 2X Xr (0) - u' -« Q( 0) 
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= u(t) — X(0) — X v (t) — 2 X^(t) — u'(t) -< Q{t). 

According to Lemma 2.8 we have u'(t) -< Q(t) — v! -« Q(t) E c € 2a , so that the paracon- 
trolled structure of u on the time interval [0, r] shows that the initial condition for vf is in 
c € 2a . As a consequence, there exists Af > 0 which is a function of ||X||^ rbe ( T+1 ), |M|^ rbe ( r ) 
such that the unique paracontrolled solution u can be constructed on [0, r(M)]. We then 
extend u and v! from [0, r] to [0, r + r(M)] by setting u{t) = u(t — r) for t ^ r and 
similarly for v!. Then by definition v! E ^ . For the remainder we know from our 

construction that on the interval [r, r + r(M)] we have u® — u' -« T Q E Jz? 2a , where 
u' -4< T Q(t ) = f 2 2 i ip(2 2 i (t — s))S'j_iu , (s V r)dsAjQ(t). 

j J — OO 

It therefore suffices to show that (v! -« T Q — u' -«Q) L r t+t (m)]E «Sf 2a . But we already 
showed that (u' -« r Q — v! -« Q)(r) E ^ 2 ", and for t E [r, r + r(M)] we get 

^(u' ^ r Q - u' ^Q){t) = {v! ^ r ^Q -u'^SeQ){t) + tf 2a 

= («' -« r DA - it -« D X){t) + tf 2a . 

Now Lemma 2.8 shows that ^ DL — uADJ E C T+T ^)^ 2a ~ 2 , and by the same 
arguments it follows that also (u'^! t D1 — v! < DA E < ^ 2 “~ 2 )|[ r t + t (m)]E C([r, r + 
r(M)], < ^ 2a ), so that the 2z? 2 " regularity of on [0, r + r(M)] follows from the Schauder 
estimates for the heat flow, Lemma 2.9. 

Reversing these differences between the different paraproducts, we also obtain the 
uniqueness of u on [0, r + f(M)]. □ 

Remark 3.6. There is a subtlety concerning the uniqueness statement. It is to be 
interpreted in the sense that there exists (locally in time) a unique fixed point (u, u' ,v0) 
in ^ r be for the map S> r \,e 3 (v,v',v^) i->- (w,w',w^) E & r be, where 

w = P.(u 0 ) + I(Dv 2 ) + X, w' = 2v Q + 4X y ', w i = w - X - X w -2X^ - w' ^Q. 

For general X E A r b e it is not true that there is a unique u E Or^" -1 for which there 
exist ufifi with (u, u',vf) E ^ r b e,x(F) and such that u = P.(uq) + /(Du 2 ) + X, where 
the square is calculated using u' and uK The problem is that v! and u$ are in general 
not uniquely determined by u (see however [HP13,FS13] for conditions that guarantee 
uniqueness of the derivative in the case of controlled rough paths). Nonetheless, we will 
see in Section 5 below that this stronger notion of uniqueness always holds as long as 
X is “generated from convolution smoothings”, and in particular when X is constructed 
from the white noise. 

Moreover, without additional conditions on X E A r b e there always exists a unique 
u E Ct so that if (9 n ) is a sequence in J£C( R, C °°) for which (0 r be(#n)) converges 
to X, then the classical solutions u n to Burgers equation driven by 9 n converge to u. 
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The solution constructed in Theorem 3.5 depends continuously on the data X and 
on the initial condition uq. To formulate this, we have to be a little careful because for 
now we cannot exclude the possibility that the blow up time T* is finite (although we 
will see in Section 7 below that in fact T* = oo for all X € T r b e , uq G c <o 2a ). 

We define two distances: write 

OO 

dx lhe (X, X) = 2 -m (l A dx rbe (m) (X, X)). (32) 

771=1 


For u G fFrbe and c > 0 we define 

t c (u) = inf{t ^ 0 : |M|® rbe(t) ^ c} 
and then for u G @ r be,x, u € %> rhe f 

drbe c(w, it) = 111/ — « , || + ||'U** — it** 11 

’ V ' 11 11 cAT C (t.)AT C (a) 11 c;At c (u)At c (u) 

and 

OO 

d r be(w, it) = 2 _m (l A (■drbe,m('U,«) + |r m («) - r m (it)|)). (33) 

771=1 

Theorem 3.7. Lei (X n ) n ,X G -T r b e fre such that lim n dx rbe (X n ,X) = 0 and let (uq),uo G 
^ 2a 6 e such that lim n ||«Q — tto||a= 0. Denote the solution to the rough Burgers equa¬ 
tion (27) driven by X n and started in X n (0) + Uq by u n , and write u for the solution 
driven by X and started in X(0) + no- Then 

lim d rhe (u n . u) = 0. 

71 

Proof. Let in G N and define M = max{||X||* rbe ( m ), ||X n ||^ rbe(m ), ||u 0 || 2 a, IK II 2 a}- Using 
the same arguments as in the proof of Theorem 3.5, we see that there exists a constant 
C > 0, depending only on m and M , such that 

drbe,m(^A ^n) ^ C(^A’ r b e (m) (^> ^n) T ||ri() ^0 11 2 o) • 

It remains to show that r m (u n ) converges to T m (u ) for all m. Write T = T m (u). Using 
the same arguments as in Lemma 3.4, we obtain that z n = u** — uh, y n = vft — Un solves 
an equation of the type 

z n = V (y 2 ) + A n (y n ) + B n (z n ) + Vh n , z n ( 0) = uq — Uq, 

where (A n ) and (. B n ) are sequences of bounded operators from 2 *?^ to CT^ 2a ~ 2 and from 
t° CtT> 2ol - 2 respectively, with operator norms bounded uniformly in n, and where 
( I(Vh n )) converges to zero in 2 z?^ a . Furthermore, we have y n = 2y n -4< Q n + z n + e(ra), 
where e(n) converges to zero in . Since also z n ( 0) converges to zero in ^ 2a , we can 
use the same arguments as above to get that (y n ,z n ) converge to zero in Jzff! x 2 z?^“, 
and therefore lim n T rn {u n ) = T. □ 
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4 KPZ and Rough Heat Equation 

We discuss the relations of the RBE with the KPZ and the multiplicative heat equation 
in the paracontrolled setting. 

4.1 The KPZ equation 

For the KPZ equation 

Jzfh = (D/i) 2 + £, h(0) = ho, 

one can proceed, at least formally, as for the SBE equation by introducing a series of 
driving terms (y T ) Tg -j defined recursively by 

jgfy* = c, c£ Y ^T2) = (Dy ri Dy r2 ) 

and such that Dy r = X T for all r G T; we will usually write Y = Y*. However, 
now we have to be more careful, because it will turn out that for the white noise £ 
some of the terms Y T can only be constructed after a suitable renormalization. Let us 
consider for example y v . If <p: M —* R is a smooth compactly supported function with 
<£>(0) = 1 and £ £ = <^(eD)^, then as £ 0 we only obtain the convergence of Y £ v after 

subtracting suitable diverging constants (which are deterministic): there exist (c £ ) £ >o, 
with lim £ _ ) .o c £ = oo, such that (y £ v (t) — c £ t)t ^>o converges to a limit y v . We stress the 
fact that while c £ depends on the specific choice of ip, the limit y v does not. To make 
the constant c £ appear in the equation, we should solve for 

Jz?/t £ = (D h £ ) 2 - c £ + £ e , h e ( 0) = p(eD)h 0 . 

In that sense, the limiting object h will actually solve 

2 &h = (Dh)° 2 + h(0) = h 0 , (34) 

where (D/i)° 2 = (D/i) 2 — oo denotes a renormalized product. The reason why we did not 
see this renormalization appear in Burgers equation is that there we differentiated after 
the multiplication, considering for example (t) = D(y v (f) — c e t) = Dy £ v (f). There 
are two other terms which need to be renormalized, y v and Y w . But on the level of 
the equation these renormalizations cancel exactly, so that it suffices to introduce the 
renormalization of yY by subtracting c £ . We still consider a fixed a € (1/3,1/2). 

Definition 4.1. (KPZ-enhancement) Let 

Tk PZ cjz^x jy 2 “ x jy Q+1 x jy 2a+1 x jy 2 “ +1 x cv 2 *- 1 

be the closure of the image of the map 

e kpz :2s?c£/ 2 (R,c°°(T)) xixK^y kpz , 

given by 

e kpz ( 0 ,c v ,c xy ) = (y( 0 ),y v ( 0 ),y x '( 0 ),y^( 0 ),y v ( 0 ),(DPoDy)( 0 )), (35) 
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where 


SfY(0) = 9, 

JSf Y v (9) = (D Y(9)) 2 - c v , 

&Y\0) = DY(0)DY v (0), 

JfY^(0) = DY^ o BY + c v /4 

JfY v (9) = BY v (9)BY v (9)-c v , 

JYP(O) = D Y{6), 

all with zero initial conditions except Y(9) for which we choose the “stationary” initial 
condition 

Y(9)(0)= [ P_ s n^ o 0( S )ds. 

Here and in the following n^o denotes the projection on the (spatial) Fourier modes 7 ^ 0. 
We call @kpz (# 5 c v , c^) the renormalized KPZ- enhancement of the driving distribution 
9. For T > 0 we define Tkpz)^) = J4pz|[o,T] and we write |^||y kpz (T) for the norm of ¥ 
in the Banach space JOjt x x 1 x JZ ^ Q+1 x 2 z ?^" +1 x Ct c ^ 2oi ~ 1 ■ Moreover, we 
define the distance dy kpz(T) (Y, ¥) = ||¥ - Y||y kpz(T) . 

YV 

Observe the link between the renormalizations for Y v and Y . It is chosen so that 
c v will never appear in the equation. Also note that for every ¥ £ 34pz there exists an 
associated X £ Aj.be> obtained by differentiating the first five entries of Y and keeping 
the sixth entry fixed. Abusing notation, we write X = DY. 

The paracontrolled ansatz for KPZ is h € f^k pz = ^k pz ,Y if 

h = Y + Y v + 2Y* + h p , h p = h'^P + h^ (37) 

where h! £ 2 z?“, and h' £ «£? 2q,+1 . For h £ £Pk pz it now follows from the same arguments 
as in the case of Burgers equation that the term 

(D hf 2 = (D h) 2 - c v 

is well defined as a bounded multilinear operator. Now we obtain the local existence 
and uniqueness of paracontrolled solutions to the KPZ equation (34) exactly as for the 
Burgers equation, at least as long as the initial condition is of the form K(0) + h( 0) with 
h( 0 ) £ ^ 2q+1 . 

Moreover, let X = DY £ A r b e , and let h £ fFkpz,Y- Then Bh = A + A V + 2A^ + D h p , 
where 

Bh p = (Bh') -« P + h' -« DP + D/i # = h' -« Q + (D/i) # 

with (D h)^ € 2z? 2a_£ for all e > 0. This follows from the fact that c C'( 2a_£ )/ 2 ^ ,1 + £ . 

In other words Bh £ * be x for all 5 < a (using an interpolation argument one can actu¬ 
ally show that (D h)$ £ Jz? 2a , but we will not need this). Moreover, if h £ fj^pz solves the 
KPZ equation with initial condition ho, then u = Bh € solves the Burgers equation 
with initial condition Bh$. 
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Conversely, let u £ IZrhe- It is easy to see that u = Dh + /, where h € ^k pz and 
/ £ CV 2a . Therefore, the renormalized product u° 2 = u 2 — c v is well defined. Note that 
u° 2 does not depend on the decomposition u = D h + f . The linear equation 

2S?h = vZ 2 + 9, h( 0) = Y(0) + h 0 , 

has a unique global in time solution h for all ho € c ^ ?2a+1 . Setting h p = h —Y —Y v —2Y^, 
we get 

2Yh p = 2S?(h - Y - Y v - 2Y X ') = «° 2 - ((DY) 2 - c v ) - 2DYDY V . 

Recalling that u° 2 = (D h)° 2 + 2/Dh + f 2 = ((DY) 2 - c v ) + 2DYDY V + J^(^“ +1 ), we 
deduce that h p £ 2z? a+1 . Furthermore, if we set h! = 2 vP + and expand u° 2 , then 

2z?(h p -h’~«P) = u° 2 - ((DY) 2 - c v ) - 2DYDY V 

- («Sf (h' -« P) - h’ -« Z£P) + h' -« T£P 
= (h' y X - h' ^ X) + .i*(2Y 2a+1 ). 

Using Lemma 2.8, we get that hf £ Jz? 2 " +1 an q therefore h £ fFkpz- By construction, D h 
solves the equation 2z?(D h) = Du 2 + Dh with initial condition Dh(0) = X(0) + Dho- So if 
u solves the Burgers equation with initial condition X(0)+Dho, then Jz?(D h — u) = 0 and 
(Dh — u)(0) = 0, and therefore Dh = u. As a consequence, h solves the KPZ equation 
with initial condition y(0) + ho- 

We can also introduce analogous distances dy kpz (Y, Y) and dkpz(h, h) as in Section 3.4, 
equations (32) and (33). We omit the details since the definition is exactly the same, 
replacing every occurrence of “rbe” by “kpz” and adapting the regularities to the KPZ 
setting. 

Theorem 4.2. For every (Y, ho) £ Xkpz x c £ 2apl there exists T* £ (0, oo] such that for 
all T < T* there is a unique solution (h, h',hf) £ £Pk pz (T) to the KPZ equation 

££h = (Dhf 2 + e, h( 0) = y(0) + ho, (38) 

on the interval [0, T], The map that sends (Y, ho) £ 34pz x < t? 2a+1 to the solution 
h € ^kpz,Y Is continuous with respect to the distances (dy kpz + IHha+i, ^kpz); and we can 
choose 

T* = inf{i ^ 0 : ||h||^ pz(t) = oo}. 

If X = DY, then X £ X T ^ e and u = Dh € f^-beCO solves the rough Burgers equation 
with data X and initial condition it(0) = X(0) + Dho- Conversely, given the solution u 
to Burgers equation driven by X and started in X(0) + Dho, the solution h £ ^k pz (T) to 
the linear PDE 2z?h = u 02 + 6, h( 0) = y(0) + ho, solves the KPZ equation driven by Y 
and with initial condition Y (0) + ho. 
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4.2 The Rough Heat Equation 

Let us now solve 

J£w = wo!;, w(0) = Wo, (39) 

where o denotes again a suitably renormalized product and where £ is the space-time 
white noise. We replace for the moment £ by some 9 £ C(R, C' 0O (T)) and let ¥ = 
@k pz (0, c v , c*) for some c v , c ^ € R be as in Definition 4.1. Then we know from the Cole- 
Hopf transform that the unique classical solution to the heat equation Jzfrc = w(0 — c v ), 
h( 0) = e h ° is given by w = exp (h), where h solves the KPZ equation with initial condition 
ho- Indeed, we have 

SYw = exp(h)Jfh — exp(/i)(D/i) 2 = exp(/i)((Dh) 2 — c v + 9) — exp(/i)(D/i) 2 = w{9 — c v ). 

Extending the Cole-Hopf transform to the paracontrolled setting, we suspect that the 
solution will be of the form w = e Y+Yw+2 ' i *+h p where h p = h! P + h!& with h! £ Jf 1 “ 
and fif € Jz? 2a+1 , or in other words 

w = e 5 +1 V+2Y ' < 'iu p , w p = w' + w^, 

with w' £ j£?“ and € J5f 2a+1 . So we call w paracontrolled and write ic £ fF r h e ,Y = ^rhe 
if it is of this form. 

Our first objective is to make sense of the term wo9 for w £ fF r h e - In the renormalized 
smooth case Y = ©k pz ($, c v , c^) for 9 £ C(R+, C' 00 (T)) and c v ,c^ £ 1 we have for 

W £ S^rhe 

&w = e y+yV+2yX ' + 2 Sf(Y + Y v + 2 Y*)w p - (DY + DY V + 2D Y*fw p 

- 2e y+yV+2yX D(y + y v + 2Y*)Bw p . 

Taking into account that 2 YY = 9, 2z?Y v = (DY) 2 — c v , and SYY* = DYDY V , we get 

w{9 - c v ) = &w- e y+yV+2yX ' [ - [4(2YY^ + Dy ^ ^ BY + BY* y DY) + &Y V ]w p 

+ [4DY V DY^ + (2D Y*) 2 ]w p 
+ ^w p -2B{Y+ Y V + 2Y*)Bw p }. (40) 

At this point we can simply take the right hand side as the definition of wo6 = w(9 — c v ). 
For smooth 9 and Y = 0 r b e ($, c v , c*), we have just seen that wo 9 is nothing but the 
renormalized pointwise product w(9 — c v ). Moreover, the operation wo6 is continuous: 

Lemma 4.3. The distribution w o 9 is well defined for all Y = © r b e ($, c v , c^) with 
9 € JY CffJ 2 (R, C 00 (T)), c v ,c V £ R, w' £ and w B £ 2z? 2q . Moreover, wo 9 depends 

continuously on (Y ,w',w^) £ 34pz x 2z?“ x Jz? 2a . 
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Proof. Consider (40). Inside the big square bracket on the right hand side, there are 

three terms which are not immediately well defined: DYDie p , 2Y > y'^ s w p , and 2z?Y^ w p . 

For the first one, we can use the fact that w p = w' P + vr and that ¥ “contains” 

v p 

DFoDF, just as we did when solving the KPZ equation. The product A Y v w p is 
well defined because AY^ G CV 2 “ _1 and w p G C ( to a . The product dtY^w p can be 
defined using Young integration: we have Y^ G C^ 1 ^ 2 L°° and w p G C^j 2 L°° and 
3cc/2 + 1/2 > 1. The same arguments show that 2z?Y^ w p is well defined. 

So taking into account that w p = w' -4< P + and using Lemma 2.8 as well as 
Bony’s commutator estimate Lemma 2.6, we see that the term in the square brackets is 
of the form 

dfiiv* + Y^ + Y V ) + {w' - 4w p DY* - 2Bw p ) P BY + C^ 2a ~ l . 

The second term is paracontrolled by DY, and since e* +2i is paracontrolled by Y, 
the product between the two is well defined as long as Y o DY is given. But of course 
Y o DY = D(Y o Y)/2. Finally, the product e 1 V+25 ^ + Y^ + Y^) can be defined 

using Young integration as described above. □ 

Besides it agreeing with the renormalized pointwise product in the smooth case and 
being continuous in suitable paracontrolled spaces, there is another indication that our 
definition of wo 9 is, despite its unusual appearance, quite useful. 

Proposition 4.4. Let ^ be a space-time white noise onlx T, and let o be its 

natural filtration, 

Tfi = a (£(1( S1]S2 ]^ : —oo < s\ < s 2 ^ t, ip G C' 00 (T,M)) , tfi 0, 

and is its right-continuous completion. Let ¥ G Ykpz be constructed from f as 

described in Theorem 9.3, and let (wt)t ^o be an adapted process with values in 2zf Q+1 , 
such that almost surely w G ^kpz,¥ with adapted w' and . Let if : M+ x T —x R be a 
smooth test function of compact support. Then 

/ ip(s, x){w o £)(s, x)dxds = / x)w(s, x))£(s, x)dxds, 

Jr+xT JR+xT 

where the left hand side is to be interpreted using our definition of w o £, and the right 
hand side denotes the ltd integral. 

Proof. Let <p be a smooth compactly supported even function with <^(0) = 1 and set 
= <^(eD)£. Theorem 9.3 states that for elf = ip 2 [x)Ax and an appropriate choice 

of cf, we have that ©k P z(£e> c^, c^) converges in probability to ¥ as e —X 0. Moreover, 
by the stochastic Fubini theorem ([DPZ14], Theorem 4.33) and the continuity properties 
of the Ito integral, 

/ (^(s, x)w(s, x))£(s, x)daxls = lim / (fi>(s,x)w £ (s,x))^ £ (s,x)dxds, 

J R + xT £ ^°Jr + xT 
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where we set 


w £ = e Y+y Y +^Y^-^p e -\-vfi), 

and where the YJ and P £ are constructed from p, , and c Y as in Definition 4.1. But 
by definition of ¥ and the continuity properties of our product operator, also w £ o 
converges in probability to w o £ in the sense of distributions on R + x T. It therefore 
suffices to show that 


t rt rt 

(w £ o£ >£ )(s,x)ds = / w e (s,x)^ e (s,x)ds = / w £ (s,x)d s B £ (s,x) 

Jo Jo 


for all (t,x) € M+ x T, where we wrote B e (t,x ) = £ £ (s,x)ds which, as a function 

of t. is a Brownian motion with covariance (27 t) _1 ip(ek) 2 = 2+ o(l) for every 
i€l. The first term on the right hand side of (40) is 2z fw £ . Let us therefore apply Ito’s 
formula: Writing w Fe = w' P £ + , we have for fixed x € T 


d t w £ 


Aw £ dt = e 


y c +yV+2y} 



Aw Ps dt 


+ w p * (dt(Y e + W v + 2Y}) - A(Y e + K V + 2Y*)dt) 

- (\)Y e + DL V + 2mf) 2 w Pe dt - 2D(K + Y? + ‘2Y^)Y)w p dt 

+ ^ e d(y £ (•,*))*, 


where we used that all terms except Y e have zero quadratic variation. But now d(y e (-, x))t = 
d (B s (-,x))t = (2+ o(l))df, and therefore the Ito correction term cancels in the limit 
with the term — w e c£dt that we get from w £ (d t Y^ — AY^dt — (Dy e ) 2 dt). So if we 
compare the expression we obtained with the right hand side of (40), we see that 


w £ o^ £ dt = w £ (d t Y £ - AWdf) + o(l) = w £ d t B £ + o(l), 


and the proof is complete. □ 

With our definition (40) of w o 9, the function w € ^ r he solves the rough heat 
equation (39) if and only if 


e y+y v +2Y x ' 


i nv'f mv i nv^ v nw i i (I'ivVnv'f' i /onv'f'^L, P 


[4 (AfY* + Dy v -< Dy + Dy v >- Dy) + A?y v + 4Dy v Dy v + (2Dy' 


w 


+ A?w p - 2D(y + y v + 2Y*)Dw ] 


= 0. 


v Xf 

Since e Y+ * +21 is a strictly positive continuous function, this is only possible if the 
term in the square brackets is constantly equal to zero, that is if 


2z? w P 


[4 (jfY^ + Dy* ^ d y + Dy* >- Dy) + jgfy’ v + 4Dy v Dy* + (2Dy*) 2 ]wP 
+ 2 D(y+ y v + 2y x ')Dw p . (41) 
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We deduce the following equation for uA 

Jgf w # = (2D w p + 4w p DY* - w') -< BY + (4+ 5fY V ')w p + jg?jgf 2 “ +1 , 

where we used again Bony’s commutator estimate Lemma 2.6. Thus, we should choose 
w' = 2D w p + 4D Y^ f w p . Now we can use similar arguments as for Burgers and KPZ 
equation to obtain the existence and uniqueness of solutions to the rough heat equation 
for every ¥ £ Ykpz with initial condition e Y ° wq with wq £ e Y 2a . Since the equation is 
linear, here we even obtain the global existence and uniqueness of solutions. Moreover, 
as limit of nonnegative functions w is nonnegative whenever the initial condition ivo is 
nonnegative. 

Theorem 4.5. For all (¥, wo) £ Y^ pz x c Y 2a+1 and all T > 0 there is a unique solution 
(w,w',w^) € ^kp Z (T) to the rough heat equation 

JYw = wod , w(0) = e Y (°)wo, (42) 

on the interval [0, T]. The solution is nonnegative whenever wq is nonnegative. IfY is 
sampled from the space-time white noise £ as described in Theorem 9.3, then w is almost 
surely equal to the ltd solution of (42). 

4.3 Relation between KPZ and the Rough Heat Equation 

From KPZ to heat equation. Let now h £ i^k pz and set w = e h . Then w = 

eJ V+21 ^w p with w p = e hP . To see that w £ 3> r \, e , we need the following lemma: 

Lemma 4.6. Let F £ C 2 (M,M), h! £ P £ ^" +1 with T£P £ Ctf a ~ l , and 0 £ 

J 2 ? 2 “ +1 . Write h p = h! -4< P + 0 £ »£f Q+1 . Then we have for all e > 0 

F{h p ) - {F’{h p )h') ^P £ ^' 2q + 1 - £ . 

Proof. It is easy to see that the difference is in C"^ 2 “ +1 , the difficulty is to get the right 
temporal regularity. If we apply the heat operator to the difference and use Lemma 2.8 
together with Lemma 2.6, we get 

_Sf( F(h p ) - ( F'{h p )h') -«P) = F' (h p )JY(h' -4! P + /i # ) - F"(h p ){Bh p ) 2 

- ( F'{h p )h') ^TYP + C^ 2a ~ l 
= F'(h p )d t h* + Ctf 2a - 1 . 

So an application of the Schauder estimates for SY shows that 

(F(h p ) - (F\h p )ti) -« P - f P._ s {F'(h p (s))d s /i # (s)}ds £ Jz? 2a+1 . 

Jo 

Since F'(h p ) £ C l ( a+1 )' /2 L°° and hf £ C^ 2a+l ^ 2 L°° , the Young integral 

I(t) = f F'(h p (s))dh\s) 

Jo 
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is well defined and in (7( 2 «+ 1 )/ 2 L°°. We can therefore apply Theorem 1 of [GLT06] to 
see that f Q P.- S dl(s ) G Cd 2a+1 ~ £ ^ 2 L°° for all e > 0, and the proof is complete. □ 

As a consequence of this lemma we have e h G for every h G @kpz and 6 < a, 
with derivative e h h! G 2zf". 

Lemma 4.7. Let ¥ G 34pz> T > 0, and let h G S>kpz,Y solve the KPZ equation on [0, T] 
with initial condition h( 0) = ho, and set w = e h G y derivative e h h!. Then w 
solves the rough heat equation on [0,T] with initial condition ic(0) = e h °. 

Proof. Recall that 

L£h p = -Y - Y w - 2Y y ') = (D hf 2 - ((DT) 2 - c v ) - 2DTDy v . 

From the chain rule for the heat operator we obtain 
2z?u; p = w p ^h p - w p (Bh p ) 2 = w p (( D/if 2 - ((DT) 2 - c v ) - 2DTDy v - (D/i p ) 2 ). 
But now 

(D hf 2 - ((DT) 2 - c V ) - 2DTDy v - (D h p ) 2 

= 4 (JfY^ + DT* -< D Y + >- Y) + F V + 4DF X 'DF V + (2D Y^) 2 

+ 2D(y + Y v + 2F X ')D/i p , 

and D h p = D(logic p ) = D w p /w p , from where we deduce that 

5?w p = [4(^y^ + Dy^ -< df + by* y y) + y v + aby^by^ + ( 2 Dy x ') 2 ]w p 
+ 2D(y + y v + 2y' < ')Bw p . 

In other words, w p satisfies (41) and w solves the rough heat equation. □ 

From the heat equation to KPZ. Conversely, let w = e Y + YV + 2Y ^ w p G £^ r h e be 
strictly positive and define h = logic. Then h = Y + Y v + 2y^ + log w p , and Lemma 4.6 
yields log w p = (ic'/ic p ) -4< P + 22 ? 2q+1_£ for all s > 0, in other words h is paracontrolled 
with derivative w'/w p . 

Lemma 4.8. Let ¥ G 34pz, T > 0, and let w G fF r he,Y be a strictly positive solution to 
the rough heat equation with initial condition w( 0) = wo > 0 on the time interval [0, T]. 
Set h = logic G with derivative w'/w p . Then h solves the KPZ equation on [0,T] 

with initial condition ho = log wo ■ 

Proof. We have 

&h p = 2^(logic p ) = \^w p + -i—(Din p ) 2 

ur (uv ) z 
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= (JzfF^ + DF* -< DF + D Y* y DF) + JZT V + 4DF V DF* + (2D Y*) 2 ]v? 

+ -^2D (f + f v + 2y x ')D^ p + (d/i p ) 2 , 

ur 

where for the last term we used that Die p = w p Dh p . Using this relation once more, we 
arrive at 

jz? h p = 4(2sfy^ + df v x Dy + Dy v >- df) + spy' v + 4 df v df* + (2DF*) 2 
+ 2D(y + y v + 2F X ')D/i p + (D/r p ) 2 , 

or in other words 2z?/r = (D/i)° 2 + 6 1 . □ 

Let us summarize our observations: 

Theorem 4.9. Let Y,wo,h be as in Theorem 4-5 and set 

T* = inf{f > 0 : minu;(t,x) = 0}. 

Then for all T < T* the function log iu|[o,t]£ ^kpz(T) solves the KPZ equation (38) 
driven by Y and started in F(0) + log wq . Conversely, if h £ @k pz (T) solves the KPZ 
equation driven by Y and started in Y(0) + ho, then w = exp (h) solves (42) with wq = 
exp (h 0 ). 

As an immediate consequence we obtain a better characterization of the blow up 
time for the solution to the KPZ equation: 

Corollary 4.10. In the context of Theorem 4-2, the explosion time T* of the paracon- 
trolled norm of the solution h to the KPZ equation is given by 

T* = inf{f ^ 0 : minexp(/i(t, x)) = 0} = infjf ^ 0 : ||/i(t)||L°°= oo}. 

zST 

5 Interpretation of the nonlinearity and energy solutions 

The purpose of this section is to give a more natural interpretation of the nonlinearity 
Drt 2 that appears in the formulation of Burgers equation. We will show that if f is a 
space-time white noise, then u is the only stochastic process for which there exists u' 
with (it, u') £ S> v he almost surely and such that almost surely 

dfU = A u + lim D((/?(eD)u) 2 + D£, u(0) = uq , 

£—>0 

whenever ip is an even smooth function of compact support with <p(0) = 1, where the 
convergence holds in probability in ^(R + x T). 

This is a simple consequence of the following theorem. 
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Theorem 5.1. Let X G T r b e and 99 € C'°°(M) with compact support and <^(0) = 0. 
Define for e > 0 the mollification Xf = ip(eT))X T , Q £ = <p(eD)Q, and then 

X £ = (X £ ,I(D[(X £ ) 2 ]),I(D[X^X e }),I(D[x}oX £ ]),I(D[(X^) 2 ]),Q £ oX £ ). 

Let T > 0 and assume that ||X £ — X||^ rbe ( T ) converges to 0 as e ^ 0. Let u G f^ r b e (T). 
Then 

lim/(D[(^(eD)«) 2 ]) = /(Dn 2 ), (43) 

£—>■0 

where Du 2 on t/ie right hand side denotes the paracontrolled square and where the con¬ 
vergence takes place in C ([0 ,T],&). 

Proof. The right hand side of (43) is defined as 
A v + 2X y ' + A V + 4X^ + I (2D(u q X) + 2D(A v (n Q + 2A X ')) + D((n Q + 2X X ') 2 )) , 
where 

u Q X = uPX + uyX + u t oX + {(u'^Q)-u , pQ)oX + C(u',Q 1 X) + u'{QoX). (44) 

On the other side we have <p{eD)u = (X e + X)f + 2X^ + <p(e'D)ifi), and our assumptions 
are chosen exactly so that the convergence of all terms of 7(D[(<^(eD)u) 2 ]) is trivial, 
except that of /(D[(<^(eD)(it / -4< Q )) o X £ ]) to 7(D[(u' -« Q) o A]), where the second reso¬ 
nant product is interpreted in the paracontrolled sense. Since / i->- 7(D/) is a continuous 
operation on C([0, T], &), it suffices to show that {ip{eD)(u' ~¥.Q))o X e converges in 
C{\f),T],@') to {v! -^.Q)oX. We decompose 

(ip(eD)u Q )o X e = (<p(sD)[(u' ~«Q)- {v! <Q)]) o X e + {ip{eD)(v! <Q)) ° X £ 

and use the continuity properties of the resonant term in combination with Lemma 2.8 
to conclude that the first term on the right hand side converges to its “without £ coun¬ 
terpart”. It remains to treat {ip(eYf){v! A Q)) o X e . But now Lemma 5.3.20 of [Per 14] 
states that 

||<^(eD)('?/ -<Q) — v! -< Q £ \\c T <^2 a -s< e & \\u'\\c T ^\\Q\\ c t ^ 

whenever <5^1. If we choose 5 > 0 small enough so that 3a — 5 > 1, this allows us to 
replace ((p(eT))(u' -< Q )) o X £ by {v! -< Q e ) o X e . Then we get 

{v! A Q e ) o X e = C{u\ Q e , Xe) + u'(Q e o X e ), 

and since C is continuous and by assumption Q £ °X £ converges to Q oX in Ct < ^ ?2 "~ 1 , 
this completes the proof. □ 
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Corollary 5.2. Let f be a space-time white noise onRx T, let X G A r b e be constructed 
from £ as described in Theorem 9.3, and let uq G YI 2ol almost surely. Assume that u is 
a stochastic process for which there exists v! with (u,u') G ^ r b e almost surely and such 
that almost surely 

dtu = A u + lim D((/?(eD)u) 2 + D£, w(0) = uq + Y(0), 

£- S >0 

whenever p is an even smooth function of compact support with </?(0) = 1, where the 
convergence holds in probability in x T). Then u is almost surely equal to the 

unique paracontrolled solution to the rough Burgers equation driven by X and started in 
uq + Y (0). 

Proof. Theorem 5.1 shows that almost surely 

ACu = Drt 2 + D£, u( 0) = uq + Y (0), 

where the square Du 2 is interpreted using the paracontrolled structure of (u,u'). This 
implies that also ( u , 2 u® + 4X^) is almost surely paracontrolled, where uQ = u — X — 
X v — 2X^ . However, another application of Theorem 5.1 shows that almost surely the 
paracontrolled square Du 2 is the same for {u,v!) and for (u, 2v9, 4X^), and therefore u 
is the paracontrolled solution of the equation. □ 

Remark 5.3. There have been several efforts to find formidations of the Burgers equa¬ 
tion based on interpretations like (43), see for example [Ass02,GJ14,Assl3,CO14,GJ13]. 
In [GJ14] the concept of energy solutions is introduced. A process u is an energy solution 
to the stationary Burgers equation if u(t,-) is (1 /\[2 times) a space white noise for all 
times t ^ 0, if it satisfies some “energy estimates”, and if for any (p as in Theorem 5.1 
we have 

u(t ) = P t u 0 + hm/(D[(^(eD)u) 2 ])(f) +I(f)(t) 

£—>■0 

almost surely for all t ^ 0, where the convergence takes place in L 2 in 3>'{T). Li [GJ14] 
it was also shown that energy solutions are a useful tool for studying the convergence of 
rescaled particle systems to the stochastic Burgers equation: it can often be shown that a 
given rescaled particle system is tight, and that any limit point in distribution is an energy 
solution of Burgers equation. However, so far it is not known whether energy solutions 
are unique and this appears to be a difficult question. On the other side, Corollary 5.2 
shows that there exists at most one energy solution which is almost surely paracontrolled. 

Remark 5.4. For simplicity we formulated the result under the assumption u( 0) = 
Y (0) + uq for some relatively regular perturbation uq. Using the techniques of Section 6 
it is possible to extend this to general u(0) G < ^’~ 1+£ . 

6 Singular initial conditions 

In this section we extend the results of Section 3 to allow for more general initial con¬ 
ditions. For that purpose we adapt the definition of paracontrolled distributions. To 
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allow for initial conditions that are not necessarily regular perturbations of X( 0 ), we 
will introduce weighted norms which permit a possible explosion of the paracontrolled 
norm at 0. Moreover, as was kindly pointed out to us by Khalil Chouk, in the case of 
the rough heat equation it is very convenient to work on 

OO 

for a £ R, rather than restricting ourselves to p = oo. For example the Dirac delta has 
regularity < io^ v 1 on T, so by working in spaces we will be able to start the rough 
heat equation from (—A)^do for any (3 < 1/4. In fact for the same reasons and in the 
setting of regularity structures an approach based on Besov spaces with finite index was 
developed recently in [HL15]. 

For p £ [l,oo] and 7 ^ 0 we define A / O r L p = {u: [0, T] —> £t'(T) : \\v\\ m i tLp < 00 }, 
where 

IMI MlLP= SUp {||t 7 u(t)|| L v}. 
te[o,T] 

In particular 11 ^ 11 ^ 7 ^ = ||u||c T i> for 7 = 0 . If further a € (0,2) and T > 0 we define 
the norm 

ll/ll^7’“(T)= nrax{||t i-> t 7 /(t)|| c «/ 2 LP , \\f\\ M l^} 
and the space 2 ^ 7 ’“(T) = {/: [ 0 ,T] —> $)' : \\f\\^’ a (r)< 00 } as well as 
= {/:R + ->• & : /| [0 , T] G if 7 ’ a (T) for all T > 0}. 

In particular, we have 2 z?^“(T) = 

6.1 Preliminary estimates 

Here we translate the results of Section 2.3 to the setting of Ftfp’ a spaces. The estimates 
involving only one fixed time remain essentially unchanged, but we carefully have to 
revisit every estimate that involves the modified paraproduct -4< , as well as the Schauder 
estimates for the Laplacian. 

Let us start with the paraproduct estimates: 

Lemma 6.1. For any (3 £ R and p £ [l,oo] we have 

Wf<9\\vP~P min {ll/IU-lbll^, II/IM 0 II/ 9 }, 

and for a < 0 furthermore 

\\f ^g\\^+P< a ,f) min{||/||^a \\g\\/3, WfWaMtfP}- 
For a + (3 > 0 we have 

\\f 0 9\\^+P< a ,t3 min{||/||^a|| 5 ||^, \\f\\ a \\g\\^}. 
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To adapt the multiplication theorem of paracontrolled distributions to our setting, 
we first adapt the meta-definition of a paracontrolled distribution: 

Definition 6.2. Let (5 > 0, a € R, and p £ [l,oo]. A distribution / £ ’tfff is called 
paracontrolled by a £ c € a and we write / £ @P(u), if there exists f £ such that 

Let us stress that u £ is not a typo, we do not weaken the integrability assump¬ 
tions on the reference distribution u. 

Theorem 6.3. Let a,/3 € (1/3,1/2). Let u £ c £ OL , v £ c £ a ~ 1 , and let (/,/') € 
and ( g,g') £ 3>p(v). Assume that uov £ ^ 2a_1 is given as iimit of (u n ov n ) in c £ 2a ~ 1 , 
where ( u n ) and (v n ) are sequences of smooth functions that converge to u in c € a and to 
v in c £ a ~ x respectively. Then fg is well defined and satisfies 

Wfg-f ~<9\\^-i< (||/ , ||/3||«||a + ||/ tt ||a+j8)(l|y , ||^||v||a-l + l|5 tt ||^+P-i) + ll/Vll^ll« 0 v||2a-l- 

Furthermore, the product is locally Lipschitz continuous: Let u £ c £ a , v £ c £ a ~ 1 with 
uov £ and let (/,/') € and (g,g') £ 3>p(v). Assume that M > 0 is an 

upper bound for the norms of all distributions under consideration. Then 

II (fg - f<g)~ (fg - f~<g)\\^- 1 ~ i 1 + m3 ) II/' - /'Ik+lk' - <?'ll<^+lk - «IU 

+ Ik - t’lla-l + ll/* 1 - / # ||a+/3 
+ Ik" - ^ll^a+zS-ilk 0 V - n°v|| 2 a-l • 

// /' = /' = 0 or g' = g' = 0, f/ien Af 3 can 6 e replaced by M 2 . 

In this setting, the proof is a straightforward adaptation of the arguments in [GIP15], 
see also [AC15]. For an extension to much more general Besov spaces see [PT15]. 

Next we get to the modified paraproduct, which we recall was defined as 

f ^<g = yZ(QiSi-if)Aig with QJ(t)= [ 2~ 2l <p(2 2l (t - s))/(s V 0)ds. 

i 

Since for / £ Mfi r L p we have in general /(0) ^ L p , this definition might lead to problems. 
Recalling that / £ 2>’“(T) if and only if 11 —>- f 7 /(t) £ it seems reasonable to replace 

/ by /(i)l t> o, so that f 0 ^ (p(2 2l (t — s))/(sV0)ds = 0. So in what follows we shall always 
silently perform this replacement when evaluating -« on elements of J£p' a (T) or Ai^L p . 

Lemma 6.4. For any (3 £ R, p £ [l,oo], and 7 £ [0,1) we have 

t 7 \\f g(t)\\^< mm{\\f\\ M -y LP \\g(t)\\p, \\f\\ M j L oo\\g(t)\\ v p } (45) 

for all t > 0, and for a < 0 furthermore 

PWf-Kgm^Z min {ll/IU/^“lkk)ll/3, Wf\\MivAW)\\^p}- (46) 
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Lemma 6.5. Let a € (0, 2), /3 € M, p £ [1, oo], and 7 £ [0,1). Then 
P\\{f ^ 9 ~ f ^ 9){t)\\^+P< \\f\\j?2’ a (t)\\9(t)\\/3 
for all t > 0, as well as 

tf \\(S?(f~«g) - /-« (Jfg)) (t)\\ v a+p -2 < \\f\\^{t)h{t)\\h- 

These two lemmas are not very difficult to show, but at least the proof for the second 
one is slightly technical; the proofs can be found in Appendix A. 

Recall the definition of the operator //(f) = fg Pt- S f(s)ds. In the singular case we 
can adapt the Schauder estimates for I as follows: 

Lemma 6.6 (Schauder estimates). Let a £ (0, 2), p £ [1, oo], and 7 £ [0,1). Then 

WtfWyf ’“(T)~ WfWMftfp- 2 ( 47 ) 

for all T > 0. If further (3 ^ —a, then 


||s e£ P s Uo\\^(p+ a )/ 2 , a ^< 11 Mo 11^-/3 • 
For all a £ M, 7 £ [0,1), and T > 0 we have 


(48) 






(49) 


To a large extent, the proof is contained in [GIP15]. We indicate in Appendix A how 
to adapt the arguments therein to obtain the estimates above. 

Just as in the “non-explosive” setting, the Schauder estimates allow us to bound 
/ -« g in the parabolic space 2i?p ,a (T): 

Lemma 6.7. Let a £ (0,2), p £ [l,oo], 7 £ [0,1) and 5 > 0. Let f £ 2£p’ 5 (T), 
geC T tf a , and ^geCrtf 0 - 2 . Then 


11/^S , ll.s? 7 ’“(T)~ \\f\\^’ s (T)(\\ g \\ c T& a - lr \\-% ? 9\\cT‘tf a - 2 ')- 

The proof is completely analogous to the one of Lemma 2.10. Finally, we will need 
a lemma which allows us to pass between different S^p' a spaces, the proof of which can 
be also found in Appendix A. 

Lemma 6.8. Let a £ (0,2), 7 € (0,1), p £ [1, 00], T > 0, and let f £ Then 

ll/ll^-e/2 II/IIt? 7 ’“(T)- 


for all e £ [0, a A 27 ). 
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6.2 Burgers and KPZ equation with singular initial conditions 

Let us indicate how to modify our arguments to solve Burgers equation with initial 
condition «o G for arbitrary f3 < 1. Throughout this section we fix a G (1/3,1/2) 
and (3 G (1 — a, 2a). For £>0we write 

/? + £ 

7 e - 2 

Definition 6.9. Let X G and 5 > 1 — 2a. We define the space ^ = ^rbe'x °f 

distributions paracontrolled by X as the set of all ( u , u', vfi) G C&P x Jzf/Z 5 '' 5 x 22?c2>“ + ‘ 5 ^ Q+< ^ 
such that 

u = X + X w + 2X* + v! ^ Q + J. 

For 6 = a we will usually write For T > 0 we set @^p ,S (T) = ^beI [o,T] > 

and we define 

IM|^e*p,s (T) = lkll^ (T) +lk # ll j?w^+ 5 \ T y 

We will often use the notation ifi = v! -4< Q + ifi, and we will also write ( uvfi) G ii^be 
or u G • We call v! the derivative and v* the remainder. 

The exp in stands for “explosive”, by which we mean that the paracontrolled 
norm of u(t) is allowed to blow up as t approaches 0. Throughout this section we work 
under the following assumption: 

Assumption (T,M). Assume that 9,0 G C ,0O (T)) and uq,uq G and that u 

is the unique global in time solution to the Burgers equation 

££u = Dri 2 + D£, -u(O) = uq. (50) 

We define X = 0 r b e ($), U® = u — X — X v — 2X^ , and u' = 2 v9 + 4X^, and we set 
u t =u Q_ u '^q Similarly we define X, u , rfi , v! , ifi. Finally we assume that T , M > 0 
are such that 

niax{||uo||_s, ||tt 0 ||-/3, l|X||* rbe( D, Ifflkbe(T)} < M. 

From the paracontrolled ansatz we derive the following equation for v$: 

= Dri 2 - DA 2 - 2DA V A - 2&f(u' -« Q), u fl (0) = u 0 - X(0) - X v (0) - 2X^(0), 

where for the initial condition we used that u' -«Q( 0) = 0. Using similar arguments as 
for Lemma 3.4, we deduce the following result. 

Lemma 6.10. Under Assumption (T,M) we have 

- &x\t) || 2a _ 2 < r^(l + Af 2 )(i + ||u|||exp (T) ). 
for all t G (0, T\. If further also ||it||^exp ||u|k ex p_('n^ M, then 

rbejXv > rbe,X'' 

- 2^X V - S£X^) - (J^u# - 2^A V - &X b ))(t)\\ 2a - 2 (51) 

<t 7,3 (1 + AL 2 )(d^ rbe ( T )(X,X) + || u' - u'\\^‘’ a (T)+\\ u ^ ~ “ # ll^ 2 “’ 2 “(t))' 
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Sketch of Proof. We use the same decomposition for - ££X V - 2 Vx\{t) as in the 

proof of Lemma 3.4. The most tricky term to bound is 

l|D(u«(()) 2 iu_ 2 < Il( u °(t)) 2 lk- 1 < ll(«+))lli«< r l ”||««|| A1?I .„) 2 . 

Now 2yo = 7^ and Lemma 6.4 and Lemma 6.8 yield 

II“ Q Ha^°L~ ^ ll u Q\\M^L x +\\ l ^\\Mf?L°°< ll u ll_M 73 L°° IIQIIcyif^ + lMlI 

~ (lMl^“’“(T) + lMll^20.2“( T ))( 1 + ll X ll* rbe (T))- 

Every other term in the decomposition of — JfX^)(t) that does not 

explicitly depend on ifi can be estimated with a factor t~ la and since 270 = 7/3 and 
a < (3 we get t _7a t~ 270 for all t G [0, T\. The last remaining term is then 


\\p>(J ox(t))y a +p-2< ^wx^Wa-^ 


where we used that « + /?>! and that 2 a > /?. 


□ 


Corollary 6.11. Under Assumption (T,M), we have 

||lt||^exp( T ) < M + T 72 “ 7,3 (1 + M 2 )( 1 + |M|^exp( T )) 
If further also IMI^ b p x (Tp INI^ exp ~(T)^ M, then 


|M ~ ^11 J ^72a>2a^ + ||'i/ — v! || jfTo,«( T ) < ||«0 ~ UQ || -/3+dx ibe (T) (X, X) 

+ r 72 “- 7 ^ (1 +M 2 )( d A . rbe{T)( X,f) + ||«' - filing.,a (T) + || U # - fi“||^ 2a . 2 « (T) ). 
Proof. Lemma 6.10 and Lemma 6.6 yield 

iMlI ^ a ’ 2 “(T)~ ll tt o||-/3+||X||^ rbe ( T )+T 72 “ 7/3 (l + M 2 )( 1 + ||«||^exp( T )) 
and similarly for — t/N. It remains to control v! and vl — v!. We get 
\\ u '\\^3o a ’ a (T) ^ 2 ll' uQ ll7f ( 3 ) a ’“(T)+4||^^||^3 ) “> a (r) 

< ||?/ -+ Q||^r > «'“(T) + lMll ^ 20 . 20 ^ + 11X11^(22), 

where we used Lemma 6.8 to bound |M||.s? t 2“’ a (T)i$ ll^ll^-rao. 2 **^)- Now Lemma 6.7 and 
then once more Lemma 6.8 yield 

||II^to,<* (T ) < T 7 “-^-«||u'+<Q|| x2£- a,a ( T ) 

< T ia-yp- a ||«'||^_o^-« (T) l|X|U rbe{ T) 
<T 72 “-^||n'||^o,. (r) ||X|U rbe(T) . 

Combining this with our estimate for ||^||^, 72 o,2o™, we obtain the bound for ||w'|| 

The same arguments allow us to treat \\u' — u '\\Cl 
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The main result of this section now immediately follows. Before we state it, let us 
recall that the distance d^ rbe (X, X) was defined in (32) and define for u £ and c > 0 


Tf P («) = inf{* ^ 0 : Nla“P( t )> 4 

and then for ri £ and u £ fF e , xp ~ 

rDe ^ rbe,X 


^rbe,c( w ’ ll U ^ HjSf(2)“’“(cAr c (u)ATc(li)) _ l _ ll W ^ ^ ||_^72a. 2 “ 


(cATc(u)ATc(tl)) 


and 

oo 

<£b£(«>«) = X] 2_m ( 1 A ( d rbe,m(“^) + I^O) “ f m(u)\))- 

m= 1 

Theorem 6.12. For every (X, uo) £ X r be x 'if - ' 9 f/iere exists T* £ (0, oo] such that 
for all T < T* there is a unique solution (u,u',ifi) € Sf) xp (F) to the rough Burgers 
equation (50) on the interval [0, T], The map that sends (X, no) £ X r be x to the 
solution (u,u',w) £ 3>^e * s continuous with respect to the distances {dx Ae + IMI-/3, d r be); 
and we can choose 

T* = inf{t ^ 0 : ||^||c t ^-' 3= °°}- (52) 

Proof. The proof is essentially the same as the one for Theorem 3.5. The only difference 
is in the iteration argument. Assume that we constructed a paracontrolled solution u 
on [0, r] for some t > 0. Then t 7 “(u < ^(t) — u' Q(r)) £ c £ 2a by assumption, and 
Lemma 6.5 shows that 


r 7 “ {v! -4! Q(t) - u\t) -< Q(t)) £ ^ 2a . (53) 

The initial condition for u$ in the iteration on [r, r + 1] is given by u®{t) — u'{t) -< Q(r), 
which as we just saw is in ( if 2a , and therefore we can use the arguments of Section 3.4 to 
construct a paracontrolled solution (u, v! € Cf c <o a ~ 1 x x T£ 2a on the time interval 

[0, f] for some f > 0. Extend now u and u' from [0, t\ to [0, r+f] by setting u(t) = u(t—T ) 
for t ^ r and similarly for v!. Since on [t, r + f] the function t H > t 7 “ is infinitely 
differentiable, we obtain that (u. u ') € C T +f ( . x . Moreover, ift — u' r Q € T£ 2a 

on the interval [r, r + f ], where 

u' -« T Q(t) = f 2 2 i(p(2 2 i(t — s))Sj-\u'(s V r)dsAjQ(t). 

j J — OO 

Using the smoothness of 1 e+ t 72 “ on [r, r+f], it therefore suffices to show that (u' -+ r Q— 
v! -+ Q)|[ TiT _|_f]£ Jz? 2a . But we already saw in equation (53) above that (u'^< T Q — 
v! -4< Q)(t) £ %? 2a , and for t £ [ r, r + f] we get from the second estimate of Lemma 6.5 
that 

2z?(u'+< T Q-u'+<Q)(t) = («'4< r DI-M'^DI)(f) + r 7 ^ 2tt - 2 . 

On [r, r + f] the factor f -7 “ poses no problem and therefore it suffices to control the 
difference between the two modified paraproducts. Now the first estimate of Lemma 6.5 
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shows that ( u' -4< DA" — v! -< DA)|[ T T+f ]€ 2 , and from Lemma 2.8 it follows that 

also (v! r DA — u' -< DA)|[ TiT+f ]G C^ 2 "^ 2 , so that the Jz? 2a regularity of vf | [ T ,r+f] 
follows from the Schauder estimates for the heat flow, Lemma 2.9. Uniqueness and 
continuous dependence on the data can be handled along the same lines. 

However, so far the construction only works up to time T* = inf{t ^ 0 : |M|^ e *P(t) = 
oo}, and it remains to show that T* = T* for T* defined in (54). Assume therefore that 
IMI°°- Then we can solve the equation starting in u(T* — e) for some very 
small e for which we can perform the Picard iteration on an interval [T* — e, T* + e], and 
in particular the solution (u, u', u$) restricted to the time interval [T* — e/2, T* + e/2] is 
in C([T* - e/2, T* + e/2], ^ a ~ l ) x 2z?“( [T* - e/2, T* + e/2]) x if 2 Q ([T* - e/2, T* + e/2]) 
(with the natural interpretation for the JY s ([a,b\) spaces). Using the uniqueness of the 
solution we get a contradiction to ||u||^ex P ^*^= oo, and the proof is complete. □ 

For the KPZ equation and Y £ Tkpz we can introduce analogous spaces ^pz Y with 
distance d^(h, h ) and obtain the analogous result: 

Theorem 6.13. For every (¥, ho) € 34pz x there exists T* £ (0, oo] such that for 

all T < T* there is a unique solution (h,h',h^) £ ^pz (T) to the rough KPZ equation 

FZh = (D hf 2 + £, h{ 0) = h 0 

on the interval [0, T], The map that sends (¥, hp) £ T r b e x to the solution 

(h,h',h^) £ is continuous with respect to the distances (dy kpz + ||-||i_^,dkpz), and 

we can choose 

T* = inf{t ^ 0 : \\u\\ CtV i- P = oo}. (54) 

6.3 Heat equation with singular initial conditions 

The rough heat equation is linear, and therefore in that case it sometimes turns out to 
be advantageous to work in if® spaces for general p in order to allow for more general 
initial conditions. More precisely, we will be able to handle initial conditions wq € ffp 13 
for arbitrary (3 < 1/2 and p £ [l,oo]. Taking p = 1, this allows us for example to start 
in (—A) 7 <5 for 7 < 1/4, where 6 denotes the Dirac delta. 

We fix p £ [l,oo], a £ (1/3,1/2), and /3 £ (0,a), and we want to solve the paracon- 
trolled equation 

htfw = wof, w(0) = Wo, 
for wo £ 'if - ' 3 . For e ^ 0 let us write 

(3 + £ 


Definition 6.14. Let ¥ £ 34pz and 5 £ (1 — 2a, 1 — a — (3). We define the space 
6 = ^rh^Y °f distributions paracontrolled by Y as the set of all (u, u',u^) £ Cffp x 
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^ 5 ' 5 x ^ 1+a + 5 ’ {1+a+S) such that 

u = e 5 +5 V+2i ^(w' -4< P + w^). 

For T > 0 we set @^p ,S (T) = ^he’ S \[o,T\i an d we define 

|MU-P.*( T )= lk / ||^^ (T) + |k"||^7 1+a+5 ,(!+«+«) (t) - 

We will often write w p = w' -4< P + . 

Recall from (40) the definition of the renormalized product wo ,£ for ¥ = ©kp Z (£> c v , cV): 

woi = F£w- e y+yV+2yX ' \ - [4(J VY* + DF* -< DF + DF* >- DF) + ^Y^]w p 

+ [4DF V DF* + (2D F*) 2 ]re p 
+ FYw p - 2D(F + F v + 2Y y, )Bw p . 

It is then a simple exercise to apply the results of Section 6.1 to see that wot; depends 
continuously on (Y,w) € 3 \ pz x ^rhe 'v ■ Moreover, w € @rhe,Y S ° lveS 

h£w = WOt;, w(0) = Wq, 


if and only if 

Jgf w p = [4(2zfY^ + DF* X DF + DF* >- DF) + ^F V + 4DF V DF* + (2DF*)V P 
+ 2D(F + F v + 2F*)Dw p , 

w p ( 0) = woe-'* (°)- yV ( 0 )- 2 >>). Since wq € % 7 ff^ and v (o)- 2 y^(o) ^ cga w j ^] 1 

a > f3, the latter product is well defined and in 

We define the distance analogously to the case of the KPZ or Burgers equa¬ 
tion. Solutions to the rough heat construction can now be constructed using the same 
arguments as in the previous section, so that we end up with the following result: 

Theorem 6.15. Let S G (1 — 2a, l — a — /3). For all (Y,wq,T) £ 34pz x x [0, oo) 
there is a unique solution {w,w',w^) € ^rheY CO rou 9h h ea t equation (50) on 

the interval [0,T], The solution depends continuously on (Y,wq) with respect to ( dy kpz + 
|M|^,-/3, d^ ). Moreover, if there exist (wq) C C°° with ||u>o — Wq \\~p-. > 0 as n —>• oo and 
such that Wq ^ 0 for all n, then w(t,x ) ^ 0 for all 


7 Variational representation and global existence of solu¬ 
tions 

7.1 KPZ as Hamilton-Jacobi-Bellman equation 

Here we show that for every Y € Frbe an d every ho € for (3 > 0 there are global 
in time solutions to the KPZ equation. The idea is to interpret the solution as value 
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function of an optimal control problem, and to “guess” an expansion of the optimal 
control. This can be made rigorous in the case of smooth data and allows us to obtain a 
priori bounds which show that no explosion occurs as we let the smooth data converge 
to Y. 

Let h solve the KPZ equation 

Sfh = (D hf - c v + 9, h(0) = h, (55) 

for 9 G C 100 ), c v G M, and h G C°°. Then by the Cole-Hopf transform we 

have h = logic, where 

J£w = w(9 — c v ), w(0) = e h . (56) 

We specified in Section 4.2 how to interpret this equation. But as long as Y (see Defini¬ 
tion 4.1) is in C(K+, C°°), we have the following simpler characterization. The function 
w: M+ xTgK solves (56) in the sense of Section 4.2 if and only if w = e Y w 1 , where tc 1 
is a classical solution to 

SYw x = (\BY\ 2 -c v )w x + 2BYBw x , ^(0) = e h ~ Y ( 0 \ 


For the rest of this section it will be convenient to reverse time. So fix T > 0 and define 
V = <p(T — t) for all appropriate ip. Then 




G 




and 


to = e^tF 1 


(dt + A) h = -((D h Y - c v ) - 9 , h (T) = h, 


with 


(dt + A + 2D i FD)tF 1 = —(|D^ | 2 —c v )tc 1 , t v\T) = e h -^ {T) . 

The Feynman-Kac formula ([KS88], Theorem 5.7.6) shows that for t G [0, T] we have 


(57) 


\=E ttX e S(7T)-^(T, 7 T)+/ t T (|D^| 2 ( S , 7s )-c v )d S 
where under P^ x the process 7 solves 

'fs = J 2BY (r, 7r)dr + B s , s > t, 

with a Brownian motion B with variance 2 (that is d(B) s = 2ds), started in Bt = x. An 
application of Girsanov’s theorem gives 

tU(t,x) = e ^(AD EtT U(BT)-V(r.B T )+/ t T (|Dt 7 | 2 (s.B s )- C v )ds e / t T Dt 7 (s,B s )dB s -/ i T |Dt 7 | 2 ( S ,B s )d/ 


= E 


t, X 


Ji{B t )-(Y (T,B T )-t (t,x)-J f T D < V(s,B s )dB s )-c v (T-t) 


Now if we formally apply Ito’s formula, we get 


V(T,B T )-^(t,x)- [ B i Y(s,B s )dB s = [ (d s + A)t p ( S ,B s )ds = - [ *6(s,B a )ds. 

Jt Jt Jt 
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Using the temporal continuity of Y (and thus of W(t,x) = f f T 9 (s,x)ds), it is actually 
possible to show that the right hand side makes sense as an L 2 (Pt i3; )-limit. But we will 
not need this and simply take the left hand side as the definition of the right hand side. 
We will however need the following generalization of the Boue-Dupuis [BD98] formula 
which has been recently established by Usttinel: 

Theorem 7.1 ([Ustl4], Theorem 7). Let B:[0,T ] —>• be a Brownian motion with 

variance a 2 and let F be a measurable functional on C([0, T]; M d ) such that F(B) € L 2 
and e~ F ( B ^ € L 2 . Then 


— logE[e i ‘ (i? ^]=infE 


F[B+ I u. s ds ) + ^ 


2 ds 


where the infimum runs over all processes v that are progressively measurable with respect 
to (Ft), the augmented filtration generated by B, and that are such that u ^ v s (uo) is 
F s -measurable for Lebesgue-almost all s € [0,T]. 


Corollary 7.2. Let h solve (55), let T > 0, and let B be a Brownian motion with 
variance (B) t = 2 1, started in Bo = 0. Then 


h(T , x) = log t v (0, x ) = log E 
= — inf E 


= supE 

V 


exp + Bt ) + J ( 9 (s,x + B s ) — c Y )ds 

(^(s,'y V s)~c v )ds + jJ^ |n s | 2 ds 
rT i 

Kl t) + J (^'(' S >7j)-c v --|n s | 2 )ds , (58) 


where we wrote = x + B t + Jjj n s ds and the supremum runs over the same processes 
v as in Theorem 7.1. 


Define now for any 7 of the form d74 = vtdt + d Bt the payoff functional 

^(7 ,v) = Htt) + l (V(s,7 a ) -c v --|n s | 2 )d.s 

:=h(~/ T )-^Y(T,'YT) + i Y(0,x) + j^ bV (s,'y s )d^y s - (c v + ^|n s | 2 )ds, 
so that h(T, x) = sup,, E[d>(7 1 ', n)]. Plugging in d7^ = vtdt + d B t , we get 

pT -1 

d*(7 ,, ,n) = h( 77.) - t 7 (T,7 j.) + ^Y (0,x) - J ^ — uD*) 7 + c v + -\v\ 2 ^j (s,jg)ds + mart. 

= %£) - V (T, 7£) + V(0, x) - £ ( - |D^ 7 | 2 +c v + ±\v- 2D^I 2 ) (s, 7 «)ds 
+ mart., 
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where we write “mart.” for an arbitrary martingale term whose expectation vanishes 
under E. Now change the optimization variable to vj = vt — 2D^Y (t, 7^), so that 


if =x + B t + /V + 2D^)(s, 1 v s )ds, 

Jo 

and the payoff becomes 

$('y v ,v) = hfiyfi) - P(T, 7^) + V (0,x) + (jD^ 7 | 2 -c V - (s, 7 J)ds + mart. 


4—V 

In the following denote X 1 = D Y 1 . We can iterate the process by considering Y , 
where Y v is as in Definition 4.1 (i.e. *Y solves (dt + A)^ 7 = — (|^Y| 

terminal condition ^Y T = 0), which allows us to represent the payoff function 


' 2 -c v ) 


with 


as 


$(7", v) = hffl) - P ( T, 7 j.) + V (0, x) + ^(0, x) 

+ J +2^^ —-lu 1 ) 2 ^ (s, 7 ^)ds + mart. 

= hffl) - V ( T , Yt) + ^ (0, X) + Y(0, x) 

+ J \ 2 +2X^X — — 2*X | 2 ^ (s,7g)ds + mart. 

We change the optimization strategy to v 2 = v\ — 2X (t, 7^), so that the dynamics of 
r y v read d 7 ^ = d B t + (v 2 + 2X + 2X )(t, 7^)dt. Now let 

2Sf Y r = \X s/ \ 2 +2XX v + 2(X + X V )DY R , Y r ( 0) = 0. (59) 

This is a linear paracontrolled equation whose solution is of the form Y R = Y^ + 
Y' -4< P + Y^ and depends continuously on Y. Moreover, another application of Ito’s 
formula yields 


*(7 W , v) = h{l v T ) ~ V (T, 7 V T ) + t 7 (0, x) + VY x) + ^ R (0, x) 

+ J (v 2 ^ — i|u 2 | 2 ^ (s,7^)ds + mart. 

Let us summarize the result of our calculation: 

Theorem 7.3. Let 8 G A?C^ 2 (M, C °°) and c v ,c^ € M, and let Y = 0k pz ($, c v , c^) G 
Tkpz Definition 4-1). Let h solve the KPZ equation (55) driven by Y, started in 
h G L°°, let Y r solve (59), and let T > 0. Then 

(h — Y — Y v — Y r )(T, x) 
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where 


= supE [h(&) - Y( 0, Ct) + £ (( x R \ 2 -\\v ~ 2% R \ 2 ) ( s , C)ds], (60) 

Ct=x + [ (2*X + 2^ V + v)(s, C)d a + B t (61) 

Jo 

and, the supremum is taken over the same v as in Theorem 7.1. 

This representation is very useful for deriving a priori bounds on h. 

Corollary 7.4. In the setting of Theorem 7.3, for all T > 0 there exists a constant 
C > 0 depending only on T and ||Y||y kpz (:r), such that 

INIc t l^C(i + INIl°°). 

In particular, if (Q^ pz (9 n ,a/(,cY))n is a converging sequence in 34pz and every 6 n is as 
above, then the corresponding solutions ( h n ) stay uniformly bounded in CtL°° for all 

T > 0. 

Proof. First, we have —\\v — 2fl | 2 (s,Cs) ^ 0, independently of v, and therefore 
supE[/i(C£) - y(0,$) + £ (l^| 2 > - 2^ i? | 2 )(s,C)ds 

< \Ml-°+\\Y(0)\\ L o O +\\X r \\ 2 L 2 t loo- (62) 

Moreover, choosing the specific control v(t) = 0 we get 

su P E[h(C^)-T(0,Cr) + ^ r (l^ R | 2 -^-2^| 2 )( S ,C)d. S ] ^ HHU~-||y(0)|U=o. 

In combination with Theorem 7.3 and (62), this yields 

IIMr)IU~< ||(y + Y v + Y R )(T)\\ L oo + \\Y(0)\\ L oo + \\h\\ L oo+\\X R \\ 2 L 2 LOO , (63) 

which gives us a bound on the solution of the KPZ equation which is linear in terms of 
the data Y € Tkpz> and quadratic in terms of the solution Y R to a linear paracontrolled 
equation that can in turn be bounded by the data Y. □ 

An immediate consequence is the global in time existence of solutions to the KPZ 
equation: 

Corollary 7.5. In Theorem f.2 and Theorem f.9 we have T* = oo. If in Theorem 3.5 
and Theorem 6.12 the initial condition is uq = Dho for some ho € c £ 2a+l (respectively 
ho € then also here T* = oo. 
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Proof. To see that T* = oo in Theorem 4.2, Theorem 4.9 and in Theorem 3.5 under 
the condition uo = D/io, it suffices to combine Theorem 7.3 with Corollary 4.10. As for 
Theorem 6.12, recall from the proof of this theorem that in order to extend a solution 
u € S>^(T) from [0,T] to [0, T'] with T' > T, it suffices to solve Burgers equation on 
[0,T 7 — T] with the initial condition vP(T) — u'(T ) -<Q(T) € < ^ >2a for the remainder. 
Moreover, if uq = Dho, then also u(T) = D h(T), so that the first part of the theorem 
tells us that u can be extended from [0, T] to [0,oo). □ 

Remark 7.6. The condition uq = D/io is equivalent to J T ito(x)d2; = 0. In case the 
integral is equal to c € M \ {0} we can consider u = u — c which solves 

T£u = Du 2 + 2cD7t + D£, u{ 0) = D/ro 

for some ho. This is a paracontrolled equation which we can solve up to some explosion 
time, and we have u = D h for the solution h to 

P£h = \Dh\ <>2 +2cDh + £, h( 0) = h 0 . 

The Cole-Hopf transform then shows that h = logic, where 

J£w = 2cDw + w o £, w(0) = e h ° , 

and based on these observations we could perform the same analysis as above to show 
that the explosion time of u (and thus of u) is infinite. We would only have to replace 
the Brownian motion B by the process ( B t + 2ct)t^o which corresponds to the generator 
A + 2cV. 

Another simple consequence of (60) is a quantitative comparison result for the KPZ 
equation. 

Lemma 7.7 (“Comparison principle”). In the setting of Theorem 6.13 let hi, hi € 
and let h\ solve 

jSf/zi = |D/ii| o2 +e, hi(0) = h 1 , 

and hi analogously with h\ replaced by hi. Then 

h\(t,x) + inf (hi{x) — hAx)) ^ hi{t,x) ^ h±(t,x ) + sup(/i2(x) — h\{x)) (64) 

X X 

for all (t,x). In particular, \\h\ — hi\\c T L°°^ ||^i — hi\\L°° for all T > 0. 

Proof. Consider regular data (Y n , fi”, hlf) that converges to (Y, h\,hi) in Tkpz x x 

, and denote the corresponding solutions by hf and hi, respectively. For every n 
the representation (60) and the decomposition lif = h\ + (hf — h\) gives 

hf(t,x) + mfthlfix) — hi(x)) ^ ^ hf(t,x) + swp(h%(x) — h'f(x)). 

X X 

Letting n tend to infinity, we get (64). □ 
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Remark 7.8. Corollary 7.5 gives a pathwise proof for the strict positivity of solutions to 
the rough heat equation started in strictly positive initial data. The classical proof for the 
stochastic heat equation is due to Mueller [Mue91], whose arguments are rather involved; 
see also the recent works [MF14,CK14], Compared to these, our proof has the advantage 
that it does not use the structure of the white noise at all, so that it is applicable in a 
wide range of scenarios. The disadvantage is that we need to start in strictly positive 
data, whereas Mueller’s result allows to start in nonpositive, nonzero data. 


7.2 Partial Girsanov transform and random directed polymers 

To formulate the optimization problem (60) for non-smooth elements ¥ of Tkpz > we first 
need to make sense of the diffusion equation (61) for £ v in that case. This can be done 
with the techniques of [DD15] or [CC15], which we will apply with a slight twist. 

Let us start by formally deriving the dynamics of the coordinate process under the 
random directed polymer measure. This is the measure given by 

d<Q>T ,2 = exp 

where £ is a space-time white noise (and thus f as well), and under P x the process B 
is a Brownian motion started in x and with variance 2. The term —oo T is chosen so 
that has total mass 1. If now h solves the KPZ equation with /i(0) = 0 and if 
h (s) = h(T — s), then we can write 

0 = ^h(0,x)+ [ (d t + A)*h(t, B t )dt + [ D*h(t,B t )dB t 
Jo Jo 

= *h(0,x) + f (— *£(t, B t ) + oo)dt + f D*h(t, B t )dB t - f \D*h (t, B t )\ 2 dt, 

Jo Jo Jo 

which shows that under Qt,x the process B solves the SDE 



B t = x + 


<r-~ 


2D h ( t, B t )dt + d Wt, 




where IT is a Q 7 ^-Brownian motion with variance 2 . 

Of course, a priori such a diffusion equation does not make any sense because D h (t, • 

is a distribution and not a function. But let , and be as in Section 7.1. Then 

we can rewrite the term in the exponential as 




oo)d t = Y ( 0 , x) + ¥ v (0, x) + T R (0, x) - Y (T, B T ) 

+ [ T (X + < X V )(t,B t )dB t - l T \Cx + < x')(t,,B t )\ 2 dt 
Jo Jo 

tT , _ 

+ / X { S ,B s ){dB s -2{X + X w ){s,B s )ds). 

Jo 
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Let us define 


Z = (X +X s 


and set 


dPf x = exp 


( j T *Z (t, B t )dB t - J T \*Z(t, R t )| 2 d#)dP. 


Then under P£ the coordinate process B solves 


B t = x+ / 2Z(s,B s )ds + B t , 


(65) 


where B is a (variance 2) Brownian motion under P t, x - The advantage of this splitting of 
the Radon-Nikodym density is that now we singled out the singular part of the measure, 
and Qt,x is absolutely continuous with respect to P^ x , with 

eM-Y(T,B T )+S^%,B,)dB,) . 

d^lT, X — - , R - dr Tx . 

E P z x [eM~Y(T, Bt) + f 0 T X (a, R s )dR s )] 

This density is strictly positive and in T p (Qr,a;) for all p G [l,oo), even for general 

Y G Tkpz (not necessarily smooth). It remains to construct the measure P^ for general 

Y G y kpz . In general P j, will be singular with respect to the Wiener measure. 

Theorem 7.9. Let Y G Tkpz and T > 0. Consider for given <pt G < ^ Q+1 and f G 
C([0, T], L°°) the solution (p to the paracontrolled equation 

(d t + A)(^ = -2^D^ + /, (p(T) = <p T . 


Then for every iGl f/iere exists a unique probability measure P on 12 = C([0, T], T), 
snc/i that Py 3.(70 = x) = 1 and /or all ip as above the process 

Tt) ~ [ f{s,'y s )ds, t G [0, T], 

Jo 

is a square integrable martingale. Here 7 denotes the coordinate process on Li. Moreover, 
assume that Y n = @k P zi^nXnxY) f or a sequence 9 n G «SfC' 1 “{ 2 (M, C 00 ) and c^,c^ G R 
is such that (Y n ) converges to Y in 34pz- Let for every n the process 7™ solve the SDE 

It = x + [ 2 < Z n (s, 7 ")ds + B t 
Jo 

for Z n = X n + X% = D (Y n + Yff), which is well posed because Z n is Lipschitz continuous 
by assumption, uniformly in t G [0, T]. Here B denotes a variance 2 Brownian motion. 
Then there exists a Brownian motion B on Li with respect to the measure P^ x , such that 
(7 n ,B) converges weakly to (7 ,B). Finally, 7 is a time-inhomogeneous strong Markov 
process under the family of measures (Pf^sT- 
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Proof. Everything follows from Theorem 2.6 and Section 3 in [CC15], provided we can 
show that Z is a “ground drift” as defined in that paper. In other words we need to show 
that Q z o Z is given, for the solution Q z to £ZQ z = DZ, Q z ( 0) = 0. But Q z = Q + Q v , 
where 22?(3 V = DX V . Since X v € do 20-1 and thus Q v € C"^ 2a , the only problematic 
term in the definition of Q z o Z is QoX. And since Qo X is contained in ¥, the proof 
is complete. □ 

Remark 7.10. We expect that for 0 < s < t < T under (F z x ) xG i the transition function 

T s ,tg{x ) = F v z[g{nt)\ls = x\ 

has a density p s j(x,y) which is jointly continuous in ( s,t,x,y). For that purpose we 
would have to solve the generator equation 


(d r + A + ZD)ip = /, r € [0,t), ip(t) = <p t (66) 

for general t € [0, T]. Taking f = 0, we would have T s j(p t (x) = <p(s,x). Then we could 
enlarge the class of terminal conditions pf, and using the techniques of Section 6.3 it 
should be possible to take pt as the Dirac delta in an arbitrary point y G T. But then 
the transition density p s ,t(x,y) ofT s j must be given by p(s,x), where p solves (66) with 
terminal condition 5(y). 

We also expect that the density is strictly positive: In the upcoming work [CFG15] 
it will be shown that if u solves a linear paracontrolled equation and u(0) € C'(T) is a 
continuous, nonnegative and non-zero function, then u(t,x ) > 0 for all t > 0, x € T. 
Combining this with the smoothing effect of the linear equation (66), which takes as 
the terminal condition at time t the Dirac delta and returns a nonnegative non-zero 
continuous function for all sufficiently large r < t, the strict positivity of the density 
should follow. 

Remark 7.11. In the setting of Theorem 7.9, we can construct the “full random directed 
polymer measure ” Qt,x by setting 

dQ T,x _ exp (~Y(T, 7 T ) + J Q r ^X (s, 7 s )dR s ) 

dF Tx Epz a [exp(—E(T, 7t ) + f 0 T 5 R (s, 7s )d B s )j ’ 

which is now a perfectly well defined expression. 

Similarly we can also construct the measure Qx,(t,x) under which the coordinate 7 
process formally solves 

'y s = x + Jt 2T)*h (r, 7r )dr + (B s — B t ), s G [t,T], 

With respect to Q r,(t,x) we have the following Feynman-Kac representation for the solu¬ 
tion w to the paracontrolled equation heat equation driven by ¥ and started in w(0) = w: 


w(t,x) = EQ T(tix) [u;( 7T )]r<;(t,x), 

where w(t,x) solves the same equation as w but started in w(0) = 1. 
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7.3 Rigorous control problem 

The control problem that we formulated in Section 7.1 worked only for “smooth” ele¬ 
ments of Tkpz; that is for ¥ = 0k pz (0, c v , c®') with 6 G C°°) and c v , c^ G R. 

Let us use the construction of Section 7.2 to make it rigorous also for general Y € 34pz- 
For that purpose fix T > 0 and set Qt = C([0, T],T), equipped with the canonical 
filtration {Ft)t&[o,T]- Recall that a function u: [0, T] x Ut —>• R is called progressively 
measurable if for all t G [0, T] the map [0. t] x Qt 9 (s,w) eG u(s,w) is &[0,t\ <8> JFt~ 
measurable. 

From now on we also fix ¥ G Tkpz and define Z as in (65). 

Definition 7.12. Let v be a progressively measurable process and let (fi, J 7 , (J r t)te[o,T\ > I 6 ) 
be a filtered probability space. A stochastic process on this space is called a martingale 
solution to the equation 


'y t = x+[ ( 2 *!z(s, 7 s ) + v(s,'y))ds + B t (67) 

Jo 

if P( 7 o = x) = 1 and whenever (pT G c (o a+l and / G C([0,T],L°°) and 99 solves the 
paracontrolled equation 

(d t + A + 2 i ZD)ip = f, <p(T) = ip T . 


on [0, T], then 

‘P&'Yt)- [ (f(s, 7s) + D(p>(s, 7 s )u(s, 7 ))ds, fG[0,T], 

Jo 

is a martingale. 

We write pm for the set of progressively measurable processes, and for a given v G 
pin and x G T we write 9Jt(v,x) for the collection of all martingale solutions to (67). 
Note that here we explicitly allow the probability space to vary for different martingale 
solutions, and also that we do not make any claim about the existence or uniqueness (in 
law) of martingale solutions. 

Theorem 7.13. Let h — Y (0) G c £ 2a+l and let h be a paracontrolled solution to the KPZ 
equation with initial condition h. Then 

(h — Y — T v — Y r )(T, x) (68) 

= sup sup E /i( 7 t) — Y(0, 7 t) + [ \ 2 -^-\v — 2X | 2 Vs, 7 )ds 

uepm jem(v,x) L Jo V 4 / . 

and the optimal v is 

v(t, 7 ) = 2D (t -*X- ^ V )(t, 7t ) = 2(A R + D t R )(t, 7 t), 
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where h R solves the paracontrolled equation 

&h R = \X R \ 2 +2(X + X W + X R )Dh R + \Dh R \ 2 , h R (0) = h-Y(0). (69) 

For this v and all ieT, the set 9Jt(v,x) is non-empty and every 7 G %Jl(v,x) has the 
same law. 

Proof. Clearly h solves the KPZ equation driven by ¥ and started in h if and only if 
h R = h — X — X v — X R solves (69). Moreover, the paracontrolled structure of h R is 

h R = 2(Dh R ) ^P + h RA 

with h R ’^ G J5f 2 “ +1 , and in particular h G CTf a+1 . Reversing time, we get 

{dt + A + 2ZD) h =-\X | 2 — 2a D h — |D h | 2 , h (T) = h-Y{ 0). 


Let now v E pm and let 7 be a martingale solution to 

rt 


7 t = x+ [ ( 2 ^(s, 7 s) + u(s, 7 ))ds + B t . 
Jo 


-R 


Since h - ¥(0) € tf 2a+1 , we can take h as test function in the martingale problem 
and get 

— 4 — R i-f-R.o 4 — R 4 —R 4 — 

/i( 7 t) — y (0, 7 t) = h (0, x) + (-\X | 2 —2a D h - |D h \ 2 )(s, ls )ds 

Jo 

f T 4-R 

+ / D h (s, 7 s )v(s, 7 )ds + mart., 

Jo 

4—R 4—R 

writing v(s, 7 ) = v(s, 7 ) — 2 A (s, 7 ^) — 2D h (s, 7 S ) we obtain 


so 


E 


Kjr) - ¥( 0 , 7 t) + (pt R \ 2 -^\v - 2 ^ R | 2 )(s, 7 )ds -^ R (0,.x) 


= E 
+ E 
= E 

This shows that 


f / G - R ^— R ^ —-^,0 ^— R / G~ R ^^ 

/ (-|A | 2 —2a D h — |D h | 2 +D h (2a + 2D h +u))(s, 7 )d4 

Jo 


0 


^|^ R | 2 -i|u + 2D^h | 2 ) (s, 7s)ds 


'0 




sup sup E 


Kit) - ¥(0,7t) + J (pl R \ 2 -^\v - 2^' R | 2 ^)(s,7)di 


^ /i fi (T, x). 
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4——R ■i—R 

On the other side, taking v = 2X + 2D h we obtain a ground drift in the terminology 
of [CC15], and therefore for all x G T there exists a 7 € 9R(v,x) and its law is unique. 
For any such 7 we obtain 


E 


h(jr)-Y(0,^ t ) + £ (\ < X B 


\ 2 --\v — 2*X 
4 


R. 


(s, 7 )ds = h R (T,x), 


and thus the proof is complete. 


□ 


Remark 7.14. We only needed h — Y (0) G < ?f 2cH ~ 1 to apply the results of [CC15]. By 
using similar arguments as in Section 6 it will be possible to weaken the assumptions on 
the ground drift (allowing a possible singularity at 0) and on the terminal condition in the 
martingale problem in [CC15], and then the variational representation of Theorem 7.13 
will extend to Ji G for /3 < 1 as in Theorem 6.13. 

Actually, the extension to h G c £ OL is immediate because we can simply start Y R in 
— ^—R ^—R d 

h — Y(0) (which still gives us a ground drift Z + 2X +2D h , where h H is now started 

in 0), and change the control problem to 


sup sup E 

'uGpm 'y£$Jl(v,x) 


(|S H p-l|„-2S^)( 5 , 7 )d S . 

Similarly, if we start Y R in —T(0) the control problem has the more appealing form 


sup sup E 

vEpm 


Hit) + 




(s, 7 )ds 


8 Convergence of Sasamoto-Spohn lattice models 


In this section we consider the weak universality conjecture in the context of weakly 
asymmetric interface models (/9/v : ®i+ x Z^r —> M (where Z^r = Z/(IVZ)) with 

dUp N (■ t , x) = A Zn ip N (t, x)dt + Vs(Bz n (D Zjv p N (■ t ), T>z n Pn (t)))(x)dt + d W N (t,x), 
<p N {0,x) = <Pq(x), (70) 


where A% N and T)% N are discrete versions of Laplacian and spatial derivative respectively, 
B% n is a bilinear form taking the role of the pointwise product, (IF/v(t, x))t£gL + , x &7, N is an 
IV-dimensional standard Brownian motion, and pff is independent of VFy. We assume 
throughout this section that 

2 t r 


e = 


N 


Equation (70) is a generalization of the Sasamoto-Spohn discretization of the KPZ equa¬ 
tion, see Remark 8.3 below. To simplify things (eliminating the need to introduce renor¬ 
malization constants), let us look at the flux D z n Pn- Assume that there exists j3 < 1 
for which (x i-g ( Dz n Pq){x/£))n converges weakly to 0 in with rate of convergence 
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e 1 / 2 . Then {(t,x) i->- T>z N <PN(t/£ 2 ,x/e)) converges to 0 (this will be a consequence of 
our analysis below), and we can study the fluctuations defined by 

u N (t,x ) = £ _ 1 / 2 D Zjv <^a i(t/e 2 ,x/e). 

This is a stochastic process on M + x TV with TV = (eZ) / (27rZ) which solves the SDE 

d UN(t,x) = x)dt + (DNBN(uN(t),UN(t)))(x)dt + d(DAr£ - 1 / 2 lTV(i, x)) 

u N (0) = u$. (71) 

where A tv, D n, and Bn are approximations of Laplacian, spatial derivative, and point- 
wise product respectively, (9 jDtv£ -1 / 2 1TV converges to D£, where £ is a space-time white 
noise and Uq (. x ) = Djve 1 / 2 ^^( x l e )■ We show that if Djve 1 / 2 ^( x / £ ) converges in distri¬ 
bution in then («at) converges in distribution to the solution of a modified Burgers 
equation involving a sort of Ito-Stratonovich corrector. 

Another way of reading our result is that (71) is a lattice discretization of the Burg¬ 
ers equation and we show that it might converge to a different equation in the limit, 
depending on how we choose Aat, D^r, and Bn- 

There are two problems that we have to deal with before we can study its convergence. 
First, it is not obvious whether un blows up in finite time, because the equation contains 
a quadratic nonlinearity. Let therefore £ be a cemetery state and define the space 

Cn = W'- ®V ®l Tjv U {C}> V 9 is continuous on [ 0 , and + t) = (,t ^ 0 }, 

where 

r c (<V = inf{t ^ 0 : ip(t) = C} and for c > 0 r c (ip) = inf{f ^ 0 : |V(t)||L°°(T JV )^ c}, 

with ||Cllz, oo (Tjv) = °°- Then a stochastic process un with values in Cn is a solution 
to (71) if t c (u n ) = sup c>0 r c (u N ) and ujv|[o,r c (iw)) solves ( 71 ) on [0 ,t^(u n )). It is a 
classical result that there exists a unique solution in that sense (which is always adapted 
to the filtration generated by and I!V, but we will not need this). 

The next problem that we face is that it/v(i) is only defined on the grid TV and not on 
the entire torus T. Since we will only obtain convergence in a space of distributions and 
not of continuous functions, some care has to be exercised when choosing an extension 
of un to T. For 5 > 0 one can easily define sequences of smooth functions (/V) and 
(<7 jv) on T such that /V and qn agree in the lattice points IV, but both converge to 
different limits in Here we will work with a particularly convenient extension of 

(un) that can be constructed using discrete Fourier transforms [HM12, HMW14], Since 
it will simplify the notation, we make the following assumption from now on: 

N is odd. 

Of course, our results do not depend on that assumption and we only make it for con¬ 
venience. We define for ip: TV —>• C the discrete Fourier transform 

& N <p(k) = £ J2 v( £ t) e ~ ieEk i k € Zj\r, 

\£\<N/2 
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(for even N we would have to adapt the domain of summation), and then 
£n<p(x) = (27t) _1 ^2 ^Nf{k)e lkx , x G T. 

\k\<N/2 

Then £np{x) = fix) for all and by construction £j\up is a smooth function with 

Fourier transform ^£]\np(k) = ^Mf{k) lifci<jv/2- If T is rea l valued, then so is £n ( P- 

Remark 8.1. A quite generic method of extending a function f.Tw —>• M to T is as 
follows: write 

= ^2 T(x)ed x , 
a;eTjv 

let tin G @'(T), and define (p(x) = if * rj^, which yields for all k G Z 
JFf(k) = JPfj(k)^r]N(k) = ^N(p(k)^rjN{k). 


Now assume that (<pn) is a sequence of functions on Tjy such that (£nTn) converges 
in S>'(T) to some limit ip oq. Convergence in @'(T) is equivalent to the convergence 
of all Fourier modes, together with a uniform polynomial bound on their growth, and 
thus we get linpv ^NfN[k) = <^foo(k) for all k G Z. So if ( tjn ) converges to 6 q in 
S>'(T), then (<pn) converges in &'(T) to ipoo. Typical examples for interpolations that 
can be constructed in this way are the “Dirac delta extension” (take tjn = <$o f or all 
N), the “piecewise constant extension” (take tjn = e _ 1 l[o,e)A or the “piecewise linear 
extension” (take tjn{x) = e^ 1 ((e _1 .x + l)l[_ e 0 ](* T ) + (1 — £- 1 -' c )I-(o,e ]( x )))■ dn particular, 
our convergence result below also implies the convergence of all these interpolations. 

For ip G Cat we then get £at<A : ^+ —>• fF'(T)U{C} and for general —»• ^'(T)U{£} 

we define 

T cW = inf {f ^ 0 : Uit)\\/3> c} 

whenever f3 G M, with ||C||/ 3 = oo. Then write dc^if,^') = \\if — if'\\c R as 

well as 

OO 

dp{'lp,f) = 2_m ( 1 A + \ T mW - T mW)\))- 

m= 1 

We will need some assumptions on the operators Aat, Dat, Bn- Let 


Aat^(.x) = £ 


f(x + ey)ir(dy), 


T) N f(x) = s 


ip(x + £y)u(dy), 


(72) 


BN(f,i’)(x) 


<p{x + £y)ip(x + £z)/i(dy, dz), 


where ir, u, and p are finite signed measures on Z and a probability measure on Z 2 , 
respectively. We define 


= 9(x) = 4 £AM, h{xi , I2)= [ e ‘i*,^» (izuiZ2) , 


-X* iX Jz 2 

and make the following assumptions on the measures 7r, p: 
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(H /) The finite signed measure 7 r on Z is symmetric, has total mass zero, finite fourth 
moment, and satisfies f z y 2 ir(dy ) = 2. Moreover, there exists c/ > 0 such that 
/(x) > Cf for all x £ [—7r, 7r]. 

(H fl ) The finite signed measure z/ on Z has total mass zero, finite second moment, and 
satisfies J z yi , (dy) = 1. 

(H h ) The probability measure y has a finite first moment on Z 2 and satisfies y(AxB) = 
y(B x A) for all A, B C Z. 

The constant c/ in (Hy) exists for example if 7r = Ylk>iPk{$k + S-k) — «>o where pk ^ 0, 
1 2 Pk = c and pi > 0. 

Theorem 8.2. Make assumption (Hf), (H g ), (H^) and let /3 £ (—1,—1/2). Consider 
for all N £ N an N-dimensional standard Brownian motion (W]y(t, x))t^>o,x£T N and an 
independent random variable (uq (x)) xG j n and denote the solution to (71) by un. If 
there exists a random variable uq such that Sn'Uq converges to uq in distribution in ^, 
then (Snun) converges in distribution with respect to the metric dp to u, the unique 
paracontrolled solution of 

A£u = Du 2 + 4cDu + D£, u(0) = uq, (73) 

where f is a space-time white noise which is independent ofuo, and where 

c = r Im(g(x)h(x)) h(x, -x)| g (x)| 2 d ^ ^ R 

1- X |/(x)| 2 

Proof. We have 

& N B N (pMk) = {2 k)~ 1 V & N <p(i)& N 1,{k-i) [ e^ +£ M*V(dy,dz), 

from where we deduce that £/v-B/v(<£, VO = rbv-B/v(£/v75 £/vV 0 for all £ M Tjv , with 

Un^p(x) = (27T)- 1 £ e ikNx ^p(k), 

k 

where k N = argmin{|7|: £ £ Z, £ = k + jN for some j £ Z} £ (—JV/2, N/2). Moreover, 
we obtain for \k\< N/2 

E [e~ 1/2 &£ N W N (t, ki)e- 1/2 &£ N W N {t , k 2 )\ 

= e Z[W N (t,££i)W N (t,££ 2 )}e~ i£{hkl+e2k2) 

\tl\,\t2\<N/2 

= £ 'Yh te~ ll{kl+k2)e = t.2ir5k 1+ k 2 =0, 

\e\<N/2 
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which shows that c>t£ _ 1 / 2 £ 7 vhFjv(t) has the same distribution as Vn£, where £ is a space- 
time white noise and VnT = l(-jv/2,iV/2)(D) = ^~ 1 {^(-n/2,n/2)^t) is a Fourier cutoff 
operator. 

We now place ourselves on a probability space where such a space-time white noise 
is given and where (uq) is a sequence of random variables with values in 0'(T), which 
is independent of £, such that Uq has the same distribution as SnUq for all N , and 
such that Uq converges to a random variable u$ in probability in (which is then also 
independent of £). Then un has the same distribution as ujy, the solution to 

(d t - A N )u N (t, x ) = D n U n B n (u n , UN)(t, x) + B N VN^(t, x), un(0) = Uq , 

and therefore it suffices to study un- The pathwise analysis of this equation is carried 
out in Section 8.2 below, and the convergence result assuming convergence of the data 
(Xjv(£)) and the random operators (An) is formulated in Theorem 8.17. The conver¬ 
gence of (Xjy(£)) and (An) in L p spaces is shown in Theorem 10.1 and Theorem 10.4, 
respectively. Thus, we get that (un) converges in probability with respect to dp to u, 
the solution of (73). □ 

Remark 8.3. If we take Bn as the pointwise product for all N, and An as the discrete 
Laplacian, Aat/(,t) = £~ 2 (f(x+e)+f(x—e) —2f(x)) andY)Nf(x) = e~ l (f(x)— f(x— e)), 
then we get c = 1/8, so the additional term in equation (73) is l/2Dw. 

However, if we take the same A n and D^v but replace the pointwise product by 

B N (<p,l/))(x) = 1 (KLp(x)i>(x) + \(<p(x)i>(x + £)+(p(x + £)i>(x)) + K(p(x + £)i>(x + e)) 

2{k + A) 

for some k,A 6 [0, oo) with k +A > 0, then one can check that c = 0. Here the Sasamoto- 
Spohn discretization [KS91, LS98, SS09] corresponds to k = 1, A = 1 / 2 . In that case one 
furthermore has 


{(p,D N B N ((p,g>)) TN = ^ <p(x)D N B N ((p,ip)(x) =0, 

which entails that already for fixed N there is no blow up in the system, i.e. un is well 
defined for all times. Moreover, now we can explicitly write down a family of stationary 
measures for un-’ For all m £ M, the evolution of un is invariant under 


N-l 


T-r exp (—exj + mxj) 

/ 4 (d*)=n- j 


3=0 


Z £ 


dx 


31 


where Z ^ is a constant normalizing the mass of fi £ m to 1; see [SS09] or simply verify 
that the L 2 iff)-adjoint of the generator of un applied to the density of n £ m equals 0 and 
then use Echeverna’s criterion to obtain the invariance of p £ m from its infinitesimal 
invariance [Ech82]. If un( 0) ~ lF m , then for all t ^ 0 the vector (uN(t, x)) X £j N consists 
of independent Gaussian random variables with variance e _1 /2 and mean m. Therefore, 
for all t > 0 the $)'- va l U ed random variable (SNUN(t, •)) converges in distribution to a 
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space white noise with mean m and variance 1/2. It is also straightforward to verify 
that sup^vE[||£jvWjv(0, Ollso ] < 00 whenever a < —1/2, and then the Besov embedding 

theorem shows that the convergence actually takes place in distribution in , for as 
required in Theorem 8.2. But if £nun is a stationary process for all N, then any limit 
in distribution must be stationary as well, and this shows that the white noise with mean 
m and variance 1/2 is an invariant distribution for the stochastic Burgers equation. 
This is of course well known, see for example [BG97] or [FQ14]. But to the best of our 
knowledge, ours is the first proof which does not rely on the Cole-Hopf transform. 

We now take 


NUN(t , x) = (d t - A n )un (t, x ) = D N n N B N (ujv,ujv)(t, x) + D N VNf(t, x) (74) 

as the starting point for our analysis, where we assume that «at(0) = Vnun( 0)- Recall 
that 

n N ip(x) = (27 r)- 1 J2 eikNx ^T(k) 


with 


k N = argmin{|£|: i € 7L,I = k + jN for some j 6 Z} € (—N/2,N/2), 


and that 

VnT = 1(-AT/2,A72)(D) = & 1 {^ L (-N/2,N/2)^ l P)- 

The operators A^v, Djv, and B jv are given in terms of finite signed measures n, u, p, as 
described in (72). 

Lemma 8.4. Let 7 t be a finite signed measure on R that satisfies (Hf). Then the function 


fix) = - 


e ixp 7r(dy) 


x 


1 - cos (xy) 2 

— 7 — 72 —y n i d y) 

(xy) 2 


is in C/ and such that /(0) = 1. 

Proof. The function ip(x) = (1 —cos(x))/x 2 = 2 sin 2 (x/2)/a; 2 is nonnegative, bounded by 
1/2, and satisfies </?(0) = 1/2. Therefore, / is bounded and /(0) = 1/2 J K y 2 7r(dy) = 1. 
Furthermore, it is easy to check that € C 2 , and thus 

|/'(x)|< f \tp'(xy)\\y\ 3 \ir\(dy) < f \y\ 3 \n\ (dy), |/"(*)|< [ |y| 4 M(dy). 

J ]R J M */ ]R 

As 7T has a hnite fourth moment, this shows that / £ C?. □ 

Lemma 8.5. Let v be a finite signed measure onM that satisfies (H g ). Then the function 


g(x) = 


f R e ixy n(dy) 


ix 


r e lxy - 1 

L ixy 


yv(dy) 


is in C/ and such that g( 0) = 1. 
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Proof. It suffices to observe that the function ip(x) = ( e lx — 1 )/(ix) is in C£ and satisfies 
93 ( 0 ) = 1 , and then to copy the proof of Lemma 8.4. □ 

The next lemma is a simple and well known statement about characteristic functions 
of probability measures. 

Lemma 8.6. Let n be a probability measure on M 2 that satisfies (H^). Then the function 

h(x, y) = [ e i ( XZl+yZ2 '>ii(dz 1 ,dz 2 ) 

Jr 2 

is in Cl and such that h{0 ,0) = 1. 

Our general strategy is to find a paracontrolled structure for (74) and then to follow 
the same steps as in the continuous setting. To do so, we need to translate all steps of 
our continuous analysis to the discrete setting. 

8.1 Preliminary estimates 

Fourier cutoff. The cutoff operator Vn is not a bounded operator on < rf a spaces (at 
least not uniformly in N) and will lead to a small loss of regularity. 

Lemma 8.7. We have 

\[Pnv\\l°°< logiVU^Hioc, 
and in particular we get for all 5 > 0 

\\Vnv- v\\ a ~5< iv- 5 iogiv|MU< |M| a . 


Proof. We have 


\\'PnV\\l°°(T:)'?S \\^ 1 1(—7V/2,A r /2) ll^i( T ) IMIl°°(t)> 
and using that N is odd we get 

2) (*) = (2»)-‘ y = - D/ 2 ) -. 

’ , z —■', 1 — cos(x) 




(~N/2,N/2) 


Now 


and therefore 


|fc|<iV/2 

cos(x(N — l)/2) — cos(x(IV + l)/2) 


1 — cos(x) 


< mi 


mm 


— cos(x) 
Nx 2 \x\ 


x 2 ’ x 2 [ ’ 


/ 7T 

\&~ 1 t ( -N/ 2 ,N/ 2 ){x) \ dx < log N. 

-7T 


To obtain the bound for VnT ~ T it suffices to note that Vn acts trivially on Aj (either 
as identity or as zero) unless 2 J ~ N. □ 
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Lemma 8.8. Let a ^ 0 and p G < & a . Then for any <5^0 

- ^|| q _ 5 < N~ s logN\\p\\ Q . 

If supp(J^V) C [— cN, cN] for some c G (0,1), then this inequality extends to general 
a€l. 

Proof. We already know that \\VnP~ N~ s log lV||<^|| a . So since UnT = Tn ((1 + 
e lN ' + e~ tN ')p) we get for a ^ 0 

II Aq,7-V((e lN + e iN )(^)||l°o ^ l 2 ?<ivll^ ? Jv(( e ,JV + e lAr )<^)||l°° 

< 2 l 2 « 3 <Af log(iV)|| (e~ iN ' + e iN ')Aj<p\\L°° 

j:2i~N 

< ]T t 2q < N log(N)2-i a \\cp\\ a 

j:2i~N 

<l 2 ,< Jv log(lV)lV- a ||^|| Q <2-^- 5 ) log(lV)Ar- <5 ||^|U. 

If also support) C [— cN, cN], then the spectrum of 2cos(Nx)u is contained in an 
annulus N s#, and therefore we can replace the indicator function 1 2 q<n by 1 2 <2 ~at in 
the calculation above, from where the claim follows. □ 

Remark 8.9. There exists c G (0,1), independent of N, such that if supp if) C 
[— N/2, N/2], then supp (^(p -< ip)) C [— cN, cN]. This means that we can always bound 
lUv (pP if) — p -< if, even if the paraproduct has negative regularity. On the other side 
the best statement we can make about the resonant product is that if p and if are both 
spectrally supported in [—N/2, N/2], then supp (J^(po if)) c [—IV, IV]. A simple conse¬ 
quence is that ifa+fi > 0, p G c € a , if G ^, and supp(^V)Usupp(^"VO C [—N/2, N/2], 
then 

II Hat MO - Ti/\\aAi3-5< log(AT)II^H q,H-011^. 


Estimates for the discrete Laplacian. 

Lemma 8.10. Let it satisfy (Hf). Let a < 1, (3 G R, and let p G c € OL and if G . 
Then for all 5 G [0,1] and N G N 

\\A N if - Aif\\p- 2 - 6 < N~ 5 \\if\\p and 

l|Ajv(<£ <if)-p< Ajv^|| q+(9 _ 2 < IMUMI/3- 

Proof. Using that n has zero mass, zero first moment, and finite second moment, we 
have 


AjA N if(x)\ 


e 2 / Ajif(x + ey)ir(dy) 

Jr 

£~ 2 / (A jif(x + ey) - Ajif(x) - DAjif(x)£y)ir(dy) 

Jr 
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< ||D 2 A i ^|| L oc / J/ 2 1vr|(d?/) < 2 j(2 /3) 


On the other side, 


\Aj(A N i/) - Afi)(x)\ = 


e / (Ajip(x + ey) - AAjV>(a;))7r(dy) 


< £||D 3 A^||loo / |y| 3 |7r|(dy) < 2"^- 3 ) e 


and thus by interpolation || Ajv'V’ — A^>|| ( g_ 2 -< 5 < e'' [ 

For the commutator between discrete Laplacian and paraproduct it suffices to control 
Ajv(<S'j_i</?AjV ; ) ~ Sj-npAjAfpip in L°°, which is given by 

\(A N (S j -npA j 'ip) - Sj-iLpAjAN^)^ 

£ ~ 2 / (Sj-i<p(x + ey) - Sj-i<p(x))Ajip(x + ey)ir(dy) 

Jr 

< ||D 2 5 j _i^||Loo||A j V’||L°° f y 2 \ir\(dy) 


+ £ 


< 


T)S j -np{x)y(A j ^{x + ey) - A j 'ip(x))Tr(dy) 

R 

2^- a -^M a \mJ y 2 \M(dy). 

Jr 


□ 


While the semigroup generated by the discrete Laplacian Aat does not have good 
regularizing properties, we will only apply it to functions with spectral support contained 
in [— N/2, N/2], where it has the same regularizing effect as the heat flow. It is here that 
we need the assumption that f(x) Cf > 0 for x € [— 7 r, 7 t] . 

Lemma 8.11. Assume that ir satisfies (Hf). Let a € I, fi 0, and let p G & with 
supp (^V) C [—N/2, N/2], Then we have for all T > 0 uniformly in t € (0, T] 

\\e tA »p\\ a+ p<t-V 2 M Q . (75) 


If a < 0, then we also have 





Proof. Let x be a compactly supported smooth function with x = 1 on [— n, tx\ and such 
that fix) Cf /2 for all x G supp(x)- Then e tAN p = e _ ^( £D ) D 2 x(eD)</? =: if E (D)ip. 
According to Lemma 2.2 in [BCD11], it suffices to show that 

max sup|D fc ^ £ (x)|t / 3 // 2 |x| /5+fc ^ C < oo, 

fc= 0 ,l, 2 xG ]R 
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uniformly in e £ (0,1]. For if £ itself we have 

\^ e {x)\tV 2 \x\P< e-^^ 2 \Vixf<l. 

To calculate the derivatives, note that 

(e ~^p{x))' = e^ x \-p'{x)p{x) + p\x)\, 

(e-^p{x))" = e-^ x) [<p\xfp{x) - 2p\x)p\x) - p" {x) p{x) + p" {x)\. 

In our case we set p £ {x) = —tf{ex)x 2 and p £ (x) = x{ £x ), and obtain 
p' £ {x) = —tx(f' (ex)sx + 2 f(ex)), P £ (x) = — t(f" (ex)(£x) 2 + Af'{ex)£x + 2 f(ex)). 

Since \exx( £ x)\< 1 and similarly for (ex) 2 x( £ x) and £xx'( £ x), we get 

\x\\p' £ (x)p £ (x)\<\Vtx\ 2 t x{ex) ^ and \x\\p' £ (x)\<t x{£x) ^ 0 , 
from where we deduce that 

|D x Vfc(®)|i /J/2 |®| /,+1 < e -¥(^) 2 (\Vtx\ 2 ^+\Vtxf) < 1. 


Similar arguments show that also |D^ £ (.'c)|t^ //2 |x|^ +2 < 1, which concludes the proof 
of (75). The L°° estimate now follows from an interpolation argument. □ 


Corollary 8.12. Let ir satisfy (Hf), let a £ (0,2), and let p £ c £ a with spectral support 
in [-N/2, N/2]. Then 

||(e* AiV — id)<£>||£,oo< t"/ 2 ||<yj|| a . 

Proof. By definition of e tAjv we have 


- idyllise ^ 


-( 2 -«)/ 2 ||A A r(p|L_ 2 ds < f "/ 2 


where we used Lemma 8.11 and Lemma 8.10 in the second step. 


□ 


Combining Lemma 8.11 and Corollary 8.12, we can apply the same arguments as in 
the continuous setting to derive analogous Schauder estimates for (e tAjV ) as in Lemma 2.9 
or Lemma 6.6 - of course always restricted to elements of SP' that are spectrally sup¬ 
ported in [—1V/2, IV/2]. 


Estimate for the discrete derivative. 

Lemma 8.13. Let v satisfy (H g ). Let a € (0, 1), f3 £ R, and let p £ and if £ . 

Then for all 5 £ [0,1] and N £ N 

IIDjvV 7 - I>0ll/3-i-5^ N ~ S W\p and 

\\B N (p-<ip) - pPT> N ip\\ a+ p-i< IMUHI/3- 
The proof is the same as the one of Lemma 8.10, and we omit it. 
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Estimates for the bilinear form. Let us define paraproduct and resonant term with 
respect to B n‘. 

B N (ip~<ip) = ^B N (Sj-np,Ajip), B N (pyip) = ^ B N (Aj(p, Sj-ii/;), 

3 3 


BnU° 9)= B N (A i f,A j g). 


We have the same estimates as for the usual product: 

Lemma 8.14. Let p satisfy (Hh). For any (3 € M and 6 G [0,1] we have 
\\B N ((p^lp) -ip^ll>\\p-8< N- S \\ip\\ L oo\\lp\\p, 
and for a < 0 furthermore 

\\B N (tp Ip) ~ ipPlp\\ a+ p-5< iV"''||^|| a ||V’|| / 3. 

For a + (3 — 5 > 0 we have 

\\B N (<poip) - ipo%j)\\ a+ p_ s < iV- <5 ||(^|| a ||^|| /3 . 

Proof. It suffices to note that ^B^(Ai(p, A ftp) and &(AnpAj'ip ) have the same support 
and that \\B^(Aitp, A^-i/Olli- 00 ^ ||Aj<^||L°o ||Aj^>||x,°o, whereas 

\\B N (Ai<p, Ajip) - Ai<pAjip\\ LO o< N~ l (2 l + 2- 7 )||A i ^|| Z/ °o ||Aj*/>|| l <x>. 


□ 


To invoke our commutator estimate, we have to pass from Bn(- -<■) to the usual 
paraproduct, which can be done using the following commutator lemma. Here we write 


B N tp = Bn( tp) = B N (ip, 1) = B n ( 1, ip). 

Lemma 8.15. Let p satisfy (Hh). Let a < 1, (3 G R, and ip G C € CL , Vn^P € . Then 

\\B N (ip^V N 'lp) - B N (V N 1p)\\ a+l3 < IM|a||7V0ll/3- 

Proof. By spectral support properties it suffices to control 


\(B N (Sj-np, AjVn'iP) ~ S j -npB N (A j V N ‘f 1 )) 0 )| 

= / (Sj-np(x + sy) - S j -iip(x))A j V N 'ip(x + ez)p(dy, dz) 
J R 2 

<e2^-<*-V\\ip\\ a \\V N 3fy. 


But now A :Vn = 0 unless 2- ? < N, and therefore we may estimate e 2 3 < 1. 


□ 
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Lemma 8.16. Let /j satisfy (Hh)- Let a + fi + 7 > 0, /3 + 7 < 0, assume that a € (0, 1), 
and let € ^ a , Vn if £ X £ ^ 7 - Define the operator 

A n X ( ( p) '■= J (nAr(niv(</? ■+ T- £y V N VO) o T- ez x) - V N {((p -< T- ey V N ifi) o T- £Z x))n{dy, dz). 

TTien /or all 5 € [0, a + (5 + 7 ) 

||njvSjv(IIjv(v? -< TVVO 0 x) - V N C(<p, V N if, x) - Vn{^B n {V n ^ 0 x)) II/ 3+7 

< AT-' 5 log(AT) 2 ||^|U||^ A r^|| /3 ||TllT+ll^ X llz. ( ^,^+- ) II^IU, 

where C is the commutator of Remark 2-4 and L(U,V) denotes the space of bounded 
operators between the Banach spaces U and V, equipped with the operator norm. 

Proof. We decompose the difference as follows: 


||riiv5Ar(nAr(^ -< V N ‘ VO ° X) ~ V N C(<p, V N fi>, x) - VN(<pB N (VNi> 0 x)) 11/3+7 


< 


-N (^Bj\f(Ylj\f(^ip -< 'P/vVO 0 x) - j T-eyV WV0) o T - ez xh{dy,dz)^ 

+ Pjv X (^) 11/3+7-5+ J\\V N (C(ip, T-eyV N fi>,T- £Z x) - C((p,V N ip,x))b+ytJ>(dy,dz). 


For the first term on the right hand side the same arguments as in the proof of Lemma 8.15 
show that for all 5 € [ 0 , 1 ] with a + /3 + 7 — <5 > 0 


n 


N^B N (Il N (<p-<VNfj))ox) - J{Dn{pP T- £ yVNfi>))°T- £Z x)h(dy,dz)) 

<N- s \og(N)^\ur N ^y\\xh 


0+/3+7—(5 


(the factor log(N) 2 is due to the operator Iljv which appears twice, see Lemma 8.8 and 
Remark 8.9). The second term is trivial to bound, and for the last term we simply use 
that 


T_ 


u - u|U_ 5 < N 5 \y\ 


e y 


whenever k € M, u € < io K , and 5 € [ 0 , 1 ]. 


□ 


8.2 Paracontrolled analysis of the discrete equation 

We now have all tools at hand that are required to describe the paracontrolled structure 
of the solution un to equation (74) which as we recall is given by 

J^nun = D n II n B n (un,un ) + D n Vn£,, mjv(0) = VnUq . 

We set 

Xjv(^) = (x n (0,x%(0,Xn(0,x}(0,Xn(0,b n (q n ox n )(0), 
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where 


&X N (g) = D N V N ^ 

= b n u n b n (x n (0,x n (0), 

&xl(Z) = D N n N B N (x N (Z),xV(Z)), 

X’xliO = B N u N B N (xl(0°x N (0), 

& X%(0 = D^n n b n (x%{Z),xV(Z)), 

^Qn(0 = B n B n (X n (£),1), 

all with zero initial conditions except X/v(£) for which we choose the “stationary” initial 
condition 

*Jv(£)(0)= [ 2 (B)B N V N tfis)ds. 

As in Section 6.2 we fix a £ (1/3,1/2) and fi € (1 — a, 2a). We will consider initial 
conditions that converge in which means that we expect the solutions (un) to 

converge in Ct^~^ whenever T > 0 and the stochastic data (X/y) converges in an 
appropriate sense. However, we have to be careful because while by now we know that no 
blow up can occur for solutions to the continuous Burgers equation, this is not so obvious 
in the discrete setting. So let Q be a cemetery state and denote by (ujsr) the unique 
solution to (74), defined up to the blow up time Tfi = inf{f ^ 0 : ||ujv(i)IU°° = °°} an d 
extended via rqv|[r^,oo) = C Note that the spectrum of un is contained in [-N/2, N/2\, 
and therefore all norms for 7 € M are equivalent to its L°° norm, so that we could 
as well have defined T // as the blow up time of ||ujv(i)ll-/ 3 - 

Theorem 8.17. Let (Xjv) be as in (85) and assume that the sequence is uniformly 
bounded in C T ^ a ~ l x C T ^ 2q_1 x if “ x x & 2a x C T ^ 2a ~ l for all T > 0, and 
converges to 

(X, X v , X^ + 2 cQ, X* 5 + cQ v + 2 cQ q 0 x , A v , Q o X + c) 

in C T ^ a ~ l x Ct^ 20 " 1 x x x a x C([5,T\^ 2a ~ l ) for all 0 < 5 < T, where 
c € M, 

X= (X,A V ,xV\x V ,QoJ) G X rbe , 

and 

^Q v = BX V , L£Q QoX = D(Qol), 

both withO initial condition. Assume also that the operator A^ = AJj N,Xn in Lemma 8.16 
converges to 0 inCT{L(fio Q ,‘^ 2 "~ 1 )) for allT > 0, a € (1/3, a). Finally, let (Vnu X ), uq € 
and assume that limjv||7 :, A^^fQ , — 11-/3= 0. 

Let (un) be the solution to (74). Then lmpy d-p(v,N, u) = 0, where u G f^beX ( see 
Definition 6.9) is the unique solution to 

2 z?u = Du 2 + 4cD« + D£, u(0) = uq. (77) 

Remark 8.18. According to Remark 7.6, equation (77) has a unique solution in ^beX 
which does not blow up. In particular lim/y Tfi = 00 , even if for fixed N we cannot 
exclude the possibility of a blow up. 
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Proof. Throughout the proof we fix a € (1/3, a) and we define 75 = (f3 + 5 )/2 whenever 
<5^0. We would like to perform a paracontrolled analysis of the equation, working in 
spaces modeled after the of Definition 6.9. For that purpose, we decompose the 
nonlinearity as follows: 

B N U N B N (u N , u n ) = ^ N (Xl + 2Xl + *Xn + x n) + 4D N IL N B N (X% P X N ) 

+ 4Dtv n^v Bn{x'n y Xn) + 20^11^5^(u^, Xn) 

+ 2Dn^-nBn(X^, 2x'n + Un) 

+ D n IIa:S v(2 x'n + u*n, 2x'n + Un)- 

The term T)nB.nBn{un, Xn) can be further decomposed as 

D n HnBn{un, Xn) = D n H N Bn{un -< X N ) + D n B-nBn(un y X N ) 

+ D n BnB n (un o X n ), 

and the critical term is of course DnBnBn(un 0 Xn). Using Lemma 8.13, Lemma 8.15, 
Lemma 8.7, and Remark 8.9 we have for all T > 0 

\\T>nB.nBn((2un + 4 x'n) -< Xn) — Un({2un + 4 x'n) -« Bn{P>nXn))\\m^^~ 2 

^ (II^TvIl <gt°L ’“(T) + ll^ivll^p)ll^ 7V H"~l' 

However, the Fourier cutoff operator n^r does not commute with the paraproduct (at 
least not allowing for bounds that are uniform in N), and in particular is not para¬ 
controlled. Rather, we have 

u< n = n n(u'n Qn) + Un 

with 

u' N = 2U Q N + 4Xl € > a (T), A G J?Z 2 *' m (T). 

This means that we need an additional ingredient beyond the paracontrolled tools in 
order to control the term DnBnBn{Bn{u' n y Qn) 0 Xn), and it is here that we need 
our assumption on the operator An- Under this assumption we can apply Lemma 8.16 
to write 

D7vn7v-R/v(njv(u(y y Qn) 0 x n) = Rn + VnC(u' n , Qn, X n ) + Vn{u' n B n {Qn 0 X n )), 

(78) 

where Rat is a term that converges to zero in 2 if u’ N stays uniformly bounded 

in Jzf/o *’ a {T). Denote now 

(X, X v , X*, X V , 77 , Q) = lim (X N ,X%, X%,x}, X%, B n {Qn ° X n ),Qn )• 
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Based on the above representation of the nonlinearity, it is not difficult to repeat the 
arguments of Section 6.2 in order to show that d—p(uN,u ) converges to 0, where 

u = A + A v + 2X* + u Q , u Q = v! ^Q + u\ 
with (vP, u', ifi) = limjv (u^,u r N ,u^ N ), where u(0 ) = no and 

2 ¥u Q = 2^A V + 4SfX* + 2B{{2X* + u Q ) 4) + 2B((2X* + u Q ) >- X) 

+ 2D(n # oi)| 2D((n' -« Q - v! < Q) o X) + 2BC(u', Q, X) + 2D {u'fj) 

+ 2B{X W {2X* + u Q )) + D(2X* + u Q f. (79) 

We also have 

v! = 2u Q + AX*. 

The fact that (Bn(Qn o Ajv)) converges not uniformly but only uniformly on intervals 
[S,T] for 5 > 0 poses no problem, because at the same time the sequence is uniformly 
bounded, so that given K>0we can fix a small 5 > 0 with supjv||u/v(<5) — Vn^o ll-/ 3 < 
and then use the uniform convergence of the data on [<5, T] and the convergence of 

(' P N Uo ) to u o- 

Now observe that 

u = X + X v + 2X* + u Q 
with vP — ifi = 2X* — 2X* = AcQ, and moreover 

u' = 2 u Q + AX* + 8 cQ = 2 u Q + AX*. 

Plugging this as well as the specific form of the X T into (79), we get 

= 4cDA + 2 ^A V + 4 j£f A^ + 4cDA v + 8 cD(Q o X) + 2D((2X* + u Q ) -< X) 

+ 2B({2X* + u Q ) yX) + 2B(J o X) + 2B((u Q - v! -< Q) o X) 

+ 2D C(vf, Q, X) + 2D[u\Q o X + c)) + 2B{X V (2X* + u Q )) + D(2A* + u Q ) 2 . 

Since 

8 cD(Q o X) + 2D o X) + 2D ((«' -« Q - v! X Q) o X) + 2DC(u', Q, X) + 2B{u'(Q o A)) 
= 8 cD(Q o X) + 2D (u Q o X) = 2D (u Q o A), 

we end up with 

2 z? u Q = Du 2 - 2 z?A - 2 z^A v - 2££X* + 4cDA + 4cDA v + 2cDh' 

= Du 2 - ££X - 2 z?A' v - 2Jgf X* + 4cDu, 

which completes the proof. □ 

It remains to study the convergence of the discrete stochastic data, which will be 
done in Section 10 below. 
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9 The stochastic driving terms 

In this section we study the random fields 

X, A v ,A* a\ A V , QoX 

which appear in the definition of X G A r b e - Our main results are the following two 
theorems, whose proofs will cover the next subsections. 

Theorem 9.1. Let £ be a space-time white noise onlxl and define 

Jz?A = D£, 

T£X y = D(A 2 ), 

Jz?A^ = D(ll v ), 

J^A^ = D^ol), 

= D(A V A V ), 

££Q = DA, 

all with zero initial condition except X for which we choose the stationary initial condi¬ 
tion 

A(0) = f P_ s D£(s)ds. 

Then almost surely X = (A, A v , A^, a\ A^, Q o A) € A r b e - If <P- R —X M is a measur¬ 
able, bounded, even function of compact support, such that 99 ( 0 ) = 1 and p is continuous 
in a neighborhood of 0, and if 

fe = = <p{e D)f 

(here & denotes the spatial Fourier transform) and X e = ©rbeC^e)? then for all T,p > 0 

linr o E[||X-X £ ||^ be(T) ]=0. 

Similarly, ift e = (X e , A £ v , x}. x), xf, Q e o X e ) for X £ = ip{eD)X and 

i?Ay = D((A e ) 2 ), 

T£X*f = D(A^(eD)A v ), 

J?x) = D(</?(eD)A^ o X £ ), 

5Txf = D((^(eD)A v v 9 (eD)A v ), 

T£Q e = DA £ , 

all with zero initial conditions, then for all T,p > 0 

limE[||X-i.||$ w 7 , ) ]=0. 
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Remark 9.2. The theorem would be easier to formulate if we assumed p to be contin¬ 
uous, of compact support, and with <^(0) = 1. The reason why we chose the complicated 
formulation above is that we do not want to exclude the function p(x) = (x). 

Theorem 9.3. Let £ be a space-time white noise on R x T. Then there exists an 
element Y € Tkpz such that for every even, compactly supported function p € C' 1 (M,M) 
with p(0) = 1 there are diverging constants c'f and cY for which 

limE[||Y-Y«||5 ta(T) ]=0 

for all T,p > 0, where Y e = @k P z( ( ^(^D)^, , cY). Moreover, we have 

* = h Jy (x)ix - 

Remark 9.4. We will only worry about the construction of X and Y. The convergence 
result then follows easily from the dominated convergence theorem, because since p is an 
even function all the symmetries in the kernels that we will use below also hold for the 
kernels corresponding to X s , X £ , Y £ . 

9.1 Kernels 

We can represent the white noise in terms of its spatial Fourier transform. More precisely, 
let E = Z \ {0} and let IT be a complex valued centered Gaussian process onKxE 
defined by the covariance 

E f f(ri)W(dr]) f g(rf)W(drf) = (2vr) _1 f g(m)f(jl-i)dVn 

_JRxE JRxE J JM.xE 

where i] a = ( s a , k a ), s- a = s a , k- a = —k a and the measure dr/ a = ds a dfc a is the product 
of the Lebesgue measure ds a on R and of the counting measure dk a on E. The functions 
f,g are complex valued and in L 2 ( R x E). Then the process f(p) = W{^p), where 
p € T 2 (RxT) and LX denotes the spatial Fourier transform, is a white noise on Lg(RxT), 
the space of all L 2 functions p with f T p(x, y)dy = 0 for almost all x. 

Convention: To eliminate many constants of the type (27r) p in the following cal¬ 
culations, let us rather work with y/2nW, which we denote by the same symbol W. Of 
course all qualitative results that we prove for this transformed noise stay true for the 
original noise, and we only have to pay attention in Theorem 9.3 to get the constant 
right. 

The process X then has the following representation as an integral 
X(t,x)= [ e ikx H t . s (-k)W( dr?), 

JRxE 

where 7] = (s, k) € R x E and 

h t (k) = e~ kH l t>0 , H t (k ) = ikh t (k). 
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This means that Ht- s (—k) = —Ht- S (k), and it will simplify the notation if we work 
with W = —W and £ = — £, which of course have the same distribution as W and £ and 
for which 

X(t,x)= f e ikx H t _ s (k)W( dr,). 

JrxE 

The space Fourier transform X(t, k ) = X t (k) of X(t, •) reads 

X(t,k)= [ e ikx H t _ s (k)W k (ds), 

Jr 


where W k '(ds) = f E S k _ k 'W(d-sdk) is just a countable family of complex time white noises 
satisfying Wfc(ds)* = W_ k {ds) and E[RV(ds)H4(d,s')] = 5 k ^ k /5(s — s')dsds'. 

Note that if s ^ t 


/ H s _ a (k)H t _ a (-k)da = 

Jr 

from where we read the covariance of X: 
E[X(t,x)X(s,y)]=E 


D —k 2 \t— s\ 


( 80 ) 




/RxE 


t- Sl (ki)W(dm) [ e ik ^H s _ S2 {k 2 )W{dm) 

JrxE 

= [ dk x e lkl{x -y> [ (-feudal 

J e Jr 

—kf\t—s\ 


= / e 
Je 


,iki(x—y) £_ 


1 


dh = -p\t- s \{x - y), 


where pt{x) is the kernel of the heat semigroup: Ptf = Pt * f = fjPt( • — y)f(y)dy. In 
Fourier space we have 

— k 2 \t—s\ 

E[X t (k)X s (k')} = 5 k+k , =0 C — - 

as expected. 

These notations and preliminary results will be useful below in relation to the repre¬ 
sentation of elements in the chaos of W and the related Gaussian computations. Recall 
that X* = X and X( T1T2 ) = B(X T \X T *). Then 


X T (t,x)= / G T (t, x, r] T ) TT W(drji) 

J(RxE)" i=1 


where n = d(r) +1, r/ T = = (rj i, G (M x E) n and dr/ T = dr/i... n = dry ■ ■ ■ dp n . 

Here we mean that each of the X T is a polynomial in the Gaussian variables W(drji), 
and in the next section we study how these polynomials decompose into the chaoses of 
W. For the moment we are interested in the analysis of the kernels G T involved in this 
representation. These kernels are defined recursively by 

G 9 (t,x,r,) =e ikx H t _ s (k), 
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and then 


G (T 1 T 2 \t, x, r? (TlT2 )) = B(G T1 (•, -,r] Tl ), G T2 (•, r? T2 ))(t, x) 

= [ daDPt- a (G T1 (a, ■,rj T 1 )G T 2 {a, •, r/ T2 ))(x). 

Jo 

In the first few cases this gives 

G v (t,x,rj V2 ) = f daDPt- a {G*(a, r]i)G*(a, i] 2 ))(x) 

Jo 

= e ik M x f daH t -„(k [ 12 ] )H t 7 - ai (k 1 )H t 7 - Sa (k 2 ), 

Jo 

where we set h = k\ + • • • + & n , and 

G^(t,2,77123) = f daBP t _ a (G v {a,-,r] 12 )G*(a r ,V3))(x) 

Jo 

= e ik W x [ daH t . a (k [123] ) (r du H a — r 7 r[k^i 2 ^H IJ '— Sl (^ki^H a '— S 2 [k 2 ) \ H a — S 3 [k^). 


( 81 ) 


In both cases, the kernel has the factorized form 

G T (t,x, VT ) = e ik M x H T (t,rj T ), 


where we further denote /q T ] = k\ i... n i =&! + ••• + k n , and this factorization holds for all 
G T . In fact, it is easy to show inductively that 

G (T1T2) (f,x,7/ (riT2) ) = e lk M x f daH t _ a (k [r] )H^(a^ T1 )H X2 (a, VT2 ), 

Jo 

from where we read 

G^(t, x, 773234) = e lk J 234 ^ x j dcrH t _ a (k[ 12M] ) f da'H a _ a ' (k [123] )x 


X [Jo d<J H v'-v" ( y k [i2})H <J »-s 1 {ki)H (jll — S2 (^ 2 ) J H a r — S 3 (ks)H a - S4 (* 4 ) 

and thus 

G^(t, 2, 773234) = e lfe [ 1234 P [ dcrl : /i_ cr (fc [ i 234 ] )ipo (^[123] 5 ^4) / 


X Wo dCr [ 1 2 ]) s 1 (_ s2 ( ) J _ S3 (fc 3 ) ^ -S 4 (*4), 

where we recall that ijj 0 (k,£) = Pi(k)pj(£)- Similarly, we have 


G V (t, 2, 773234) = e* fc [ 1234 l X / da da da" i 7 i-a(% 234 ])#a-<r^[ 12 ])#a-cr"(% 4 ])>< 

JO Jo Jo 

X JJcr' — si 1 )Hcr'— S2 0^2) H a "_ S3 (fc 3 )-fI (7 »_ S4 (fel) . 
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9.2 Chaos decomposition and diagrammatic representation 

The representation 


X T (t,x)= / G T (t,x,r] T )T\W(drn) 

J(RxE)« i=1 

is not very useful for the analysis of the properties of the random fields X. It is more 
meaningful to separate the components in different chaoses. Denote by 


/ f (m-n)W (dr]i... n ) 

J(RxE) n 

a generic element of the n-th chaos of the white noise W on (M x E). We find convenient 
not to symmetrize the kernels in the chaos decomposition. If we follow this convention 
then we should recall that the variance of the chaos elements will be given by 


E 


lxE) n 


f(m-n)W(dr]i... r 


CTg5n d(R x E)» 


where the sum runs over all the permutations of {1 ,,n} and where we introduce the 
notation rj-i = (si,— k±) to describe the contraction of the Gaussian variables. By the 
Cauchy-Schwarz inequality 


E 


lxE) r ' 


/(??!.. ■n)W(drii... r 


^ (n!) 


lxE) r ‘ 


\f(vi-n)\ 2 drn-r. 


so that for the purpose of bounding the variance of the chaoses it is enough to bound 
the L 2 norm of the unsymmetrized kernels. The general explicit formula for the chaos 
decomposition of a polynomial 


ixEy 


f(Vi-n)Y[W(dr}i 


2=1 


is given by 


ixEy 


f(m-n)Yl w ( dr k) 


i= 1 


E 

e=o ' 


lxE) k 


fo(r)i-e)W(dr) 1 ... e ) 


with fe{r]i...£) = 0 if n — i is odd. To give the expression for fp with n — i even, let 
us introduce some notation: We write S(£,n) for the shuffle permutations of {1,... ,n} 
which leave the order of the first £ and the last n — t terms intact. We also write 


7]Uv = T)i_£ U Vi ...m = (771, ...,r)e,v 1, .. .,v m ) 
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for the concatenation of r] and v. So if n — £ = 2m for some m, then 

ft (Vi...e) = 2 _m V V /(c r [ 7 ?L .. £ Ur( 7 ^ + i)... (£+m )Uf(? 7 _(^ + i)..._(^ +m )))])dr 7 (£+1) 

aeS(£,n) re5(m,2m), *^( Rx£: ) m 
feS(m) 

where cr(r/i... n ) = J7 <T ( 1 )... (T ( n ). This is just a formal manner of writing that we sum over 
all ways in which 2m of the 2m + t variables can be paired and integrated out. For 
example, we have 

W {(hj-i) = W (dr/ 12 ) + < 5(771 + 7?_ 2 )d77i dr/ 2 - 

In general we will denote by G T ( the kernel of the 7-th chaos arising from the decompo¬ 
sition of X T : 

n . 

^(M) = v / GJ(t,x,r] 1 ... i )W(dri 1 ... i ). 

(= 0 J(RxeY 

Terms X T of odd degree have zero mean by construction while the terms of even degree 
have zero mean due to the fact that if d(r) = 2n we have 

E[X T (t,x)\ =2~ n ^2 [ G T (t,a;,r(?7i... n Uf(77(_i)...(_ n ))))d?7i... n . 

re<S(n,2n) (® x ®) n 
feSn 

But now = k\ + • • • + k n — kf^ • • • — kft n ) = 0 and we always have 

G T (t,x, 7 ) 1 ... 2 n) oc k\ i...( 2 n)l> which implies that 

G T (t , x , T(rji... n U f(? 7 ( _ 1) ... ( _ n) ))) = 0. 

This is a special simplification of considering the stochastic Burgers equation instead of 
the KPZ equation. Later we will study the kernel functions for the KPZ equation to 
understand some subtle cancellations which appear in the terms belonging to the 0-th 
chaos. 

Applying these considerations to the first nontrivial case given by A v , we obtain 

X v (t,x)= I G v (t, 2 , 7712 ) IF (df/idr/ 2 ) + Go(t,x) 

J(RxE ) 2 

with 

Go(t,x)= f G v {t,x 1 r) 1( _ 1) )dr ]1 . 

J(RxE ) 2 

But as already remarked 

G v (t,x,r) 1{ _ 1) ) = e^C-DP H t ^, T (kn/i')i)H a - Sl (ki)H ( j_ S2 (k-i)dcT = 0 

Jo 

since Ht- a ( 0) = 0. Consider the next term 

X y '(t,x)= f G X '(t,2,77i23)IF(d77id?72d?73) + / G'f (t,x,rn)W(dr)i). 

J(RxE ) 3 JRxE 
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In this case we have three possible contractions contributing to , which result in 


Gi(t,x,rfr) = / (G^(t,x,r 7 12( _ 2 )) + G^(t,x,r/ 21( _ 2) ) + G Xr (t,x, 77 2 (- 2 )i))dr 72 . 


/RxE 


Y 

But note that G v (t, x, rj 2(_2)i) = 0 since, as above, this kernel is proportional to fc 2( _ 2 ) = 

Y Y 

0 . Moreover, by symmetry G v (f, 2, 7 / 12 (_ 2 )) = G v (f, 2 , 7 / 21(^2)) and we remain with 


Gi(t,x,rj 1) = 2 G yj (t,x,r ] i) = 2 x, r/ 12 ( _ 2 ) )dr/ 2 , 

JRxE 

where we introduced the intuitive notation G^(f, 2,7/1) which is useful to keep track 
graphically of the Wick contractions of the kernels G T by representing them as arcs 
between leaves of the binary tree. 

Now an easy computation gives 

rt r& 


G V (f, 2,7/1) = e 


— Jk ix 


dcr / da'H t _ r 7 (ki)H CT ,_ sl (k 1 )V' <? {a - a', /c 4 ), 


where 


V' v (a,k 1 )= I H a (k [ 12 ] )H r 7 _ S 2 (k 2 )H^ S 2 (k. 2 )d m = I dk 2 H a {k [12] )- 


=-M k l 


We call the functions V,£ vertex functions. They are useful to compare the behavior of 
different kernels. 

By similar arguments we can establish the decomposition for the last two terms, that 
is 


X^(t, x) = 


ixE ) 3 


G^it, x, 7/1234)^(d 7 /i 234 ) + 


IxE ) 2 


G}(t,x,r h 2 )W(dr 112 ) 


and 


with 


X xy (t,x)= [ G Xy (t, 2 , 7 /i 234 )W(d 7 /i 234 ) + I GY(t,x,ij 12 )W(dr]i 2 

J(RxE ) 3 J (RxE ) 2 


G%(t,x,rji 2 ) — f (G^(t, 2, ?/ 123 (_ 3 )) + 2 G^(t, 2, ?/ 13 (_ 3 ) 2 ) + 2 G^(t, 2, 7 / 132 (_ 3 )))d 7/ 3 
JRxE 


= (G^(i, 2, t/i 2 ) + 2 G^(t, 2,7/12) + 2 G^(t, 2,7/12)) 


and 


Gf ( t , 2, 7 /i 2 ) =4 / G V ( t , 2, 7 /i 3 (_ 3 ) 2 )d 7/ 3 = 4 G V (t, 2,7/12) 


,v 


■V/ 


/ExS 


Here the contributions associated to G H (f, 2, t/i 2 ) and G v (t,x,r/i 2 ) are “reducible” since 
they can be conveniently factorized as follows: 






G v {t,x,r]i 2 ) = / G s (t, 2 ,?/ 123 ( _ 3 ) )d 7/ 3 

JRxE 
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fdo f do' r d o''H t _ a {k [l 2 ] )H a ,_ a ,,{k [V 2 ] )H all _ sl {k 1 ) 
Jo Jo Jo 


x H a ff_ S 2 (k 2 )V^(o - o',k[ V 2}) 


V V (o,ki)= / 'i/j 0 (k[ 12 ],k-2)H ( 7 (k[i2])H a - S 2 (k2)H^ S 2 (k-2)dri2 


= / dk 2 Mk[i2],k-2)H a (k [1 2 ] ) 


e ~W\ k 2 


v f \ 

G v (t,x,rj 12 ) = / G°(t, x, 77 13 ( - 3 ) 2 )df ?3 

JRxE 


= e ik M x Mki,k 2 ) / d oH t _ a {k [l 2 ] )H a _ S 2 {k 2 ) / da' / da"iW(fci) 

JO Jo Jo 

xH^^ih^io' -o"M) 

= e ik M x Mki,k 2 ) [ doH t _ a (k [ 12 ] )H a _ S 2 (k2)e- ik ^G V (o } x, m ). 

Jo 

On the other side, the term G^(t, x, 7712) cannot be reduced to a form involving or 

\Q 

V °, and instead we have for it 

G^(t,x,rj 12 ) = / G^(t, x, ??i32(_3) )dr)3 

JRxE 

rt pa pa' 


— e *fc[i2]a; 


pt pa pa 

/ dcr / da / da"ift_ 0 -(A:[i2])i^ 0 -' —S 2 (^ 2 )Ha"—si (k\)x 

Jo Jo Jo 


x V v (a — 0,0 — o', k\ 2 ) 


v (o - o' ,0 - o",k 12 ) = / dri3'ip 0 (k [ 13 2],k- 3 )H a ^ S 3 (k- 3 )H a _ (Tl (k [ 13 2])x 

JRxE 

X Ha'-a"(k[ 13 })H a »- S 3 (k 3 ) 

f e - k lW~a"\ 

= dk 3 ip 0 (k [123] , k- 3 )H a _ a ,(k [ 123 ])H al - a " (k [13] )--- 

Similarly we have for G^ 


G*(t,x,ri 12)= / G v (f, x, r/ 13( _3 )2 )d773 
^Rxi? 


_ e *fc[i2]a; 


/*£ pa pa 

do da / da"i 7 i_ (T (A:[i 2 ] )H a r_ Sl (ki)H a "_ S2 (k 2 )x 
Jo Jo Jo 
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x v' <} '(a — a , a — a", k\ 2 ) 


with a vertex function which we choose to write in symmetrized form: 

V (a-a ,a -a ,k 12 ) = - / (b] 3 H a ^(k [13 ])H a _ a »(k[ 2 (- 3 )])H a ^ S3 (k 3 )H^_ S3 (k- 3 ) 

^ JRxE 


1 

+ 2 


'IxB 


drjsHfj—a" (k^i^)H(j— a i 3)] )H a i— S3 {k^)H a ff _ S3 (^-3) 


/ d/csi^-o-/ (fe[ 13 ] )H a _ a n (h 2(_3)] ) 

Je 


e -kl\a'-a"\ 


+ - / dk?,H a _ a ii(k\ i i3\)H a _ a i (fc[ 2 (_3)]) 


g-fcllo-'-o-"! 


9.3 Feynman diagrams 

While here we only need to analyze few of the random fields X T , we think it useful to give 
a general perspective on their structure and in particular on the structure of the kernels 
G T and on the Wick contractions. Doing this will build a link with the theoretical physics 
methods and with quantum field theory (QFT), in particular with the Martin-Siggia- 
Rose response field formalism which has been used in the QFT analysis of the stochastic 
Burgers equations since the seminal work of Forster, Nelson, and Stephen [FNS77]. 

The explicit form of the kernels G T can be described in terms of Feynman diagrams 
and the associated rules. To each kernel G T we can associate a graph which is isomorphic 
to the tree r and this graph can be mapped with Feynman rules to the explicit functional 
form of G T . The algorithm goes as follows: consider r as a graph, where each edge and 
each internal vertex (i.e. not a leaf) are drawn as 

ww— and 

To the trees V, X', \ V we associate, respectively, the diagrams 




These diagrams correspond to kernels via the following rules: each internal vertex comes 
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with a time integration and a factor (ik) 


m 



k l 12] 



d(T. 


Each external wiggly line is associated to a variable r/i and a factor H a _ Si (ki), where 
a is the integration variable of the internal vertex to which the line is attached. Each 
response line is associated to a factor h a _ C7 i(k), where k is the moment carried by the 
line and a, a' are the time labels of the vertices to which it is attached: 


k 

(J ^■ CJ 


h,j—a' (&) • 


Note that these lines carry information about the casual propagation. Finally, the out¬ 
going line always carries a factor ht- a (k), where k is the outgoing momentum and is a 
the time label of the vertex to which the line is attached. For example: 


m m 

f R 2_ dcrdcr' (ifc[i23]) ^t-a 0[123]) (^[12]) K-cj’ 0[12]) X 

x H a /_ Sl Q-f— S 2 (^ 2 )^( 7 —S 3 (^ 3 )• 

k[ 123] 

Once given a diagram, the associated Wick contractions are obtained by all possible 
pairings of the wiggly lines. To each of these pairings we associate the correspond¬ 
ing correlation function of the Ornstein-Uhlenbeck process and an integration over the 
momentum variable carried by the line: 

* , [ ik £^l 

a WWWW cr J 2 

So for example we have 

Vi 

r r e -klW-<r'\ 

2 / dcrdcr Hf—(j(k \) / d/t2 ~ 2 )]si(^i)? 

J r2_ Je 1 

ki 
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which coincides with the expression obtained previously. Note that we also have to take 

into account the multiplicities of the different ways in which each graph can be obtained. 

\ vy 

The contractions arising from G v and G w result in the following set of diagrams: 



The diagrammatic representation makes pictorially evident what we already have 
remarked with explicit computations: G ^ and G^ are formed by the union of two graphs: 



while the kernel G® 
to that of . 


cannot be decomposed in such a way and it has a shape very similar 


9.4 Bounds on the vertex functions 

Up to now we have obtained explicit expressions for the kernels appearing in the Wick 
contractions. These kernels feature vertex functions. Let us start the analysis by bound¬ 
ing them. Consider for example the first non-trivial one, 

G V (t,x,r] 1 )= f do f do'H t - (T {ki)H (T ,_ Sl (ki)V V (o - a', ki). 

Jo Jo 

We have 

|G V (f,£,r?i)|^ f d of do' M t - a (ki)M a /_ sl (ki)\v'^ (o — o' ,k\)\, 

Jo Jo 

where M t (k) = \H t (k)\= |£:|exp(— k 2 t)tt^o- Since the integrand is positive, we can 
extend the domain of integration to obtain an upper bound: 

[ d of do'M t - C 7 (ki)M CT '- sl (ki)\V^(o — o' ,ki)\= Z(t,x,rii). 

J M J M 

The quantity Z(t,x,rji ) is given by a multiple convolution, and since we are interested 
in L 2 bounds of G^(t,x,r) i), we can pass to Fourier variables in time to decouple the 
various dependencies: 

[ \G^(t,x,r]i)\ 2 dr)i ^ f \Z(t, x, ??i)| 2 dr/i = [ \Z(t, x, 0i)| 2 d6»i, 

JrxE JwlxE JrxE 


78 



where 6 \ = (uji : k\) and 


Z(t,x, 9 i) = [ 

Jk 


= I e- l “ lSl Z{t,x,rn)d Sl = M(-6»i)M(-6»i)V V (-o;i,A;i) 


with 


Then 


V V (u Jl ,k 1 )= [ Ida. 

JR 


\Z(t,x, 6 i)\ z & 0 i ^ / |M(0 1 )M(0 1 )| 2 |V V (o;i,fei)|^ 1 


/Rx£; 


/RxS 


< 




/RxS 


yV, A^OIda) d<?i, 


where Q(0) = #(0) is the time Fourier transform of H : 

««) = - * 




+ k 2 


This computation hints to the fact that the relevant norm on the vertex functions is 
given by the supremum norm in the wave vectors of the L 1 norm in the time variable. 
Summing up, we have in this case 


\G V (t,x, m )\ 2 d m ^ I \Q( 9 1 )Q(e 1 )\ 2 \\V v (a,k 1 )\\l 1 d 9 1 . 


2 mt rQi 


!RxE 


/R xE 


For the other contraction kernels we can proceed similarly. Consider first 




x 


G v (t, x, 7712 ) = / G v (t,x,r] i 23 (- 3 ))dj ?3 


/RxE 


_ e *fc[i2]* 


da / dcr' / da // ifj_ (T (A;[i2]) Hcr'—cr" (k[i2])Ha"—s± (fei)x 


x H^_ S 2 (k 2 )V^{a - a',k[ 12 ]) 


for which we have 

G^(f,X,7 ?12 )| 2 d7 ?12 < / |Q(% 2 ])Q(%2])e(0l)Q(02)| 2 ||^( C T,A ;[12] )||i 1 d0 12 . 

J(KxE ) 2 

The next term is 

= e <fc [“l x ^ 0 (*:i,fc 2 ) f daH t _ a (k [ 12 ] )H a _ S 2 (k 2 )e- ik ^G V (a,x, m ), 

Jo 

for which 


IxE ) 2 


IxE ) 2 


x, ??i 2 )| 2 d??i 2 ^ 


\^o{ki, k 2 )Q( 8 ^ 2 ])Q( 9 i)Q( 9 i)Q( 9 2 )\ 2 \\V^(a, k\)\\ 2 L i d# 12 


IxE)' 
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And finally 


G*{t,x, m2 ) = e ik W x da da' da // ^_ (T (/c [12] )^_ S2 (fe 2 )^_ Sl (A;i): 
Jo Jo Jo 

x V^(a - a', a - a ", k\ 2 ), 


for which 


J(RxE) 2 
Similarly, we get 


| G^(t,x, m2 )\ 2 d Vl2 ^ [ \Q(e [12] )Q(e l )Q(6 2 )\^{a,a',k l 2)\\ 2 L i L i d0 12 . 


[ |G V (t,x,7 ?12 )| 2 dr/ 12 ^ / |g(% 2 ])Q( 0 i) ( 3(^)| 2 r V (^^^i2)|lii L i d 012 . 

J(RxE) 2 J(RxE) 2 CT CT ' 

For these computations to be useful it remains to obtain explicit bounds for the norms 
of the vertex functions. 


Lemma 9.5. For any e > 0 we have 


I da\V yj (a,k 1 )\+ I da\V*(a, < \ki\ e . 

J M J M 

Proof. Some care has to be exercised since a too bold bounding would fail to give a finite 
result. Indeed, a direct estimation would result in 


% 


da 


V V {a, ki) < f dk 2 [ da|fc [12 ]|e a( ' k2+k [^] ) ^ f dfc 2 J ^ 12] | = +oo 

JE Jo JE ^2 mi 2 l 


due to the logarithmic divergence at infinity (recall that f E stands for ^z\{o})- To 
overcome this problem, observe that 


F V (a,0) = 2 / d k 2 H a {k 2 ) 


e -M k 2 


= 0 


since the integrand is an odd function of k 2 . So we can write instead 

pq, . . . r .. e ~\ a \ k i 


V v (a,h) = 2 / d k 2 [H a (k [12] ) - H a (k 2 )\- 
Je 

and at this point it is easy to verify that 


da|F V (a,fci)|< \h\ e 


for arbitrarily small e > 0. Indeed, 

p p poo 

/ da|I/ V (a,fci)|< / dk 2 / da\H a (k [12] ) - H a (k 2 )\e- ak2 

Jm. J e J o 
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and a first order Taylor expansion gives 


I H a (k [l2] ) - H a (k 2 )\< | k \| / dre-^ 2 ^ 1 ) 2 ". 


Therefore, 


/*1 poo 

da\V V (a,fa)\ < |fci| / dr dfa dae~ c{k2+Tkl)2(r - k ^ 
: Jo Je Jo 

~ |fcll i dT / E (fe+rtW*r 


but now 


dfc 2 


< 


d/c 2 


< 1. 


7s (fa + rki) 2 + kl~ J E k% 

On the other side, the sum over T is bounded by the corresponding integral over M, so 
a change of variables gives 


dfa 


< 


dfa 


< 


1 


d£: 2 


< 


1 


Je (fa + rfa) 2 + ~ J R (fa + rk i) 2 + k\ ~ r\fa\ J R (. k 2 + l) 2 + k% ~ r|fci|' 

By interpolating these two bounds, we obtain 


do-|y V (cr,A:i)|< |fci| £ 




Y? 

for arbitrarily small e > 0. The same arguments also give the bound for f R da| V °(a, fa )|. 


□ 


Lemma 9.6. For any e > 0 we have 

[ V^(a,a',ki 2 ) dada'<\k [l2] \~ 1+£ . 

Jr 2 


Proof. We get 


V^(a, a’,fa 2 ) 


da da’ < 


< / dfe 3 


dAfc^o(^[ 123 ], k-3)H a (k [132] )H (J/ (fc [13] )e 

l%L23]ll%3]l 


-k%(a'+cr) 


dada' 


E (*% + k\ m )(kl + kf m ) 


</ [13] > 


1 


~ Je dki (kl + ^fl23]) 1/2 ( fc 3 + fc fl 3 ]) 1/2 
~ j E dh (k 2 + kf 12 3] ) 1 / 2 |fc 3 | ~ |fc[12]l ^ 


for arbitrarily small e > 0. 


□ 
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Lemma 9.7. For any e > 0 we have 


y V (cr,cr',fci 2 ) dcrdcr' < |fe [12 ]| 


— 1+£ 


Proof. We can estimate 


V^(a,a',k 12 ) 


dada' 


dada' 


f f dk 3 H a (fc [13] )H a / (k [{ _ 3)2] )e- fc 3 2 l-'-l 
J J e 

poo poo p 

/ dcr / da' / dfc 3 1^(fc [13] )(^[(-3)2])Ie- fc i CT ' 
JO Jo JE 

poo poo p 

+ d a / da' / dk 3 \H a+a \k [l3] )H a (k [{ _ m )\e~ k ^ c 

Jo Jo Je 


< / dk 3 


l%3]ll%-3)2]l 


+ / dk 3 - 


Ifyl3]ll%-3)2]l 


(fcg + —.^21)(^ni 3 1 + ^Y-3'l2l) Je (&3 + ^fl3])(^f(-3)2] + ^M.3l) 


v 3 T- n, [(-3)2]A n '[i3] t<v[(_3) 2 ] 

1 

E ^ {k 2 z + £f ( _ 3)2l ) 1/2 (£f 13l + tf ( - 3 )2]) 1/2 


*'[13]> 


< / d k 3 


v [(—3)2] ^ V n '[l3] ' 3)2] > 

1 

E ^ (*3 + k [13 ]) 1/2 ( fe [(-3)2] + fc [13]) 1/2 


+ / dfe 3 


< / dA: 3 - 


1 


/e I ^ 31 (^ 3 ] + ^ f (_ 3) 2 ])^ 2 
whenever e > 0 . 


< 


I ^[12] I 


—1+£ 


□ 


9.5 Regularity of the driving terms 

In this section we will determine the Besov regularity of the random fields X T . Below 
we will derive estimates of the form 

sup 2 ‘*nlr) q \t - sr 72 ^E[(AgX T (t, x) - A q X T (s, x)) 2 ] < 1 (82) 

q^0,xET,s,t^0 

for any 7 i(r), 7 2 (r) ^ 0 with 7 i(t)+ 7 2 (t) = 7 (r). Each A q X T {t 1 x) is a random variable 
with a finite chaos decomposition, so Gaussian hypercontractivity implies that 

sup 2™^ q \t - s\-™W 2 E[(A q X T (t, x ) - A q X T (s, x)f] < 1 

q^0,xET,s,t^0 

for any p Js 2 , and then 

sup 2^ q \t - s\~™W 2 n\\A q X T (t, •) - A q X r (s, ■ )||*} < 1. 

q^0,s,t^0 v ’ 


From here we derive that for all e > 0 


sup |t - s\-™W 2 E[\\X T (t, 0 - * T (s, Oil* 


s,t^0 


B, 


71(t)— eJ 
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< 1 , 


= sup|t- S |-^^/ 2 E[ Y, 2^-^ q \\\X T {t,-) - \X T {s,-)\\% m 

sJ^O 1 

and by the Besov embedding theorem we get 

II X T (t, •) - X(S, OIL,! W-2e< || X T (t, •) - X T ( S , -)|| 7 i W -e, 

00,00 u p,p 

where e > 0 can be chosen arbitrarily small since p is arbitrarily large. Thus an appli¬ 
cation of Kolmogorov’s continuity criterion gives 

J S 1 

whenever e > 0, 7 i(t)+ 72 (t) ^ 7 (r), andp, T > 0. This argument reduces the regularity 
problem to second moment estimations. 

Next, observe that every X T solves a parabolic equation started in 0 (except for 
r = • which is easy to treat by hand). Thus we get for 0 ^ s ^ t 

X T (t, •) - X T (s, •) = ( P t - r u T (r)dr- f P s - r u T (r)dr 
Jo Jo 

= (P t - S - 1 )X r (s) + [ P t _ r u T (r)dr. (83) 

J S 

We estimate the first term by 

|A,(P t _, - l)A r (s,x)|< \e< t - s ^ 2q - l||A (? X r (s,x)|< 1 1 - s\^/ 2 2™^\A q X T (s,x)\ 

whenever 72 (r) G [0,2]. The first estimate may appear rather formal, but it is not 
difficult to prove it rigorously, see for example Lemma 2.4 in [BCD11], So if we can 
show that 

E[\A q X^s,x)\ 2 }<2~ q ^ T \ 

for 7 (r) ^ 2 , then the estimate (82) follows for the first term on the right hand side 
of (83). For the second term, we will see below that when estimating E[| A q X T (t, x)| 2 ] 
we can extend the domain of integration of Pt~ r A q u T {r)dr until — 00 , and this gives 
us a factor 2~ 2<? . If instead we integrate only over the time interval [s, t], then we get 
an additional factor 1 — e c b-s) 2 2<3 ^ w j 1 i c ] 1 we can then treat as above. We will therefore 
only prove bounds of the form E[| A q X T (t, x)\ 2 } < 2~ 2qi ^ T \ 

Let us start by analyzing Q o X, whose kernel is given by 

GQ° x (t,x,T) 12 ) = e ik M x Mh,k 2 )H t _ sl (k 1 ) f d aH t ^(k 2 )H a _ S2 (k 2 ). 

Jo 

Now note that, by symmetry under the change of variables k\ —>• —k\ we have 


G\ 


QoX 


(t,x) = f G q oX (t, x, ?7 1( —!) )dryi = 0 
J(RxE) 
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since 


G <3 °- Y (t,x,r/ 1 (_ 1) ) = il) 0 (ki,-ki)H t - sl (ki) [ daH t _ a (-ki)H a ^ S2 (-ki) 

Jo 

and H t -a(—ki) = —H t - a {ki). So only the second chaos is involved in the following 
computation: 

E[(A g (QoI)(t,s)) 2 ] </ 2 d^( % 2 ])^(.i,. 2) 2 J J dada , H t _ a (k 2 )H t _ a/ (k 2 )^—^ - 

< J dk 12Pq (k [12] ) 2 Mh,k 2 ) 2 \k 2 \ 2 J* dada'e-^-^-"'^-^ 

< [ dk 12Pq (k [12] ) 2 Mh,k 2 ) 2 \k 2 \- 2 < [ ^ dkdk'p q (k) 2 Mk - k\ k') 2 \k'\~ 2 

J E ^ J E ^ 

< 2 ? V2 ,:_2i < 1. 

rsj / j rxj 

i>q 

Similarly we see that E[(A 5 (Q o X)(t, x) — A ? (QoI)(s,i)) 2 ] < 2 qK \t — s\ K / 2 whenever 
k € [0, 2], from where we get the required temporal regularity. 

Next we treat the X T . As we have seen in the case of the vertex functions it will be 
convenient to pass to Fourier variables. In doing so we will establish uniform bounds for 
the kernel functions ( G T ) r in terms of their stationary versions (r r ) r : that is the kernels 
which govern the statistics of the random fields X T (t, •) when t —>• +oo. 

Let recursively T* = G* and 

r (riT2) (f,x,77 (riT2) ) = f dsd x Pt- s (r T1 (s,-,ii T1 ),r T2 (s,-,7lT 2 ))(x). 

J —oo 

Like the G kernels they have the factorized form T T (t, x, r) T ) = e ik ^ x 'yj (?y T ). The advan¬ 
tage of the kernels T is that their Fourier transform Q in the time variables (si,..., s n ) 
is very simple. Letting 0 t = (, Ui,ki ) we have 

Q T (t, x, 0 T ) = f ds T e iu} ^T T (t, x, rj T ) = e ik M Y(0 r ), 

J(RxE) n 

where 

q'(0) = -^ 2 , q {T1T2) (0 {T1 r2) ) = q-(0 [{T1T2)] ) q ^(0 T1 )q X2 (0 T2 ). 

The kernels for the T terms bound the corresponding kernels for the G terms: 

| G^(t,x, V{Tin) )\= | H x (t, VT )\< W(vt)\= |r (riT2 )(t,x,7 7(Tir2) )|. 

This is true for r = • and by induction 

|G (rir2) (Lx,77 (riT2 ))| = |iL (T1T2) (L??( riT2 ))|= |fe[( ri r 2 )]l [ h- s (k [{TlT2]] )\H T1 (s,r] T1 )\\H T2 (s,7j T2 )\ds 

Jo 
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< l fe [(rir 2 )]l / h t -s(k [{TlT2) ])\g T s 1 {r] T1 )\\gl 2 (7j T2 )\ds 

< l fe [(rir 2 )]l I h-s(k[( T1 ^)])hl 1 (r] Tl )\\'y^(r] T2 )\ds 

J — OO 

= l7t (TlT2) (^r 2) )l=|r (T1 " 2) (t,X,7 7{riT2) )|. 

These bounds and the computation of the Fourier transform imply the following estima¬ 
tion for the L 2 norm of the kernels G uniformly in t ^ 0: 


\G T (t, x, fy r )| 2 ds r ^ / \T T {t,x,p T )\ 2 ds r ^ / \Q T (t,x,0 T )\ 2 du T% 


that is 


\G T (t,x,r] T )\ 2 ds T ^ / \q T (9 r )\ 2 duj T 


(84) 


This observation simplifies many computations of the moments of the A T ’s and gives 
estimates that are uniform in t ^ 0. Actually it shows that the statistics of the A T ’s are 
bounded by the statistics in the stationary state. 

Now, note first that due to the factorization (81) we have that 


A q G T {t,X,7lr) = p q (k [r ])G T {t,X,7Jr), 
so the Littlewood-Paley blocks of X T have the expression 

« n 

A q X T (t,x)= / p q (k [T] )G T (t,x,7j T )Y\w(d7ji) 

J(RxE )" i=1 

which we rewrite in terms of the chaos expansion as 

d(r) 

A q X T (t,x) = Y^ / Pq{k[r])G}(t,X,p 1 ...i)W(d7]l...e). 

^_q J (M.xE) 71 

By the orthogonality of the different chaoses and because of the bound (84) we have 

d{r) 

E[(A <? A r (t,x)) 2 ] < ^2 / Pq(k[ T ])\G T £ (t,x,p 1 ... e )\ 2 dp 1 ... e 

d (j) 

/ Pq(k[ T ]M(vi-t)\ 2 dvi-e- 

J(M.xE) n 

By proceeding recursively from the leaves to the root and using the bounds on the vertex 
functions that we already proved, the problem of the estimation of the above integrals 
is reduced to estimate at each step an integral of the form 



e\~ fi de, 


85 



where we have a joining of two leaves into a vertex, each leave carrying a factor propor¬ 
tional either to \9\~ a or \9\~P with a,f3 ^ 2 and where the length |(9| of the 9 variables 
is conveniently defined as 

|#|= |w| 1 / 2 +|A;|, 

so that the estimate |g(0)| |0| 1 holds for q = q*. 

Lemma 9.8. For this basic integral we have the estimate 

[ | o\- a \o'-e\-Pdd<\o'\- p , 

jRxE 

where p = a + (3 — 3 if a, f3 < 3 and a + (3 > 3, and where p = a — e for an arbitrarily 
small e > 0 if (3 ^ 3 and a € (0,/?]. Similarly, if /3 ^ 3 and a € (0, /3\, then 


if 0 [k, k! - k) 2 \9\~ a \9’ - 6\-Pd9 < \9'\- a+£ \k'\ 3 - 0 . 


/RxE 


Proof. Let t € N such that |0'|~ 2 e : 

\e\- a \e' - 6\-Pde < £ 2 —« [ 

/ ; -\n JRxE 


IRxE 


-0 ~2J 


d 9. 


i,j>0 


Then there are three possibilities, either i < i ~ j or i < j ~ £ or j < i ~ In the first 
case we bound 


JRxE 

and in the second one 


^|0j~2b|0|~2%|0'-0|~2.jd# ~ / l|0|~2*d# < 2 3t , 


/Kx_E 


JlxE 

and similarly in the third case 


l|0j~2b|0|~2%|0'-0|~2.jd# < / l|0|~2*d# < 2 3t , 


/Kx_E 


JlxE 

So if a + f3 > 3 we have 


l|0 , |~ 2 h|0|~2 i ,|0'-0|~2.)d# < 


/RxE 


1\e h 2id9 < 2 3j . 


f -0\~ P &e < ^2 2~ ai - (3j+Si + ^ 2 -a»-«+3i + 2 -°d-t3j+Zj 

< 2 —^(or-l-/3—3) _|_ 2 — /3t+(3—«)+t _|_ 2 —Q'^-l-(3—< ^—pt- 

where p can be chosen as announced and where we understand that (5)+ = e if 5 = 0. 

Let us get to the estimate for the integral with if a . Let i' ^ i be such that |&:'|~ 2 f 
and write 


\ifo(k,k' - k)\\9\~ a \9' -9\-P&9 


'RxE 
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E E 

JRxE 

J2 J2 [ l»'>^l|0f 

i,3>0i'<i/\j ^ RxE 






Now we have to consider again the three possibilities £ < i ~ j or i < j ~ £ or j < i ~ t. 
In the first and second case we bound the integral on the right hand side by 2 2 * +i , and 
in the third case by 2 2 - ? ~*~ i . Then we end up with 

[ \Mk,k' -k)\\6\- a \0' - 9\-Pd6 

JRxE 

< y ^ 2 _a * — /3j+2i+i' _|_ 2 _ai_ ^+2i+i' _|_ y ^ 2~“* _ /3j+ 2 J+*' 




i<j 




< 2 -dcH-/3-3) + 2~^+( 3 - a )+ £ +( 3 -")- £ ' + 2- to +( 3 -< s )+^K 3 -/ 3 M' 
with the same convention for (<5) + as above. 

We will also need the following simple observation: 


Lemma 9.9. Let a > 2. Then 


dw(M 1/2 +ifcir a < ifci 5 


Proof. We have 


/ dc^M^+lfel)- 0 = |fc|- a / dw(|w|A;| _2 | 1//2 +l) _ “ = |fc| 2_a / dwdwl^+l)- 4 *, 

J ]R */ IR «/IR 

and if a > 2 the integral on the right hand side is finite. □ 

Consider now A v . Combining the bound (84) with Lemma 9.8 and Lemma 9.9, we 
get 

E[(A q X v (t,x)) 2 ] < [ p q (k {1 2]fW [ i2 ] )q(0i)q(e 2 )\^e l2 

J (RxE ) 2 

= [ d9[12]d02Pq(k[l2]) 2 \q(0[l2])q{0[12] - O2)q{02)\ 2 

J(RxE ) 2 

< / d0 [12] d0 2/ 9,(A: [12] ) 2 | g (0 [12] )| 2 |0 [12] - 0 2 |- 2 |0 2 |- 2 

J(RxE ) 2 

< [ d% 2]P(? (/c [12] ) 2 | (? (% 2] )| 2 |0 [12] |- 1 < [ dk [12] Pq{k[1 ^ < 1. 

IlxE JE l«[12]| 

Y 

As far as A v is concerned, we have 


E[(A ? xVx)) 2 ]< 


2 V, ^ 

A?(%23]) G (t,x,rj 123 ) d?yi 23 
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+ [ Pq(k l) 2 clr /1 

Jm.xE 

< / P« ? (fc[123]) 2 k(^[123])g(^3)9(%2])9(^l)9(^2)| 2 d6'i23 

i(RxE ) 3 

+ [ Pq(h) 2 \q(0i)q{0i)\ 2 ( [\V V (a,ki)\da] d0i. 

JRxE \Jr / 

For the contraction term we already know that J k |F^(cj, fci)|dcr < |A:i| e < |$i| £ , so 

[ pq(h) 2 \q(0i)q(0i)\ 2 ( [\V^(a,ki)\da) d6»i < f p q (k 1 ) 2 \9i\~ A+2e d9 1 

JlxB ' tl ' JlxE 

< f p q (ki) 2 \ki\- 2+2£ dk\ < 2^ 2e ~ 1 ^ >q . 

JE 

The contribution of the third chaos can be estimated by 



P< ? (%234]) 2 |V’o(A;[234],^l)g , (6 , [1234])9(^l)g(^[234])9(^2)g(%4])g(^3)Q(^4)| 2 d01234 

% [ Pq{k [1 23]) 2 \ipo(k [ 23],k 1 )q(e [1 23])q{9 1 )q(e [23] )q(e 2 )q(9 3 )\ 2 \9 3 \^ 1 de 1 23 


P q {k[i 2 ]) 2 \A(k 2 , ki)q{9 [1 2 ] )q(9 1 )q(9 2 ){ 2 \92\ 2+£ d6»i 2 
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< 

r-^i 


ixE 


Pg(>i) 2 k(#i)| 2 N 2+£ \h\ 1+e d6»i < / p q (k 1 ) 2 \ki\ 2+£ \h\ 1+£ dki 


< 


For the contractions we have 


IxE ) 2 


P q {k m ) 2 \G*(t : x : in2)\ 2 din2 


< 


Pq(k [1 2]) 2 \q(0 [12] )q{d [V2] )q(e 1 )q{e 2 )\ 2 \\V^< y a,k [12] )\\ 2 Ll de 1 2 


< 


IxE ) 2 


IxE ) 2 


P q (h)' 2 \q(0i)\ 4 \0i ~ 0 2 \ 2 \9 2 \ 2 |/ci| 2 e d 6 »i 2 


f Pg(h) 2 \9i\ 4 |6»i| 1 \h\ 2£ ddi< [ p q {k\) 2 \ki\ 3+2£ dfei<2^ 
Jrxe Je 


— 2+2e) 


The following term is 

P q (k [12] ) 2 \G^(t,x,pi 2 )\ 2 dp 12 


IxE ) 2 


< 


2 iitA?/ 


IxE ) 2 


^(fe [ 12])>o(fel,fe2) 2 k(%2])9^l)9(0l)^2)r||^ B ( ( T,fc [12] )||i 1 d0 1 2 


< 


[ p q (k i)Vo(fci - k 2 M?\q(0i)\ 2 \<h - 0 2 \- 4 \e 2 \- 2 \h\ 2£ d9 12 

J(RxE ) 2 


< 

r^i 


/RxE 


p g (k 1 ) 2 \q(9 1 )\ 2 \9 1 \- 2+£ \k 1 r 1+2£ d9 1 < / p q (ki) 2 \ki\~ 3+3e dki < 2«< 


/RxE 


Next, we have 


Pq(k[i 2 ]) 2 \G^{t, x , r/i 2 )| 2 dr/i2 


IxE ) 2 


< 


Pg(%2]) 2 k(^[12])g(^l)9(^2)| 2 ||^(o-, 0- , ,fcl2)Hil i l ; d6'i2 


!xB) ; 


< 


PqihfWiM 2 ^ - e 2 \- 2 \e 2 \- 2 \k l \~ 2+2£ de l2 


IxE ) 2 


< 

r^j 


[ pq(h) 2 \0i\ 3 \h\ 2+2e d9 l < f p q {k 1 ) 2 \k 1 \ 3+2£ dk 1 <2 q( 

JrxE Je 


—2+2e) 


and therefore E[(A g X^(t, x)) 2 ] < 2 q ( 2 3e \ 
The last term is then 


< 


E[(A q X^(t,x)) 2 ] 

f P q (k[i 234] ) 2 \q{9 [1234] )q(9 [12] )q(9 1 )q(9 2 )q(9 [34] )q(9 3 )q(9 4 )\ 2 d9 1234 

J(lxB ) 4 

+ / Pq(k[i 2 ]) 2 \GY(t,x,p 12 )\ 2 dp 12 . 

J»xE) 2 


2 <?(—2+2s) 


-2+3e) 
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The first term on the right hand side can be bounded as follows: 


IxE) 4 


Pq{k [12 34]) 2 \q^[i234])q^[i2])q{0i)q^2)q{0[34])q^3)q^4)\ 2 dei234 


< 


< 


IxE) 4 


IxE) 4 


p q {k[ 12 ]) \q(0[i2])q{0i)q(0 2 )\ |6»i| \e 2 \ d e 12 
P q (k l ) 2 \q{9 l )\ 2 \9 l \^ +£ d9 l <2^ q . 


On the other side, the contraction term is given by 
/* 

Pq (*[12] ) 2 I GY (k X , 7/12) | 2 d ??12 


IxE ) 2 


< 


< 

r^j 


IxE)'- 


IxE) 2 


Pq(k [12] ) \q(O[i2])q(8i)q(02)\ 2 \\V xy (a,a',k 12 )\\ 2 LlLl d9 u 

G o' 

Pg{ki) 2 \q{0i)\ 2 \9i - 9 2 \~ 2 \e 2 \~ 2 \k 1 \~ 2+2£ ddi 2 


< 


[ p q (k 1 ) 2 \9 1 \- 3 \k 1 \~ 2+2s d9 1 < [ p q (ki) 2 \ki \ _3+2e dAq < 2«(- 2 + 2£ ), 

JRxE JE 

so we can conclude that E[(A g X^(i, x)) 2 ] < 2~ q< ' 2 ~ 2£ ' > . 


9.6 Divergences in the KPZ equation 

The data we still need to control for the KPZ equation is 

y,y v ,y*,y\y xy 

since q o X was already dealt with. The kernels for the chaos decomposition of these 
random fields are given by 

Gt(k) = l t ^oe~ tk \ G^ T1 ’ T2 \t,x,ri T ) = [ daP t - a (G T1 (a,-,7j T1 )G T2 (a, •,fy T2 ))(x), 

Jo 

so they enjoy similar estimates as the kernels G T and therefore all the chaos components 
different from the 0-th are under control. The only difference is the missing derivative 
which in the case of the X T is responsible for the fact that the constant component in 
the chaos expansion vanishes. The 0-th component it given by 

tc T = E[y T (t, x)\. 

Some of these expectations happen to be infinite which will force us to renormalize Y T 
by subtracting its mean. 
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For K v we have 

E[K v (t,x)] = f E[P t . a (X £ (a,-) 2 )(x)}da 
Jo 

= [ dcr f dine lkl< ~- 1)X e~ (t ~ CT)kl, ~- 1 '>(p(£k 1 )H CT _ sl (ki)^(£k- 1 )H CT _ sl (k- 1 ) 
Jo JRxE 

= f dcr [ drjnp 2 (£ki)H a _ sl (ki)H a ^ sl (k^i) = *- [ dkup 2 ^), 

Jo Jrxe ^ Je 


and since (p € C* 1 we get 


|(/ p(x)dx-sj dkip(Ek 1 )'} = (^ 2 (x) - ^ 2 (£fc))d£ 


7E 


1 /*£ 



D(<^ 2 )(eA: + Ax)xdxdA, 


which converges to D</? 2 (x)dx = 0 as e —)• 0. Now recall that here we are dealing with 
(27T) 1 / 2 times the white noise, so that the constant of Theorem 9.3 can be chosen as 

j = ikly ix)dx - 


The next term is T , which has mean zero since it belongs to the odd chaoses and 
thus = 0. 


What we want to show now is that a special symmetry of the equation induces 

vy y, 

cancellations which are responsible for the fact that while c v and c v are separately 
divergent, the particular combination 


YY 

is actually finite. If this is true, then we can renormalize Y H and Y w as announced 
in Theorem 9.3. In terms of Feynman diagrams (which have the same translation into 
kernels as for Burgers equation, except that for the outgoing line the factor ihp is 
suppressed), this quantity is given by 

/ = c v + 4c' = 2x 

because all other contractions vanish since they involve contractions of the topmost 
leaves (e.g. x, ?h(-i) 2 (- 2 )) = 0). Moreover, writing explicitly the remaining two 
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contributions and fixing the integration variables according to the following picture 



we get 


"U,2 
[ 12 ] 


1 = 2 [ d hdk 2 [da f da' / CT da // (zA: [12] )(-ifc [12] )e- |CT '- a '' l{ ^ +fc22) - {<T - CT ' )fc ^]- (<T - CT " )fc : 

J Jo Jo Jo 

+ 8 [ dhdk 2 f da ["da' r da"(ik [l2] )(ik 2 )e- {a '- a " )ik2 i+k 2 2 )-(<r-a')kf 12] -(a-a’’)^ 12] _ 

J Jo Jo Jo 

Now note that the second term can be symmetrized over a ', a" to get 
1 = 2 [ dhdk 2 f da ^ da' r d ^"(i/fc [12] )(-ifc [12] )e- |CT '- CT/,|(fc?+fci) - (<r - CT,)fc ^ ] - ( — 

J Jo Jo Jo 

+ 4 [ dhdk 2 f da [" da' ^ da // (^ [12] )(zA: 2 )e- |CT '- CT '' l(fe?+fc22) - {a - CT ' )fc [^]- (CT - <T '' )fci 

J Jo Jo Jo 

At this point the second term can still be symmetrized over k 2 . k\ to hnally get 
ft 


"U,2 
[ 12 ] 


[ 12 ], 


1 = 2 [ d hdk 2 I da r da' / drT , '(zfc [12] )(-ifc [12] )e- |fT, - CT,/|(fc?+fci) - (CT - fr ' )fc f^]-(—")fc [12] 

J Jo Jo Jo 

+ 2 [ dhdk 2 f da ^ da' f a da // (iA; [12] )(zfc [12] )e _l<T,_<T,,|{ ^ + ^ )_{<T_CT,)fe t L2 ] _ ^ <T_<T,, ^i 

J Jo Jo Jo 


"H.2 
[ 12 ] 


and conclude that 1 = 0. For the moment this computation is only formal since we did 
not take properly into account the regularization. Introducing the regularization cp on 
the noise, the integral to be considered is 

/ rt rcr pa 

dk\dk 2 ip 2 {£k\)ip 2 (ek 2 ) da da' / da "x 

Jo Jo Jo 


rt 


[ 12 ] 


/ pt pa pa 

dkidk 2 (p 2 (£h)(p 2 (£k[i 2 ]) da da' / d<j"(i/c[ 12 ])(?'A: 2 )x 

Jo Jo Jo 


x e -(°' , - 0 '")( fc l+ fc 2)-( <7 - 0 ' , ) fc fi2]-( 0 '- (T ") fc fi2] 

Here the symmetrization of a' and a" can still be performed, giving 

p pt pa pa 

Is = 2 dhd k 2 da da' da"(ik[ 12 ])[(p 2 hh)ip 2 (£k 2 )(-ik[ V2 ]) + 2(p 2 hh)ip 2 (£k[ V2 ])(ik 2 )} x 

J E 2 Jo Jo Jo 
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x e 


-|0-'-<T"|(fcf+fc|)-((T-(T , )fcn 2 ]_(0--<T")fc?i2] 

which is equivalent to 

4 = 4 / dhdk 2 [ da f da' [ da"(ikt 12 i)ip 2 (£ki)(ik 2 )[<p 2 (£k[ 12 \) - <p 2 {£k 2 )\ x 
4 b 2 Jo Jo Jo 


x e -| CT '- ff "|(fc 2 +fc 2 )_a'fc 2 2 ] - ( T"fef 12] _ 
Now perform the change of variables ex' —> £ 2 <x', a" —> e 2 a" to obtain 
4 = 4s 2 [ dk\dk 2 [ dcr 


' E 2 


'0 Jo 


/£ 2 ro-4 2 

dcr' / do-"(iefe[ 1 2])^ 2 (e/ei)(iefe2)[^ 2 (eA:[i2]) - <£ 2 (£fc 2 )]x 


'o 


x e 


— | cr / —cr" | (e 2 A: J+e 2 /c2)—(cr , U-cr /, )e 2 fcj 2 


[121 


By taking the limit £ —>• 0, the two sums over fci and k 2 become integrals: 

coo poo p 

lirn 4 = 4t / da' da" dkidk 2 (ik[ 12 ])(p 2 (ki)(ik 2 )[(p 2 (k[ 12 ]) - tp 2 (k 2 )]x 
l ~ .1 0 .1 0 J R 2 


x e 


— I <t'— a" I (fc J+fc|) — (<x' +cr") kf. 


[ 12 ] 


= 4t / dAqdA:2(%2])^ (ki)(ik 2 )[<p (k { 12 ]) - (p {k 2 )j 2 2 2 2 

JR 2 rC[i2] + K 2 “T K [12] 

r \ \ 

= 2 1 dk 1 dk 2 (ik [12] )[2cp 2 (k 1 )(p 2 (k [12] ){ik 2 ) - V 2 (ki)(ik [l2] )ip 2 {k^}-^— 

JR 2 ^[12] "'l ~r rv2 ' "'f] 

= 2t [ d/qd4 (ik [rzi ( ik[12] ) V k[12] ^ “ ^ 2 ^ 2 ^ + ^ 2 (fci)^ 2 (%2])*(^2 - h) 


[ 12 ] 


*1121 (k\ + k 2 + kf 12] ) 


= —21 / dk\dk 2 
J R 2 

which is indeed finite. 


[12] V^l w n-2 t fvj 12 p 

<^ 2 (^i)(^ 2 (fe[i2]) - ^(h)) + y 2 (fci)y 2 (fc [ i2]) ^ 2 fc[i ^ fc2 + fc 2 + fc 2 > 


10 Stochastic data for the Sasamoto-Spohn model 


10.1 Convergence of the RBE-enhancement 

Here we study the convergence of the data 


x N ,x%,xl,x},x%,Q N ,B N (Q N ox N ) 


for the discrete Burgers equation (74). We will pick up some correction terms as we pass 
to the limit, which is due to the fact that in the continuous setting not all kernels were 
absolutely integrable and at some points we used certain symmetries that are violated 
now. We work in the following setting: 
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Assumption (f,g,h). Let f,g,h € C£(M, C) satisfy /( 0) = g(0) = h(0, 0) = 0 and 
assume that f is a real valued even function with f(x) ^ Cf > 0 for all |.x|^ ir, that 
g(—x) = g(x)*, and that h(x,y ) = h(y,x ) and h(—x,—y) = h(x,y)*; we write h(x) = 
h(x, 0). Define for N € N 

J^Aa np(k) = -\k\ 2 f £ (k)^(p(k), JPD N ip(k) = ikg £ (k)&ip(k), 

&B N (<p,rl>)(k) = (27 -£)!*(£, (k-£)), 

e 

where £ = 2n/N and f e (x) = f{ex) and similarly for g and h. Let f be a space-time 
white noise and set 

Xjv(0 = (X n (0,Xn(0,xM),Xn(0,Xn(0,B n (Q n oX n )(0), 


where 

&x n {£) 

-2^(0 

^5(0 

-2^(0 

^Qiv(£) 


DjV^jvC, 

DjvnjvBjv^jv (0 ,Xjv(0), 

DjvnArBjv^J^oXjvCO), 
D 7V n 7 v J B 7 v(A)'(e),AV(0), 
Djv-Bjv(.Xjv(£), 1), 


(85) 


all with zero initial conditions except X/v(£) /or which we choose the “stationary” initial 
condition 

Xjv(0(0) = [ e- sf ^ 2 (D)D N V N f(s)ds. 

J —oo 


We are ready to state the main result of this section: 


Theorem 10.1. Let £ be a space-time white noise and make Assumption (f,g,h). Then 
for any 0 < 5 < T and p > q the sequence = Xtv(£) converges to 


X = (A, A v , A V + 2cQ, A V , + cQ v + 2cQ Q oX , QoX + c), 


in LP[Lt,C T ^ a - 1 x CW 2 "" 1 xif T “x J^ 2 " x J^ 2 " x C(\5,T\,K 2a ~ l )\ and is uniformly 
bounded in Cr^" -1 x Cr^ 2 " -1 x x x 2z?^ a x Cr*^ 2 " -1 ]. Lfere we wrote 


1 f n Im (g(x)h(x)) h(x,—x)\g(x)\ 2 

4tt J 0 x \f{x)\ 2 


and 

2z?Q v = DX V , ATQ QoX = D(QoX), 
both with 0 initial condition. 
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The proof of this theorem will occupy us for the remainder of the section. Define the 
kernel 

H t N (k) = ttzot ]kl<N/2 ikg £ (k)e- tk2 MV 

which satisfies 

r AT AT \ a (k)\ 2 e - fc2 /e(fc)l*-d 

/ r da = % |<w/2 ^y- - - (87) 

for all k G E. Let us start then by analyzing the resonant product Bn(Qn o Xjy). 


Lemma 10.2. For all 0 < 5 < T, the resonant products (Bn(Qn o Xn)) are bounded 
in L P (D, and converge to Q o X + c in L P (D, C'([<5, T], < if 2a ~ 1 )). 

Proof. We have the chaos decomposition 

B N (Q N oX N )(t,x) = / K N (t,x,7j 12 )W(dg 12 ) + (2-iry 1 Xjv(i,a;,T?i(-i))d? 7 i, 

l(IxE ) 2 JRxE 


where 

K N (t, x, 7/ 12 ) = [' d(reik[ 12 } x ijj 0 (/ci, k 2 )h £ (/ci, k 2 )H^_ a (k\)h £ [k\)H^_ s ^ [k\)H^_ S2 {k 2 f 

Jo 

Now Kj\f satisfies the same bounds as the kernel in the definition of Q o X and converges 
pointwise to it, so that the convergence of the second order Wiener-Ito integral over Kn 
to Qo X in L P (D, Ct^ 201-1 ) follows from the dominated convergence theorem. For the 
term in the chaos of order 0 we use (87) to obtain 


K N (t, x, ??!(-!))dr/i 


/RxE 


= J dh J dah e (k\, -fc 1 )g^ g (fc 1 )/i e (fc 1 ) ™ ---- 

r _ \n r (k-l )\ 2 e — /e(fci)(*—o-) 

= / d/ci / derh e (ki,-ki)t\ kl \ <N / 2 ikig £ (ki)h £ (ki y - --- 

Je Jo Je{Kl) z 


= e 


^ .g(£ki)h(eki) h(eki, —£k\)\g{ek \)\ 2 ^ - 2 k ?fBki)u ^ 0 o^ 
2^ e fci A.fVrkA 1 6 ’■ 


|fci|<JV/2 


4/ 2 (eA:i) 


It is not hard to see that the sum involving the exponential correction term converges 
to zero uniformly in t € [<5, T] whenever 5 > 0. For the remaining sum, note that 

, g(-eki)h(-eki) _ . f gjek^hfeh) \* 

— ek\ \ ek\ / ’ 


while the second fraction on the right hand side of (88) is an even function of k\. and 
thus 


£ 


E 

\ki\<N/2 


.g(eki)h(£ki) h(eki, -£ki)\g(£ki)\ 2 
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— £ 

0<fci<iV/2 


lm(g(eki)h(eki)) h(ski, -ski)\g(eki)\^ 


ski 


2 P(eh) 


Now lm(gh) € Cl with lm(g^)(0) = 0, and therefore lm(g(x)h(x))/x is a bounded 
continuous function. Since furthermore f{sk\) Cf > 0 for all |fci|< N/ 2, the right 
hand side is a Riemann sum and converges to 2irc which is what we wanted to show. □ 


Next, define 
and the vertex function 


E n = E n (-N/2, N/2) 


( \ n (U n \|2 p— k2fe(k2)\cr\ 

VZ(a,ki) = l, fcl | <JV/2 / ^ £ (A;f 2] ,A;_ 2 )/ l£ (A ;i ,fc 2 )^(A;f 2] )^y4- tt -dfc 2 . 




(89) 


We also set 




V N (a,ki) = t lkl \^ N/2 / dfc 2 ^o(fc [12] , /c_ 2 )/i £ (/c [12] , fc_ 2 ) 

J E n 

Ar 1(7 f/to')I 2 e — fc i/e( fc 2)|o-| 

x / l£ (fc 1 ,fc 2 )^(^ 2] )M^-— -. 

Lemma 10.3. Let cGRk the constant defined in (86). Then V^(-, k ) converges weakly 
to V^(a, k) + 2irc5(cr) for all k € E, in the sense that 

roo poo 

lim / ip(a)V^(a,k)da = / (p(a)V^(a,k)da + 27rap(0) 

N^ooJq J o 


R with \ip{cr) — </>(0)|< \a\ K for some k > 0. Similarly, 

V) yn 

Vy(-. k) converges weakly to V°(a,k) + 2ecS(o) for all k € E. Moreover, for all 5 > 0 


for all measurable ip: [0, oo) 

*#(•,*) 
we have 


POO /-OO 

sup / |V)v((j,/c)|dcr + sup / |V)v(cj, t)|dcr 

JV Jo TV Jo 


< 

r-j 


Proof We write V$(a,k) = {V^(a,k) - t\ k ^ N/2 V%{a, 0)) + °)- Let us 

hrst concentrate on the second term: 


rVt 




PQt 


datp(a)V^(o, 0) 

= f dk 2 h £ (k 2 , ~k 2 )h £ (k 2 )(ik 2 )g £ (k 2 ) [ 

Je n 2,j £ [k 2 ) J o 


(90) 


daip{a)e- 2k ^ k ^ c 


Y h e (k 2 , -k 2 )k 2 lm(h £ (k 2 )g £ (k 2 )) 
0<\k 2 \^N/2 


\de{k 2 )\ 

fe{k 2 ) 


2 r °o 


daf(a)e- 2k i Mk2)rT . 


Now add and subtract <p(0) in the integral over a and observe that 

POO POO POO 

/ da(ip(a)-f(0))e- 2k2 ^ k ^ < daa K e~ 2k ^ < \k 2 \~ 2K / dae~ ck ^ <\k 2 \ 

Jo Jo Jo 


—2—2 Hi 
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and then 


0 <\k 2 \^N/2 


<S K ' 


^2 M fe 2 , -fc 2 )fc 2 Im(fe e (fc 2 )ff e (fc 2 ))IM 2 2k 

^ Je{k 2 ) 

' (e ^2 \H £k 2,-£k2)lm(h(ek2)g(£k 2 ))\ 

r> I _ I ^ tv r /o J \ 2,) 


0<\k 2 \^N/2 


whenever n' G [0,1 A 2k). The term in the brackets is a Riemann sum and since 
\Ini(hg)(x)\< |x| and k' < 1, it converges to a finite limit. Thus we may replace (p(a ) 
by <£>(0) and end up with 


-¥>(0) ^2 h £ (k 2 ,-k 2 )lm(h £ (k 2 )g £ (k 2 )) 


0 <\k 2 \^N/2 


\9e(k 2 )\ 2 

VKk2)k 2 ' 


which converges to 27tc</?(0). 

It remains to treat the term 


V^{a,k) - t\ k ^ N/2 V^(a,0) = t\ k \^ N / 2 I dk 2 (W^(a,k,k 2 ) - W%(a,0,k 2 )), 

J e n 

where the right hand side is to be understood as the (indirect) definition of W^. Since 
for fixed ( k , k 2 ,a) the integrand converges to (H a (k+k 2 ) — H a (k 2 ))e~ k ^ /2, it suffices to 
bound \W%(a, k , k 2 ) — W^(a, 0, k 2 )| uniformly in N by an expression which is integrable 
over (a,k 2 ) G [0,oo) x E. For the remainder of the proof let us write i/j £ (£,m) = 
h £ ((£ + m ) N , — m)h s {£, m)g £ ((£ + m) N ), which satisfies uniformly over \£\, \m\< N/2 

\ip £ (£, m) - ip £ ( 0, m )|< e\£\< (e\£\) 5 

for all 5 G [0,1] (to bound \{£ + m) N — m\ note that (£ + m) N = £ + m + j(m,£)N for 
some j(m,£) G Z and that if \£\< N/2 , then \£ + jN\£\ for all j G Z). We now have 
to estimate the term 


\'ip e (k, k 2 )(k + k 2 ) N e~ af ^ k+k ^ N ^ k+k ^ N ^ - Vg( 0, C fe2 >|i| u | fc2 1 < e -)o 

^ |^ c (fe, fc 2 ) - ^ e (0, fe 2 )||^|e-^«C fc2 ^l| fcU | fca|<JV/2 e- fc 5 c ^ 

+ fc 2 )||(fc + fc 2 ) Jv e —^« fe+fc2 ) JV )« fe - < - fc2 ) JV ) 2 - fc 2 e—^^ 2 > fe i|l| fc |,| fc 2 | <JV/ 2 e- fc ^/- 

\N\ 


\JV)2 

\K\2 


<((e\k\) 5 \k 2 \+\(k + k 2 ) N e-^ k + k ^ )« fc+fc2 ) > -fe 2 e-^( fc2 ) fc 2 |)l| fcUfc2 | <7V/2 e- fc 2 c / c 


The first addend on the right hand side is bounded by \k\ 8 \k 2 \ 1 ~ 8 e~ crc f k 2 and for 6 > 0 
this is integrable in (cr,k 2 ) and the integral is < \k\ 8 . To bound the second addend, let 
us define ip £ ^(x) = X e~ a f e ^ x and note that for |x|< N/2 


WeA x )\= I 1 - xcrf (ex)£x 2 


xaf £ (x)2x\e~ af ^ x)x2 < (1 + a\x\ 2 )e~ c f° x2 


< e ~ acx2 
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for some c > 0. From here we can use the same arguments as in the proof of Lemma 9.5 
to show that 


_|_ fo'jN e -'rfe({k+k2) N )((k+k 2 ) N ) 2 


\ m<N/2 e~ k *cf° < F(k,k 2 ,a) 


for a function F with f E dk 2 / 0 °° daF(k, k 2 , a) < Ifcl* 5 , and thus we get the convergence 

yo V7 

and the bound for Vff. The term can be treated using the same arguments. □ 

Proof of Theorem 10.1. We introduce analogous kernels as in the continuous setting: 
Define G 9 N (t,x,r)) = e lkx H^f_ s {k) and and then inductively for r = (t\t 2 ) 


• i JV 

G T N (t,x,r ) T ) = e H x H T N (t,ii T ) 

= e ik M x £ daH t N _ a (k^ T] )h £ (kf Ti] , k^)H^ (a, r, T2 ). 

A first consequence of this recursive description is that the contribution to the 0-th chaos 
always vanishes: just as for the kernels G T of Section 9 we have 


G T N(t,x,a(ri 1 ...n(-i)...(-n))) = 0 

whenever a G S 2n , because G T N {t,x,r]\... 2n ) oc [k 2 2n \. 

We decompose every G T N into two parts: 


G T N (t,x,r ] T ) 


G Er {t,X,r]T)t\k 1 \,...,\k d(T) \^N/(2d(T)) +G T N (t,X,7 j T )(l - l|fc 1 |,...,|fc d(T) |^V/(2d(T))) 
K T N {t,x,r) T ) + Fff(t,x,r] T ). (91) 


Let us Hrst indicate how to deal with which gives a vanishing contribution in the 
limit. Define 

„ _ ii ikg e {k) 

Qn( 0) - hkHN/2 ico + fAk)k2 , 

which satisfies |(?at( 0)|< |#| _1 , uniformly in N (recall that we defined |0|= |w| 1/,2 -|-|A;|). 
Similarly as in Section 9.5 we can bound every integrand that we need to control by 
a product of terms of the form </jv(0j^,i), discrete vertex functions, and factors like 
ipo(kr^n, fcjy,]), where d N = (u,k N ). Moreover, every integrand contains a factor \qN{@i)\ 
for each its integration variables 0i, and we can decompose 1 — H|| \k d(T) \^N/(2d(r)) into 

a finite sum of terms that each contain a factor ^-\ 0 i \>N/{ 2 d{r)) for some i. We can there¬ 
fore estimate |? 7 v(^:)|l| 6 » i |>iv/( 2 d(r)) 1$ N~ 5 \0i\ 1 ^ s for an arbitrarily small 5 > 0, which 
gives us a small factor. We now only have to show that the remaining integral stays 
bounded. To do so, we need to control the basic integral 


'Rx En 


\0\- a \{0' - 0) N \~Pd01 


\k'\^N/2- 
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But (kJ — k) N = k' — k + jN for some |j|^ 1 and if N/2, then \k' + jN[^ \k'\ for 
all j € Z. So by Lemma 9.8 we have 

[ \o\-«\{e' -e) N \-P&ot! \o\- a \o' + (o,jN)-e\-Pdei lk ,^ N/2 

«/R xEn |j|<1 ^ Bn 

< ^2 P l|fc'|<Af/2 < 1^1 P l|fc'|^iV/2! 

\M 1 

where a, f3, p are as in the announcement of the lemma. The second integral in Lemma 9.8 
is estimated using the same argument, which also allows us to show that all discrete 
vertex functions apart from and V° T satisfy the same bounds as their continuous 

Y) . 

counterparts. Since we estimated Vjy and in Lemma 10.3, we now simply need 
to repeat the calculations of Section 9.5 to show that the Wiener-Ito integral over FJj 
converges to 0 in L p in the appropriate Besov space. 

We still have to treat the term KJ^ in the decomposition (91) of G T N . In the descrip¬ 
tion of I\]y we can now replace every kF, j by the usual sum k^ (where t' is an arbitrary 
subtree of r). Then K] y satisfies the same bounds as G T , uniformly in N, and converges 
pointwise to it. So whenever G T is absolutely integrable we can apply the dominated 
convergence theorem to conclude. However, in the continuous case we used certain sym- 
metries to derive the bounds for QoX, V v , and V°, and these symmetries turn out to 
be violated in the discrete case. This is why we separately studied the convergence of 
Bn{Qn ° -Xjv) in Lemma 10.2, and in Lemma 10.3 we showed that and satisfy 

the same bounds as V ^ and V°, uniformly in N, and converge to V ^ + 27rce>(cr) and 
Y? 

V° + 2-7rcS(a) respectively. It thus remains to see which correction terms we pick up 
from the additional Dirac deltas. 

Let us start with the contraction of GJj. Here we have 


X N(t,x)=[ G%(t,x,r}\2z)W(&r]i2z) + [ G\ 1 (t,x,r]i)W(drii) 

J (RxE ) 3 JRxE 


with 

W 


g n, i(t,x,rn) = tt 1 Gjf (t, x, T]i) = tt 1 e lkix I da I da'H^k^H^^k^V^(a-a' , h) 

Jo Jo 

for as defined in (89). Now by Lemma 10.3, the right hand side converges to 

ir~ 1 G^(t,x,r]i) + 2ce ikix f d*H t - <T (k 1 )H a - ai (h), 

Jo 

and we have 


rVt 


IxE 


2ce iklX I daH t _ a {k 1 )H a _ sl (k 1 ))W(dr ll ) = 2cQ(t,x). 
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Y Y 

In conclusion, X^ converges to X v + 2cQ in L p (12,). 

The only other place where appears is in the contractions of X£. We have 




-A 


.'fe 


Xfi(t,x)= / G^(t,x,77i 234 )Vh(dr/i234) + / G$ j2 (t, x, 7/i 2 )VF(d7yi 2 ) 


5xE) 4 


xE) 2 


with 


g1 2 (M,»1i2) = (2tt) 1 (G^(t,x,r/i 2 ) + 2G^(t,x,??i 2 ) + 2G'|(t,x,?7i 2 )). 


Now Gjy(t,x,r) 12 ) can be factorized as 

rt ra pa 


Gl(t,x, Vl2 ) = e ik M x h £ (k 1 ,k 2 ) [ da f da' f da' 1 H t N _ a (k [l2] )H^_ a „(k [12] ) 

Jo Jo Jo 

x H%,_ 81 (k{)H$,_ a2 (k 2 )vg{a - a',k [12] ), 


and Lemma 10.3 shows that the right hand side converges to 

rt r<* 


G v (t,x,ri 12 ) + 2irce lk { 12 l x 


Je 

7 

with 


/ da / da Hi— a {k^i 2 ^)H a — (T ii(k^\ 2 ^)H a "— Sl {k\)H (T ii— S2 (k 2 ) ) 

r o Jo 


IxE) 2 


(e ik M x J da jT d o"H t _ a (k [12] )H a _ a „ ( k [12] )H a "_ sl (fci)iL CT »_ S2 (fc 2 )) W(d m 
= Q v (t,x). 

yj 

Similarly, we factorize G^(t,x,rf i 2 ) as 


G%(t,x, m2 ) = e*M x h e (k 1 ,k 2 )MkiM f da [da 1 [ da"H t N _ ff (k [12] )H?_ 82 (k 2 ) 

Jo Jo Jo 

X ( J „n 


x H?_ ff ,(ki)HZ,_ ai (k 1 )Vj(<T' - a", ki), 
and Lemma 10.3 shows that the right hand side converges to 

G^{t,x,r] 12 )+2irce lk l 12 y x 'iJj 0 (k 1 ,k 2 ) [ da ( da'H t _ <J {k^ 2 ])H fT ^ S2 {k 2 )H (T _ a r(k{)H a r_ sl {k 1 ) 

Jo Jo 


with 


IxE) 2 


e ik M x MkiM [da [ da , H t - a (k [12] )H a - aa (k 2 )H < ,- tT ,(k 1 )H IT ,-8 1 

Jo Jo 


= Q QoX (t,x). 

In conclusion, converges to X^ + cQ v + 2cQ^ oA " in L p (fl,2z?^ Q ). □ 
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10.2 Convergence of the random operator 

The purpose of this subsection is to prove that the random operator 


AnU) :=nAr( J(H N (f -< T- £y Q N ) o ( T - £Z X N )n(dy, 

-V N (J ({f^Q N ))o(T- e ,X N Mdy,dzj) 

converges to zero in probability. 

Theorem 10.4. Let a € (1/4,1/2). Then we have for all 5, T > 0 and r ^ 1 

n\M\c T L^-^' r s jv«- i/2+< . 

We will prove the theorem by building on three auxiliary lemmas. 
Lemma 10.5. The operator Ajy is given by 

A N (f)(t,x ) = Y A q A N (A p f)(t,x) = Y j 9p, q (t,x,y)A p f(y)dy 

q,p q,p JT 


with 


where 


^g pq {t,x,-){k) = Y ^p tq {x;k,ki,k 2 )^ r QN(t,ki)^ r X N (t,k 2 ), 

ki,kz 


( 92 ) 


(93) 


(94) 


Fp q (x;k,ki,k 2 ) = (2tt) 2 p p {k)if^{k, ki)h £ (ki, k 2 ) 

x [e* k M- k)Nx Pq {{k [12] -k) N )M(h ~k) N ,k 2 ) 

- e l(k[12] ~ k)x Pq(k [12] - k)-if 0 {ki - k,k 2 )t lk[12] _ k ^ N/2 \ 

and where p p is a smooth function supported in an annulus 2 V stf such that p p p p = p p . 
Proof. Parseval’s formula gives 

A q A N (A p f)(t,x) = [ g p (t,x,y)A p f(y)dy 
J T 

= (2vr) _1 Y ^ 9 p, q {t, x , -){-k)p p {k)^f{k). 

k 

It suffices to verify that the same identity holds if we replace & g p q {t,x , -){k) with the 
right hand side of (94) and if we fix some k and consider f(x ) = (27r)~ 1 e zkx , i.e. & f(£) = 
Sf k . Let us look at the first term on the right hand side of (92) 

J A,n A r[(n J v(A p ((27r)- 1 e ife -)^r_ £ , 1 ((27r)- 1 e ifel -)))o(r_ £Z2 ((27r)- 1 e ife2 -))]/x(dzi,dz 2 ) 
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= (2 7 r)- 3 A g n 7V [(n i v(p p (A;)^(A:,fc 1 )e^ +fel )-))o(e^-)] J e ik ^ z ' +ik2£Z2 f i(dz 1 ,dz 2 ) 

= (2Tr)~ 3 A q U N [M(k + h) N , k 2 )p p (k)^(k, kJeW+W+^hefakz) 

= {2-k)~ 3 e l ^ k ^ +k)N 'p q {{k [l2 ] + k) N )if 0 ((k + ki) N , k 2 )p p (k)p p (k)^(k, ki)h e (ki, k 2 ), 

where in the last step we used that ((& + k\) N + k 2 ) N = (k + k\+k 2 ) N . The second term 
in (92) is treated with the same arguments, and writing Qn and Xn as inverse Fourier 
transforms equation (94) follows. □ 

Lemma 10.6. For all r ^ 1 and all a, f3 € R we can estimate 


n\\ A N (t) - A N {s)\\ r L{<)fct ^j 


<E 2 : ^ 2 - P -( sup ^ E[ | ( ^ (t , X) . ) _^ (S)X) . ))(A;) | 2 ] ) r/2 . 

q,p X&T k 

Proof. From the decomposition in Lemma 10.5 we get 

\\A N (f){t) - A N (f)(s)\\ r B p = ^2^11 Y f (g£ q (t,x,y)-g£ q (s,x,y))A p f(y)dy 


< 


Y 2 qrf) 2 


v 

—pra 


r 

£5( T) 




\ 9 p,q(t, x, y) - 9 p >q {s, X , y)\dy 


N 


Lli T) 


so that 


E[||Ajv(t) - Ajv(s) \\ r LiVa > B P r) 

< Y 2 qr ^2~ pra K 

q,p 

< Y 2 qr ^2~ pra K 

q,p 


1 9p, q (t,x,y) - g£ q (s,x,y)\dy 


p,q\ 


L S (T)J 


1 9p, q (t,x,y) - g£ q (s,x,y)\ 2 dy) 


S 1/2 

r 

) 

m t). 


Y 2qrP 2~ pra I E [(El^*( i ’ a; ’-)-^&(s,*,-))W| 2 ) r/2 ]dx, 


q,p 


where we applied Parseval’s formula. Now Ylk\(^9p,q(ti x ->') — ^9p,q{ s i x i '))(^)| 2 a 
random variable in the second inhomogeneous chaos generated by £, and therefore our 
claim follows from Nelson’s estimate. □ 


Lemma 10.7. For all p, q ^ — 1, all 0 ^ t\ < t 2 , and all A G [0,1] we have 
Ym^ r gp,q( t2 ’X,-)(k)-^ r g P ,q(ti,x,-)(k)\ 2 ] < l 2 p :2 q< N 2 pl ' 1 ~ x ^ +q N~ 1+3X \t 2 — ti\ x . (95) 

k 
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Proof. We write 

^Q N (t 2 , ki)^X N (t 2 , k 2 ) - ^ r QN(ti,ki)^ r X N (ti,k 2 ) 

= (&Q N (t 2 , h) - &Q N (h, ki))&X N (t 2 , k 2 ) 

+ &QN(h,ki)(&X N (t 2 , k 2 ) - &X N (ti,k 2 )) 

and estimate both terms on the right hand side separately. We only concentrate on 
the second one, which is slightly more difficult to treat, the first one being accessible to 
essentially the same arguments. We have the following chaos decomposition: 


&QN(ti,ki)(&X N (t 2 , k 2 ) - &X N (ti,k 2 )) 

[ (V ) 2 5 ilM Se 2M r <_ r (^i)^ Sl (^i)dr«_ S2 (^ 2 ) - Hi[_ aa (e 2 ))) W(d Vl2 ) 

i(BxE w ) 2 V JO J 

'stukj-tuk* f 1 - H*[_ n (-e r))) d m . 


lxE N p 

+ 2vr [ 

J Ex E n 

Using the same arguments as in Section 9.5, we can estimate the second term on the 
right hand side by 


txE]\[ 


j ' - <_„(-«i))l dm 




for any A € [0,2]. Based on the decomposition of Lemma 10.5, the contribution from 
the chaos of order 0 is thus bounded by 


y, | y: Tp !q (x\k,ki,k 2 )E[&Q N (ti,k{)(&X N (t 2 ,k 2 )-&X N {t 1 ,k 2 ))\ 

k k!,k 2 €E N 


(96) 


I rN („. u u u \ l^ 2 — ^1 ~ 2 

k ki£E N 1 11 

~ X] \pq( k )Pp(. k ) ^2 1 p-<(k,k l)(^o((fcl - k) N ,h) - ipoih -fc,fci)) ^ 2 2\ x , 

k fciSEjv 1 


where we used that |fc|< N/2 on the support of Now r f 0 ((ki — k) N ,k\) = 

i(j 0 ( k i — k,ki) unless \k\ — k\> N/2 and |fci|~ N, and there are at most \k\ values of 
h with |L|< |fei|< N/2 and |fci — k\> N/2. Therefore, the right hand side of (96) is 
bounded by 


< 1 1 2 - h\ x J2 hk\<N/2\p q {k)p p {k)\ 2 \k\ 2 N- 2+2X < t p ~ g t 2q < N \t2 - t\\ x 2 3q N~ 2+2X 

k 

< t 2 P t2 q< N 2 P ( 1 ~ X ' >+g N~ 1+3X \t 2 - h\ X . 
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Next consider the contribution from the second chaos, which is given by 


(2vr) 2 [ (T^ q {x-kMM) [ L H^MH^k^dr 
J(RxE ) 2 V Jo 

x (H t N 2 _ S2 (k 2 ) - 

Taking the expectation of the norm squared, we obtain (up to a multiple of 27 t) 

[ r£,(z; k, h ,k 2 ) r - <_ S 2 (/c 2 )) 

J(RxEjv ) 2 Jo 


d?7i2 - 
(97) 


As in Section 9.5 we can show that 
ft i 


/o 


<-r(fe 1 )^ ai (fc 1 )dr«_ i 2 (fc 2 ) - <- S 2 (fc 2 )) 


dsi 2 < 


|fc 2 | 2 A |t 2 -ti| A 

M 2 


whenever A € [0,2], and plugging this into (97) we get 


(97) 


< 


T^(x-,k,k 1 ,k 2 )\ 


’E% 


2 |fc 2 | 2 A |t 2 -ti|' 


dfei 2 


<N 2 *\t 2 - tl \' 


>E 2 n 


[fell 2 

dki 2 pl(k)i>%(k, A;i)|A:ir 2 

e m^-V N * Pq ((k ll2] - k) N )M(h - k) N ,k 2 ) 

- e i(fc l 12 l- fc) >,(fe [12] - fctyo(fci - k, k 2 )t lk[12] _ k] ^ N/2 


First observe that the difference on the right hand side is zero unless |fei|~ N, so that 
we can estimate |&i|~ 2 < fV - 1 +A |/c| _1 ~ A , and also we only have to sum over |fci|> \k\. 
Moreover, the summation over k 2 gives 0(2 q ) terms which leads to 

(97) < l 2 P , 29 < A r 2 < ?N- 1 + 3 A |t 2 - t 1 | A p 2 (fe)|fcr A , 

and the sum over k of the right hand side is bounded by 

< l 2 p, 2 9<Ar2 p ^ 1-A ) +<? N -1+3A |f 2 - ti| A , 


which completes the proof of (95). 


□ 


Proof of Theorem 10. f. Let a € (1/4,1/2), r ^ 1, and write 7 = 2a—1+1 /r. Combining 
Lemma 10.6 and Lemma 10.7, we get for all A ^ 0 with (1 — A)/2 > a 

E|M*(f) - < E[||. 4 N (/)(t) - A N (f)(s) ir i( ^, B ;, | ] 1/r 

^ 2 9 ’' 7 2 -f ’ r “(2 p * 1-A * +9 iV -1+3A |( — s| A 7 2 + r 

Q,pS°S2 N 
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~ s |A/2 jya—l/2+l/r+3A/2 


Since .Ajv(O) = 0, Kolmogorov’s continuity criterion gives 

nuNm r c T L^,^)\ llr 5. v“- 1/2+1 / r+ "/ 2 

whenever A/2 > 1/r. Choosing 1/r + 3A/2 ^ 5, the estimate (93) follows. 


□ 


A Proofs of some auxiliary results 

Proof of Lemma 6.f. We prove the second claim, the first one follows from the same 
arguments. We also assume that / € M^L P , the case / € MfL°° is again completely 
analogous. Let t € [0, T] and j ^ 0 and write M = II/IIju^lpIIpWII/?- Then 

f 2 2 i<p(2 2j (i — s))S/_i/(s V 0)dsAjg(t) ^ 2~^M f 2 2j \p(2 2 ^ (t — s))|s _ 7 d-s. 
o LP Jo 

We now split the integral at i/2. The Schwartz function p satisfies |</?(r)|< r -1 , so using 
the fact that 76 [ 0 , 1 ) we get 


/•i/2 rt/Z r-t/Z 

/ 2 2j |^(2 2 j (i-s))|s _7 ds< / (t-s)~ 1 s~' Y ds<t^ 1 s _7 ds < i -7 . 

fo Lo j 0 

The remaining part of the integral can be simply estimated by 

f 2 2j |< j c( 2 2j (t — s))|s _7 ds ^ (i/2) -7 / 2 2 >(2^(i - s))|d* < r 7 , 

Jt/2 Jr 

which concludes the proof. 


rt/2 


i/2 


□ 


Proof of Lemma 6.5. Let us begin with the first claim. By spectral support properties 
it suffices to control 


(J 2 2 V(2 2j (i - s))Sj-if (s)ds - 5,_i/(i))A, 5 (i) 


LP 


(98) 


< 


2 2 V(2 2j (i - s))Sj-if (t)dsAjg(t) 


lp 


+ 


2 2] p(2 2 \t - s))Sj-ift, s dsAjg(t) 


Lp 


where we used that tp has total mass 1. Using that \p(r )|< |?’| 1 a / 2 , we can estimate 
the first term on the right hand side by 

roo 

< 2 -J/> r 7+o/2ll/ll || 9 ( t )|| (i / | r |-W2 d r < 2 -««+P) t -T|l / || || s ,( i )|| (3 . 

J 2 2 ii 
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As for the second term in (98), we split the domain of integration into the intervals [0, i/2] 
and [t/2,t]. On the first interval we use again that |</?(r)|< | r |-i -<*/2 and then simply 
estimate ||/t, s ||LP^ (t 0 / 2 ^ + s a / 2 ~ 7 )||/||and (t — s) -1-a / 2 < i _1_a / 2 , and the 
required bound follows. On the second interval, an application of the triangle inequality 
shows that 


\\ft,shp= ll« 7 (« 7 /(s))-i 7 (* 7 /(*))l|i><* 7 |f — s| Q/ 2 ||/|li?j' c 


(O’ 


from where we easily deduce the claimed estimate. 

Let us now get to the bound f || ^ g) ~ f ^ (Jfg))(t)\\^+p -2 < ||/||^'.“(t)|b(0ll/3- 


Given Lemma 6.4, only the estimate for 

rt 


dt (J ftvm — s))Sj-iu(s)ds S j A jv(t). 

is non-trivial. For fixed j we obtain 

f 2^ip'(2 2 ^[t — s))Sj-\u{s)ds = 2 2 ^ f 2 2 ^(p'(2 2 ^(t — s))Sj-iu(s)t s ^ods. 

Jo Jr 

Since pJ integrates to zero we can subtract 2 2 - 7 f R 2 2j p'(2 2 ^ (t — s))Sj-\u(t)ds. The result 
then follows as in the first part of the proof. □ 

Proof of Lemma 6.6. We start by observing that exactly the same arguments as in 
Lemma 2.9 (replacing L°° by L p at the appropriate places) yield 

ll^/ll^f p 00 ’“(T)^ II/IIct^p'~ 2 ’ II s ^ INIk«- 

It remains to include the possible blow up at 0. The estimate (49), 


\m 




MfVp 2 


(99) 


is shown in Lemma A.9 of [GIP15] (for p = oo, but the extension to general p is again 
trivial). For the temporal Holder regularity of If that we need in (48), we note that the 
estimate 

for /3 ^ —a is again a simple extension from p = oo to general p, this time of Lemma A.7 
in [GIP15]. An interpolation argument then yields 


II PM\lp< t P/2 \\u\L-p 


(100) 


whenever (3 < 0, and from here we obtain 

rt 


|| (P f - id)it||z, P = f d s P s uds = f 

Jo ^ J 0 


PcAuds 


lp 
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< »- 1 + “ / 2 ll«lk»d» < i “ /2 IMI*. ( 101 ) 

for a € (0, 2). To estimate the temporal regularity, we now have to control 
P [ Pt-rfAr -s' 1 [ Pt-rfAr 

J 0 J 0 LP 

[ \\Pt-rfr\\LAr + s' 1 [ WPt-rfrhAr + s' 1 (P t . s - id)/ f(s) 

Jo Js 


< (^(P — s' 1 ) f ( t — ?’) Q//2 V 7 dr + s 7 f (t — r )"/ 2 V 7 dr^ 

*/ 0 J s 


Al Cj?OL-2 


where we used (100) for the first two terms and (101) for the last term. Now (99) allows 
us to further simplify this to 

< ((C - »’)C /2 -T /"(l- r )«/2-l r -7 dr+ /'‘( f - r )« /2 -l dr+ |« _ S |»/ 2 ) 

Since a > 0 and 7 < 1, the first time integral is finite. Moreover it is not hard to see that 
(P — s 7 )t Q / 2-7 < \t — s | Q / 2 (distinguish for example the cases s ^ i/2 and s € (i/2,i]), 
and eventually we obtain 


||C//(t) - iV/Mlli-S 1 1 - »p /2 ii/ii aw -». 

It remains to control the temporal regularity of s 1 —)• s 7 P s uo, but this can be done using 
similar (but simpler) arguments as above, so that the proof is complete. □ 

Proof of Lemma 6.8. For the spatial regularity, observe that 

t 1 11 Aj u (t) 11 lp ^ min{ 2~ : > a , t a/2 } 11 u \ \ #g> a ( T ) • 

Interpolating, we obtain 


or in other words i 7 _£ / 2 ||u(i)|| a _ £ ^ IMIjs?j 7 ’ a (T)- The statement about the temporal 
regularity is a special case of the following lemma. □ 

Lemma A.l. Let a € (0,1), e € [0, a), and let f: [0, 00 ) X be an a-Holder continuous 
function with values in a normed vector space X, such that /(0) = 0. Then 

\\t^t- £ f{t)\\ c «-e x < \\f\\c«x- 
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Proof. Let 0 ^ s < t. If s = 0, the required estimate easily follows from the fact that 
/(0) = 0, so assume s > 0. If t > 2s, then we use again that /(0) = 0 to obtain 


£ f(s)\\x 




Now t a £ ^ (|f — s|+s)“ £ , and s/\t — s|< 1 by assumption, thus the result follows. If 
s < t ^ 2s, we apply a first order Taylor expansion to t~ £ — s~ £ and obtain 

ll^~ £ /W - s~ £ f(s)\\x < t~ £ \\f{t) - f(s)\\x £r~ £ ^ 1 (t-s)||/(s)||,Y 

\t-s\ a ~ £ ^ |f-s|" £ |f-s| a \t-s\ a ~ £ 

for some r € (s,t). Now clearly the hrst term on the right hand side is bounded by 
\\f\\c<*x- For the second term we use ||/(s)||x< s a \\f\\ C zx and get 


er £ 1 (t-s)||/(s)|| A ' 


,-£-l+au „|£+1 —a 


1 1 — s|’ 




using t — s ^ s in the last step. 
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